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Abstract 

This thesis is an attempt to enhance understanding of the following 
questions 

A- Given a multipartite quantum state (possibly mixed), how to find out 
whether it is entangled or separable? (Detection of entanglement.) 
B- Given an entangled state, how to decide how much entangled it is? 
(Measure of entanglement.), in the context of multipartite quantum states. 

We have explored two approaches. In the first approach, we assign a 
weighted graph with multipartite quantum state and address the question 
of separability in terms of these graphs and various operations involving 
them. In the second approach we use the so called Bloch representation 
of multipartite quantum states to establish new criteria for detection of 
multipartite entangled states. We further give a new measure for entan- 
glement in A-qubit entangled pure state and formally extend it to cover 
A-qubit mixed states. 

We give a method to associate a graph with an arbitrary density matrix 
referred to a standard orthonormal basis in the Hilbert space of a finite di- 
mensional quantum system. We study related issues such as classification 
of pure and mixed states. Von Neumann entropy, separability of multi- 
partite quantum states and quantum operations in terms of the graphs 
associated with quantum states. In order to address the separability and 
entanglement questions using graphs, we introduce a modified tensor prod- 
uct of weighted graphs, and establish its algebraic properties. In particular, 
we show that Werner's definition (Werner 1989 Phys. Rev. A 40 4277) of 
a separable state can be written in terms of graphs, for the states in a real 
or complex Hilbert space. We generalize the separability criterion (degree 
criterion) due to Braunstein et al. (2006 Phys. Rev. A 73 012320) to a 
class of weighted graphs with real weights. We have given some criteria for 
the Laplacian associated with a weighted graph to be positive semidefinite. 
We settle the so-called degree conjecture for the separability of multipar- 
tite quantum states, which are normalized graph Laplacians, first given by 
Braunstein et al. [Phys. Rev. A 73, 012320 (2006)]. The conjecture states 
that a multipartite quantum state is separable if and only if the degree 
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matrix of the graph associated with the state is equal to the degree matrix 
of the partial transpose of this graph. We call this statement to be the 
strong form of the conjecture. In its weak version, the conjecture requires 
only the necessity, that is, if the state is separable, the corresponding de- 
gree matrices match. We prove the strong form of the conjecture for pure 
multipartite quantum states, using the modified tensor product of graphs 
defined by Ah S. M. Hassan and P. S. Joag [J. Phys. A 40, 10251 (2007)], 
as both necessary and sufficient condition for separability. Based on this 
proof, we give a polynomial-time algorithm for completely factorizing any 
pure multipartite quantum state. By polynomial-time algorithm we mean 
that the execution time of this algorithm increases as a polynomial in m, 
where m is the number of parts of the quantum system. We give a counter- 
example to show that the conjecture fails, in general, even in its weak form, 
for multipartite mixed states. Finally, we prove this conjecture, in its weak 
form, for a class of multipartite mixed states, giving only a necessary con- 
dition for separability. 

We give a new separability criterion, a necessary condition for separa- 
bility of A^-partite quantum states. The criterion is based on the Bloch 
representation of a A-partite quantum state and makes use of multilinear 
algebra, in particular, the matrization of tensors. Our criterion applies to 
arbitrary A-partite quantum states in 7-^ = 7i^^ TC'^'^ • • • ® 7Y^". The 
criterion can test whether a A-partite state is entangled and can be applied 
to different partitions of the A-partite system. We provide examples that 
show the ability of this criterion to detect entanglement. We show that this 
criterion can detect bound entangled states. We prove a sufficiency condi- 
tion for separability of a 3-partite state, straightforwardly generalizable to 
the case A > 3, under certain condition. We also give a necessary and suf- 
ficient condition for separability of a class of A-qubit states which includes 
A-qubit PPT states. We present a multipartite entanglement measure for 
A-qubit pure states, using the norm of the correlation tensor which occurs 
in the Bloch representation of the state. We compute this measure for 
several important classes of A-qubit pure states such as GHZ states, W 
states and their superpositions. We compute this measure for interesting 
applications like one dimensional Heisenberg antiferromagnet. We use this 
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measure to follow the entanglement dynamics of Grover's algorithm. We 
prove that this measure possesses almost all the properties expected of a 
good entanglement measure, including monotonicity. Finally, we extend 
this measure to A^-qubit mixed states via convex roof construction and 
establish its various properties, including its monotonicity. We also intro- 
duce a related measure which has all properties of the above measure and 
is also additive. 



Chapter 1 

Introduction and Overview 



Charak Samhita (First book of Aurveda) 
( Whole is not known by knowing its parts) 



■ ■ ■ But I CAN SAFELY SAY THAT NOBODY UNDERSTANDS QUANTUM MECHANICS- ■ ■ 

Richard Feynmann 

Entanglement is a subtle and eluding property of quantum systems com- 
prising many parts. Entanglement induces correlations between the mea- 
surable properties of different parts of a quantum system which cannot 
be reproduced by any procedure involving only the local operations (LO) 
on and classical communication (CC) between various parts of the system 
11]. In consonance with this, entanglement in a quantum system cannot 
increase (or be created) via LOCC. This principle is connected to another 
intriguing property of entanglement: a multipartite quantum system can 
get entangled in various inequivalent ways, which cannot be transformed 
into each other via LOCC. However, the most challenging aspect of entan- 
glement is that it cannot be 'built in parts', that is, the entanglement of 
N parts is not a sum or a simple function of the entanglement of M(< A^) 
partite subsystems j2]. 

The concept of entanglement has played a crucial role in the devel- 
opment of quantum physics. In the early days entanglement was mainly 
perceived as the qualitative feature of quantum theory that most strikingly 
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distinguishes it from our classical intuition. The subsequent development 
of Bell's inequalities made this distinction quantitative, and therefore ren- 
dered the nonlocal features of quantum theory accessible to experimental 
verification P, IH, [3]. Bell's inequalities may indeed be viewed as an early 
attempt to quantify quantum correlations that are responsible for the coun- 
terintuitive features of quantum mechanically entangled states. At the time 
it was almost unimaginable that such quantum correlations between dis- 
tinct quantum systems could be created in well controlled environments. 
However, the technological progress of the last few decades means that we 
are now able to coherently prepare, manipulate, and measure individual 
quantum systems, as well as create controllable quantum correlations. In 
parallel with these developments, quantum correlations have come to be 
recognized as a novel resource that may be used to perform tasks that are 
either impossible or very inefficient in the classical realm. These devel- 
opments have provided the seed for the development of modern quantum 
information science. 

Given the new found status of entanglement as a resource it is quite 
natural and important to discover the mathematical structures underlying 
its theoretical description. We will see that such a description aims to 
provide answers to three questions about entanglement, namely (1) its 
detection and classification, (2) its manipulation and, (3) its quantification. 

In this thesis, we deal with the first and the third problem. We have 
used two approaches, the combinatorial and the geometric (Bloch repre- 
sentation) approaches for studying the detection problem and we give a 
geometric measure for quantifying the entanglement of multipartite pure 
states, we extend it to mixed states by convex roof construction. Our mea- 
sure satisfies all properties expected of a good measure of entanglement. 

In order to fathom and use entanglement and its role in various quantum 
phenomena, we must be able to say what entanglement is, and how we 
actually use it. In any quantum communication experiment we would 
like to be able to distribute quantum particles across distantly separated 
laboratories. Perfect quantum communication is essentially equivalent to 
perfect entanglement distribution. If we can transport a qubit without 
any decoherence [B], then any entanglement shared by that qubit will also 
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be distributed perfectly. Conversely, if we can distribute entangled states 
perfectly then with a small amount of classical communication, we may use 
teleportation jTj to perfectly transmit quantum states. However, in any 
experiment involving these processes, the effects of noise will inevitably 
impair our ability to send quantum states over long distances. A way of 
trying to overcome this problem is to distribute quantum states by using 
the noisy quantum channels that are available, but then to try and combat 
the effects of this noise using higher quality local quantum processes in the 
distantly separated laboratories. However, there is no reason to make the 
operations of separated laboratories totally independent. It turns out that 
the ability to perform classical communication is vital for many quantum 
information protocols - a prominent example being teleportation. 

We have loosely described entanglement as the quantum correlations 
that can occur in many-party quantum states. This leads to the question 
what differentiates quantum correlations from classical correlations? The 
distinction between 'quantum' effects and classical effects is frequently a 
cause of heated debate. However, in the context of quantum information a 
precise way to define classical correlations is via LOCC operations. Clas- 
sical correlations can be defined as those that can be generated by LOCC 
operations. If we observe a quantum system and find correlations that can- 
not be simulated classically, then we usually attribute them to quantum 
effects, and hence label them quantum correlations. The entanglement is a 
resource because it hfts the so-called LOCC constraint, i.e. entanglement 
and LOCC together can perform tasks that cannot be accomplished by 
LOCC alone. Using LOCC-operations as the only other tool, the inherent 
quantum correlations of entanglement are required to implement general, 
and therefore nonlocal, quantum operations on two or more parts P, . As 
LOCC-operations alone are insufficient to achieve these transformations, 
we conclude that entanglement may be defined as the sort of correlations 
that may not be created by LOCC alone. 

Entanglement has proved to be a vital physical resource for various kinds 
of quantum-information processing, including quantum state teleportation 
fWl [TT] , cryptographic key distribution |[12] , classical communication over 
quantum channels [ISl [H], [15], quantum error correction p!B], quantum 
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computational speedups pJTj, and distributed computation |[TH1 [12]. Fur- 
ther, entanglement is expected to play a crucial role in the many particle 
phenomena such as quantum phase transitions, transfer of information 
across a spin chain [201 IZIj etc. Therefore, quantification of entanglement 
of multipartite quantum states is fundamental to the whole field of quan- 
tum information and in general, to the physics of multicomponent quantum 
systems. 

Whereas the entanglement in pure bipartite state is well understood, 
the understanding of entanglement in mixed bipartite state is far from 
complete. In section 1.1, we review the entanglement of bipartite quantum 
systems. We will state the available measures and criteria for detecting 
entanglement for both bipartite pure and mixed states. In section 1.2, we 
deal with multipartite entangled states. In section 1.3, we briefly summa- 
rize graph theory and density matrix of a graph. In section 1.4, we discuss 
a geometric approach i.e. Bloch representation of quantum states. In sec- 
tion 1.5, we close the chapter by giving some basic multilinear algebra. 
The material in section 1.3, 1.4 and 1.5 forms a background for chapters 
2, 3, 4, and 5. 



1.1 Bipartite Entanglement 

In this section, we define the entanglement in bipartite quantum states. We 
review the work that has been done in the bipartite systems in connection 
with the detection and quantification of entanglement. 

Consider a system consisting of two subsystems. Quantum mechanics 
associates to each subsystem a Hilbert space. Let Ha and Hb denote these 
two Hilbert spaces, let \i)A (where i = 1,2,3, •••) represent a complete 
orthonormal basis for Ha and Ij)^ (where j = 1,2,3, •••) a complete 
orthonormal basis for Hb- Quantum mechanics associates with the system, 
i.e. the two subsystems taken together, the Hilbert space spanned by the 
states \i)A® In the following, we will drop the tensor product symbol 
and write \i)A 'S> as \i)A\j)B, and so on. Any linear combination of 
the basis states \i)A\j)B is a state of the system, and any state \i^)AB of 
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the system can be written [22 



AB = 



where the Cij cire complex coefficients, we take \i/j)ab to be normahzed, 
hence Xl^j IQjI^ = 1- 

If we can write \iP)ab = \i^^^^)A\i^^^'')B, we say the \iP)ab is product 
state (separable state). If \i/j)ab is not a product state, we say that it is 
entangled. 

By using local operators and classical communication (LOCC) any state 
\i^)AB of two subsystems A and B can be transformed to the form [23, E] 

k 

\iP)ab = ^ k < dim{HA ® 'Hb), 

where the positive coefficients di are called Schimdt coefficients. The state 
is entangled if at least two coefficients do not vanish. Pure entangled state 
contains quantum correlation which can not be simulated by any classical 
tools. A fundamental Theorem was proved by Bell |[3|, who showed that 
if the constraint of locality was imposed on the hidden variables, then 
there was an upper bound on the correlations of results of measurements 
that could be performed on the two distant systems. That upper bound, 
mathematically expressed by Bell's inequality [3], is violated by some state 
in quantum mechanics, thus the state contains quantum correlation which 
is Non-local property of quantum state [21] . 

However, in real conditions, owing to interaction with the environment, 
called decoherence, we encounter mixed states rather than pure ones. A 
mixed state is a classical mixture of pure quantum states [22]. These 
mixed states can still possess some residual entanglement. A mixed state 
is considered to be entangled if it is not a mixture of product states [23 • 



In mixed states the quantum correlations are weakened and hence the 
manifestations of mixed state entanglement can be very subtle [26]. It is 
difficult to apply directly the above definition of entanglement of mixed 
states to know whether a quantum state is entangled or not, because the 
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mixed state contains both classical and quantum correlations, and can be 
prepared using infinite possible ensembles. 

For pure states, it is easily shown that the CHSH inequality is violated 
by any nonfactorable state [271, EB] , while on the other hand a factorable 
state trivially admits a (contextual) LHV model jlj . 

For mixed states, Werner [25] constructed a density matrix p^j for a 
pair of spin-s particles. Werner's state can not be written as a sum of 
direct products of density matrices, ^jCjpf pf, where A and B refer 
to the two distant particles and j runs over the states in the ensemble. 
Therefore, genuine quantum correlations are involved in p^. Nevertheless, 
for any pair of ideal local measurements performed on the two particles, 
the correlations derived from p^j^ not only satisfy the CHSH inequality, but, 
as Werner showed [25], it is possible to introduce an explicit LHV model 
that correctly reproduces all the observable correlations for these ideal 
measurements [24]. Thus for mixed states entanglement and nonlocality 
are two different resources. 



1.1.1 Quantification and Detection of Bipartite Entanglement 

Given the wide range of tasks that exploit entanglement, one might try to 
define entanglement as 'that property which is exploited in such protocols'. 
However, there is a whole range of such tasks, with a whole range of possible 
measures of success. This means that situations will almost certainly arise 
where a state pi is better than another state p2 for achieving one task, but 
for achieving a different task p2 is better than pi . Consequently using a 
task-based approach for quantifying entanglement will certainly not lead 
to a single unified perspective. However, despite this problem, it is possible 
to assert some general statements which are valid regardless of what your 
favorite use of entanglement is, as long as the key set of 'allowed' operations 
is the LOCC class. This guides us as to how to approach the quantification 
of entanglement, and so we will state some of this statement [I] : 

i) Separable states contain no entanglement 

ii) All non-separable states allow some tasks to be achieved better than 
by LOCC alone, hence all nonseparable states are entangled. 
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iii) The entanglement of states does not increase under LOCC transfor- 
mations. 

iv) Entanglement does not change under Local Unitary operations. 

v) There are maximally entangled states. 

These properties give us some hints for the quantification of entangle- 
ment from the perspective of LOCC transformations in the asymptotic 
limit. However, one can try to salvage the situation by taking a more ax- 
iomatic approach. One can define real valued functions that satisfy the 
basic properties of entanglement that we outlined above, and use these 
functions to attempt to quantify the amount of entanglement in a given 
quantum state. 

We will now discuss and present a few basic axioms that any measure 
of entanglement should satisfy |T1, [29] . 

1- A bipartite entanglement measure E{p) is a mapping from density 
matrices into positive real numbers, p E{p) G M""", defined for states of 
arbitrary bipartite systems. A normalization factor is also usually included 
such that the maximally entangled states of two qudits has E{\iIj^)) = log d. 

2- E{p) = if the state is separable. 

3- E{p) does not increase on average under LOCC. 

4- For pure state |'?/^)('?/^| the measure reduces to the entropy of entan- 
glement defined bellow. 

Any function E satisfying the first three conditions is called an entan- 
glement monotone. The term entanglement measure will be used for any 
quantity that satisfies axioms 1,2 and 4, and also does not increase under 
deterministic LOCC transformations. Frequently, some authors also im- 
pose additional requirements for entanglement measures: [i29] convexity, 
additivity and continuity. 

The study of the LOCC transformation of pure states has so far en- 
abled us to justify the concept of maximally entangled states and has also 
permitted us, in some cases, to assert that one state is more entangled 
than another. However, we know that exact LOCC transformations can 
only induce a partial order on the set of quantum states. The situation is 
even more complex for mixed states, where even the question of when it is 
possible to LOCC transform one state into another is a difficult problem 
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with no transparent solution at the time of writing. 

All this means that if we want to give a definite answer as to whether 
one state is more entangled than another for any pair of states, it will 
be necessary to consider a more general setting. In this context a very 
natural way to compare and quantify entanglement is to study LOCC 
transformations of states in the so called asymptotic regime. Instead of 
asking whether for a single pair of particles the initial state p may be 
transformed to a final state a by LOCC operations, we may ask whether for 
some large integers m, n we can implement the 'whole-sale' transformation 
p®n — ^ (j*^"^. The largest ratio m/n for which one may achieve this would 
then indicate the relative entanglement content of these two states. In 
this setting we consider imperfect transformations between large blocks of 
states, such that in the limit of large block sizes the imperfections vanish. 

Such an asymptotic approach will alleviate some of the problems that 
we encountered in the case of manipulation of single bi-partite states. It 
turns out that the asymptotic setting yields a unique total order on bi- 
partite pure states, and as a consequence, leads to a very natural measure 
of entanglement that is essentially unique. We will start by defining first 
entanglement measure - the entanglement cost Ec{p) fSO, EH, [321 E31 134] . 
For a given state this measure quantifies the maximal possible rate r 
at which one can convert blocks of two-qubit maximally entangled states 
into output states that approximate many copies of such that the ap- 
proximations become vanishingly small in the limit of large block sizes. 
Ec{p) measures how many maximally entangled states are required to cre- 
ate copies of p by LOCC alone, we can ask about the reverse process: at 
what rate may we obtain maximally entangled states from an input supply 
of states of the form p. This process is known as entanglement distillation 
[ j5Ul ES] ( usually reserved for the pure state). Ed{p) tells us the rate at 
which noisy mixed states may be converted back into the singlet state by 

LOCC. Given these two entanglement measures it is natural to ask whether 
? _ 

Ec = Ejjj i.e. whether entanglement transformations become reversible in 
the asymptotic limit. This is indeed the case for pure state transforma- 
tions where E^ip) and Ec{p) are identical and equal to the entropy of 
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entanglement [3D]. Ejj{p) and Ec{p) are not equal for mixed states and 



also difficult to compute for mixed states, except for some simple but very 
special states jMl E 



Thus, we need a related measure of entanglement, which is the entan- 
glement of formation [22] • For a mixed state p this measure is defined as 



Ef{p) = inf{Y,iPtE{\ilJi){'ilJi\) : p = (V^i |}- Given that this mea- 

sure represents the minimal possible average entanglement over all pure 
state decompositions of p, where E = S{trB{\'ilj){ij\} is entropy 
of entanglement for pure states [32|. The variational problem that de- 



fines Ep is extremely difficult to solve in general and at present one must 
either resort to numerical techniques for general states [38], or restrict 
attention to cases with high symmetry (e.g. [3Sl HDl E]), or consider 
only cases of low dimensionality. Quite remarkably a closed form solu- 
tion is known for bipartite qubit states [331 EH EB]- This exact formula 
is based on the often used two-qubit concurrence. For general bi-partite 

qubit states it has been shown that [M] Ef{p) = S{ ^^'^^^^ ^—^), with 
S{x) = —xlog^x — (1 — x)log2{^- — x), C being the concurrence. For higher 
dimensional systems this connection breaks down, in fact there is not even 
a unique definition of the concurrence [121 33] • Another important class of 
measures is entanglement measures from convex roof constructions. The 
entanglement of formation Ep is an important example of the general con- 
cept of a convex roof construction. The convex roof / of a function / 
is defined as the largest convex function that is bounded from above by 
the function / for all arguments [42|, HI] . The importance of the convex 
roof method is based on the fact that it can be used to construct entan- 
glement monotones from any unitarily invariant and concave function of 
density matrices [45]. Also various such quantities have been proposed 



over the years, such as the relative entropy of entanglement [461 ESI 147] , 
the squashed entanglement [i48| and Logarithmic Negativity [12| . Mintert 
et al. [5D| found a lower bound on I-concurrence [13| which is simpler to 
estimate than the I-concurrence itself. Another attempt of generalizing the 
concurrence for mixed states in higher dimensions was made by Badziag 
et al. in [51]. Yet another proposal to deal with mixed states in higher 
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dimensions is presented in [52] • For more details we refer the reader to |[I] . 
Unfortunately, all these measures are difficult to implement experimentally 
and they require substantial efforts to estimate. 

On the other side, there are attempts to understand the separability 
problem, which is to identify the states that contain classical correlations 
only (or no correlations at all). These states are termed separable states, 
and their mathematical characterization has been formulated by Werner 
[25]. We call this the problem of entanglement detection. 

A major step in the characterization of the separable states was done 
by Peres [53] and the Horodecki family [M] • Peres provided very powerful 
necessary condition for separability. Later on, Horodeckis demonstrated 
that this condition is also sufficient for composite Hilbert spaces of dimen- 
sion 2x2 and 2 x 3. A density matrix that verifies Peres criterion is termed 
"PPT" for positive partial transpose. In general, there exist PPT states p 
which are not separable in TiA^HB spaces {dim{7iA) > 2, dim{HB) > 4 or 
dim{HA) > 3) [55j. The PPT entangled states have been termed "bound 
entangled states" to distinguish them from the "free entangled states". 
"Bound entangled states" are entangled, however, no matter how many 
copies of them we have, these states cannot be "distilled" via local opera- 
tions and classical communication to the form of pure maximally entangled 
states [5^. We encounter thus new problems such as: How can one distin- 
guish a separable state p from a PPT entangled state pi Are all non-PPT 
states (NPPT states) "free entangled" i.e. distillable? 

From the end of last century, there has been a growing effort in searching 
for necessary and sufficient separability criteria. Several necessary condi- 
tions for separability are known: Werner has derived a condition based on 
the analysis of local hidden variables (LHV) models and the mean value 
of the, so-called, flipping operator [25], the Horodeckis have proposed a 
necessary criterion based on the so-called a-entropy inequalities [57j, etc. 
A general necessary and sufficient condition for separability was discov- 
ered by the Horodecki family in terms of positive maps [M]. A map is 
defined positive if it maps positive operators into positive operators. Later 
on the reduction criterion of separability was introduced [EE, [59] . Viola- 
tion of this criterion is sufficient for entanglement to be free. Sufficient 
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conditions for separability are also known. In [60] it was proved that any 
state close enough to the completely random state ctq = I/NM is separa- 
ble. In ||55], in which the first explicit examples of entangled states with 
PPT property were provided. Another necessary criterion of separability 
was formulated which is the so-called range criterion. The analysis of the 
range of the density matrices, initiated by P. Horodecki, turned out to be 
very fruitful, leading, in particular, to the algorithm of optimal decompo- 
sition of mixed states into the separable and inseparable part jMl E2] , and 
to systematic methods of constructing examples of PPT entangled states 
with no product vectors in their range, using either so-called unextendible 
product bases (UPB's) p3l Ell, or the method described in jHS]. Also, 
considerable progress in the study of PPT entangled states has been made 
[ 1661 [67j . Lewenstein et al. employ the idea of "subtracting projectors on 
product vectors" [SH E2]- They introduced the 'edge' state, which has a 
property that no projection onto the product state can be subtracted from 
it, keeping the rest positive definite and PPT. They mentioned a different 
approach to the entanglement problem, based on the so-called entangle- 
ment witnesses. An entanglement witness is an observable W that reveals 
the entanglement of an entangled density matrix p. Rudolph [SB]; Chen 
and Wu [69j discovered new criterion called computable cross norm (CCN) 
criterion or matrix realignment criterion. Quite remarkably the realign- 
ment criterion has been found to detect some of PPT entanglement. It 
also provides nice lower bound on concurrence function [[70j. General sep- 
arability criteria based on local uncertainty relation valid both for discrete 
and continuous variables have been introduced in [71], [72] • Further it has 
been shown [i73j that PPT entanglement can be detected by means of local 
uncertainty relations introduced in [72]. This approach has been further 
developed and simplified by Giihne [7^ and developed also in entropic 
terms [75] . Recently, Braunstein et al. [73, [77] have initiated a new ap- 
proach towards the mixed state entanglement by associating graphs with 
density matrices and understanding their classification using these graphs. 
Hildebrand et al. [TS] testified that the degree condition [7B] is equiva- 
lent to the PPT criterion. Sixia Yu et al. [29} have given a new family of 
entanglement witnesses and corresponding positive maps that are not com- 
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pletely positive based on local orthogonal observables. de Vicente [i80j has 
introduced a new approach to study the separability of bipartite quantum 
systems in arbitrary dimensions using the Bloch representation of their 
density matrix. He has obtained analytical lower bounds on the concur- 
rence of bipartite quantum systems in arbitrary dimensions related to the 
violation of separability conditions based on local uncertainty relations and 
on the Bloch representation of density matrices [BI] • Very recently Githne 
et al. |[H2] have proposed a unifying approach to the separability problem 



which uses a representation of a quantum state by a covariance matrix of 
locally measurable observables and they have proposed nonlinear witness 
[ IH5] . Despite many efforts and seminal results obtained in the recent years, 
the problem of separability of bipartite mixed states remains essentially 
open. 



1.2 Multipartite Entanglement. 

In this section, we deal with entanglement in multipartite quantum states. 
We state the main difference between multipartite and bipartite entangled 
states. We review the work that has been done in connection with the 
detection and quantification of multipartite entanglement. 

Multiparticle entanglement is genuinely different from entanglement in 
quantum systems consisting of two parts. To understand what is so dif- 
ferent consider, say, a quantum system that is composed of three qubits. 
Each of the qubits is to be held by one of three laboratories distantly 
separated. It may come as quite a surprise that states of such composite 
quantum systems may contain tripartite entanglement, while at the same 
time showing no bi-partite entanglement at all. In contrast to the bipartite 
setting, there is no longer a natural "unit" of entanglement, the role that 
was taken by the maximally entangled state of a system of two qubits. 
Quite strikingly, the very concept of being maximally entangled becomes 
void. Instead, we will see that in two ways there are "inequivalent kinds 
of entanglement". Consider multipartite entanglement of pure quantum 
states. A theory of entanglement should not discriminate states that dif- 
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fer only by a local operation. Here, "local operation" can mean merely 
a change of local bases (LU operations) or, else, general local quantum 
operations assisted by classical communication, that are either required to 
be successful at each instance (LOCC) or just stochastically (SLOCC). For 
each notion of locality, local unitary operation (LU) or Local operations 
and classical communication (LOCC) or just stochastic-LOCC (SLOCC), 
the questions that have to be addressed are how many equivalence classes 
exist, how are they parameterized and how can one decide whether two 
given states belong to the same class? 

For the case of bi-partite qubit states, two quantum states are LU- 
equivalent if and only if their respective Schmidt normal forms coincide. 
All classes are parameterized by only one real parameter. Some simple 
parameter counting arguments show that in the case of N-qubit systems 
the situation must be vastly more complex. Indeed, disregarding a global 
phase, it takes 2^^^ — 2 real parameters to fix a normalized quantum state 
in 7"^ = (C^)*^^. The group of local unitary transformations SU{2) x • • • x 
SU{2) on the other hand has 3N real parameters |j2]. Therefore, one needs 
at least 2^+^ — 3N — 2 real numbers to parameterize the sets of inequivalent 
pure quantum states [84j. This lower bound turns out to be tight [85]. It 
is a striking result that the ratio of non-local to local parameters grows 
exponentially in the number of systems. In particular, the finding rules out 
all hopes of a generalization of the Schmidt normal form. A general pure 
tripartite qubit state, say, cannot be cast into the form sin6\000)+cos0\lll) 
by the action of local unitaries jS]- Considerable effort has been undertaken 
to describe the structure of LU-equivalence classes by the use of invariants 
or normal forms [841, [851, [861, [87|, [88] . Acm et al. [89] have proved for any 
pure three-qubit state the existence of local bases which allow one to build 
a set of five orthogonal product states in terms of which the state can be 
written in a unique form. This leads to a canonical form which generalizes 
the two-qubit Schmidt decomposition. It is uniquely characterized by the 
five entanglement parameters. When one deals with SLOCC operations, 
the group of SLOCC, SL{C^) x • • • x SLiC^) has 2^+^ - 6iV - 2 parameters 
that are necessary to label SLOCC equivalence classes of qubit systems. It 
turns out that the three-qubit pure states are partitioned into a total of Six 



14 



1. Introduction and Overview 



SLOCC-equivalent classes [HD]- The picture is complete for three-qubit : 
any fully entangled state is SLOCC-equivalent to either \GHZ) or \W) [90]. 
Three-qubit W-states and GHZ-states have already been experimentally 
realized, both purely optically using postselection [SH [22] and in ion traps 
[93]. The two states behave differently, however, if a system is traced 
out. Specifically, tracing out the first qubit of the GHZ state will leave the 
remaining systems in a complete mixture. For \W) will leave the remaining 
systems in a mixed entangled bipartite state. Thus, the entanglement 
of \W) is more robust under particle loss than the one of \GHZ) [90| . 
From point of view of asymptotic manipulation of multipartite quantum 
states, there is no longer a single essential ingredient as in bipartite the 
maximally entangled state or EPR-state, but many different ones. In the 
multi-particle case, however, it is meaningful to introduce the concept of a 
minimal reversible entanglement generating set {MREGS). An MREGS 
5 is a set of pure states such that any other state can be generated from 
S by means of reversible asymptotic LOCC. It must be minimal in the 
sense that no set of smaller cardinahty possesses the same property j94l 
^51 ES]- Yet, it can be shown that merely to consider maximally entangled 
qubit pairs is not sufficient to construct an MREGS P3] . To find general 
means for constructing MREGS constitutes one of the challenging open 
problems of the field: as long as this question is generally unresolved, the 
development of a "theory of multi-particle entanglement" in the same way 
as in the bi-partite setting seems infeasible. 

Regardless of whether there is a unit of multipartite entanglement, re- 
searchers have tried to find some measure of multipartite entanglement of 
pure states. Recently , Meyer and Wallach [i97j have defined a polynomial 
measure which is scalable, i.e. which applies to any number of spin-1/2 
particles. Wong and Christensen [i98j have proposed a potential measure 
of a type of entanglement of pure states of N-qubits, the N-tangle. For a 
system of two qubits the N-tangle is equal to the square of the concurrence, 
and for systems of three qubits it is equal to the "residual entanglement" . 
The geometric measure of entanglement [99] makes use of a geometric dis- 
tance to the set of product state: Eceo = 7Tim|||'?/^)(?/^| — cr||2, where ||.||2 
is the Hilbert-Schmidt norm, and the minimum is taken over all product 
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states. The relative entropy of entanglement in the multipartite setting 
is defined as the minimal distance of a given state to the set of fully sep- 
arable states, quantified in terms of the quantum relative entropy jlOO) . 
There are also many measures are defined for multipartite entanglement 
of pure and mixed states, as in IM], dOl [ffl MM, ME, W, [Ml W ■ 



Recently Lamata et al. fl09j have proposed an inductive procedure to 
classify A^-partite spin-1/2 entanglement under stochastic local operations 
and classical communication provided such a classification is known for 
N — 1 qubits. The method is based upon the analysis of the coefficient 
matrix of the state in an arbitrary product basis. For mixed state the 
classification scheme is based on separability properties |110f . At the low- 



est level there is the class of states that can be prepared using LOCC 
alone. Its members are called fully separable and can be written in the 
form p = YliiPi{Pi ® P2® ' " ® Pn)- Evidently, states of this kind do not 
contain entanglement. A state is referred to as k-separable, if it is fully 
separable considered as a state on some k-partite split. By the use of this 
terminology, the separability classes can be brought into a hierarchy, where 
/c-separable classes are considered to be more entangled then /-separable 
ones for /c < /. States that are not separable with respect to any non-trivial 
split are fully inseparable. The number of all splits of a composite system 
grows exorbitantly fast with the number of its constituents. One is natu- 
rally tempted to reduce the complexity by identifying redundancies in this 
classification. After all, once it is established that a state is fully separa- 
ble, there is no need to consider any further splits. For three systems. The 
five possible splits (1-2-3, 12-3, 1-23, 13-2, 123) have already been known 
for pure state in the above discussion. It is a counter-intuitive fact that 
there are mixed states that are separable with respect to any bi-partite 
split but are not fully separable j63]. An analogous phenomenon does not 
exist for pure states. The following sub-classes of the set of bi-separable, 
i. e. 2-separable, states are all non-empty fllOj . 



• 1-qubit bi-separable states with respect to the first system are sepa- 
rable for the split 1-23 but not for 12-3 or 13-2. 



• 2-qubit bi-separable states with respect to the first and second system 
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are separable for the split 1-23 and 2-13, but not for 12-3. 

• 3-qubit bi-separable states are separable with respect to any bi-partite 
split but are not fully separable. 

Together with the inseparable states and the fully separable ones, the 
above sets constitute a complete classification of mixed three qubit states 
modulo system permutations jlllj . The quantification of entanglement of 



multipartite mixed states is void. The most measures are taken as convex 
roof constructions PHI 11071 1105] . There is some progress in the detection of 
entanglement of multipartite mixed state. Kai Chen and Ling- An Wu jll2l 



have generalized partial transposition separability criterion for the density 
matrix of a multipartite quantum system. This criterion contains as special 
cases the famous Peres-Horodecki criterion [53] and the realignment crite- 



rion p2l EBI- Xiu-Hong Gao et al. fll3j have derived an analytical lower 



bound for the concurrence of tripartite quantum mixed states [107]. A 
functional relation is established relating concurrence and the generalized 
partial transposition [112j . Chang-Shui Yu and He-shan Song [114] have 



presented a method to construct full separability criterion for tripartite sys- 
tem of qubits. Later on, they have generalized it to the higher-dimensional 
systems. The above criteria need a complete quantum state tomography, 
this can be a costly procedure. It may be desirable to detect entangle- 
ment without the need of acquiring full knowledge of the quantum state. 
This is where entanglement witnesses come into play. Recently, Korbicz 
et al. [115j have derived spin squeezing inequalities that generalize the 



concept of the spin squeezing parameter and provide necessary and suf- 
ficient conditions for genuine 2-, or 3-qubit entanglement for symmetric 
states, and sufHcient condition for A^-qubit states. The inequalities have 
a clear physical interpretation as entanglement witnesses. Also Usha Devi 
et al. [116] have shown that higher order inter-group correlations involv- 



ing even number of qubits are necessarily positive semidefinite for separable 
symmetric A-qubit states. Toth and Giihne [117] have presented entangle- 



ment witnesses for detecting genuine multiqubit entanglement. The gener- 
alized Bell-type inequality is used to characterize and detect multipartite 
entanglement. Dili |118j has studied the relation between distillability 
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of multipartite states and violation of Bell's inequality. He proved that 
there exist multipartite bound entangled states that violate a multipartite 
Bell inequality. This implies that (i) violation of Bell's inequality is not a 
sufficient condition for distillability and (ii) some bound entangled states 
cannot be described by a local hidden variable model. Later on, Acm |119l 



has proved that for all the states violating this inequahty there exists at 
least one splitting of the parties into two groups such that some pure-state 
entanglement can be distilled. We saw that, for bipartite systems, bound 
entanglement is clearly defined as it involves only two spatially separated 
parties and a necessary and sufficient condition for distillability of bipartite 
quantum states is known jSS]- In a multiparty setting , however, due to 
several distinct spatially separated configurations, the definition of bound 
entanglement is not unique. A multipartite quantum state is said to be 
bound entangled if there is no distillable entanglement between any subset 
as long as all the parties remain spatially separated from each other. When, 
however, one also allows some of the parties to group together and perform 
local operations collectively, two qualitatively different classes of bound en- 
tanglement arise: (a) activable bound entangled (ABE) states. The states 
that are not distillable when every party is separated from every other but 
become distillable, if certain parties decide to group together [ |120l 11211 . 



This implies that there is at least one bipartite partition/cut where the 
state is negative under partial transposition (NPT) [53] • Such states have 
been also referred to as unlockable bound entangled (UBE) states in the 
literature, (b) Nonactivable Bound Entangled states- states that are not 
distillable under any modified configuration as long as there are at least 
two spatially separated groups. In other words, such states are always 
positive under partial transposition across any bipartite partition [6^ . For 
unlockable bound entangled states, we refer the reader to fl20l 11211 11221 . 



Zukowski et al. |123j have derived a single general Bell inequality which 



is a sufficient and necessary condition for the correlation function for 
particles to be describable in a local and realistic picture, they also de- 
rived a necessary and sufficient condition for an arbitrary A-qubit mixed 
state to violate this inequality. Later on, in |124j it was shown that there 
exist pure entangled A > 2 qubit states that do not violate any Bell in- 
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equality for N particle correlation functions for experiments involving two 
dichotomic observables per local measuring station. Laskowski et al. fl25l 



have shown that the generalized Bell-type inequality, explicitly involving 
rotational symmetry of physical laws, is very efficient in distinguishing be- 
tween true A^-particle quantum correlations and correlations involving less 
particles. This applies to various types of generalized partial separabilities. 
Very recently, Badziag et al. jl26j have presented an intuitive geometrical 
entanglement criterion. It allows formulation of simple and experimentally 
friendly sufficient conditions for entanglement of A^-qubits. Li et al. [ jl27l 



have investigated the separability of arbitrary dimensional tripartite sys- 
tems. By introducing a new operator related to transformations on the 
subsystems a necessary condition for the separability of tripartite systems 
is presented. 



1.3 Graphs and density matrix of a graph 

In this section, we give a brief summary for graphs, which is necessary for 
chapter 2 and 3. We will also introduce the definition of density matrix of 
a graph. 

A graph G consists of a vertex set V{G) and an edge set E{G), where 
an edge is an unordered pair of distinct vertices of G. If vu is an edge, 
then we say that v and u are adjacent or that w is a neighbour of -y, and 
denote this by writing -u ~ w. A loop is an edge whose endpoints are equal. 
Two graph G and H are equal if and only if they have the same vertex 
set and the same edge set. Two graphs G and H are isomorphic if there 
is a bijection, (p say, from V{G) to V{H) such that -u ~ w in G if and 
only if (p{v) ~ ^{u) in H. We say that (p is an isomorphism from G to 
H. Since (f is bijection it has an inverse, which is an isomorphism from 
H to G. If G and H are isomorphic, then we write G = H. A graph 
is called complete if every pair of vertices are adjacent, and the complete 
graph on n vertices is denoted by Kn- A graph with no edges (but at least 
one vertex) is called empty. The graph with no vertices and no edges is the 
null graph. A subgraph of a graph G is a graph H such that V{H) C V{G), 
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E{H) C E{G). If V{H) = V{G), we call H a spanning subgraph of G. 
Any spanning subgraph of G can be obtained by deleting some of the edges 
from G. The number of spanning subgraphs of G is equal to the number 
of subsets of EiG). 

A subgraph H oi G is an induced subgraph if two vertices of V{H) are 
adjacent in H if and only if they are adjacent in G. Any induced subgraph 
of G can be obtained by deleting some of the vertices from G, along with 
any edges that contain a deleted vertex. The number of induced subgraphs 
of G is equal to the number of subsets of V{G). 

A clique is a subgraph that is complete. It is necessarily an induced 
subgraph. A set of vertices that induces an empty subgraph is called an 
independent set. A path of length r from -u to w in a graph is a sequence of 
r + 1 distinct vertices starting with v and ending with u such that consec- 
utive vertices are adjacent. If there is a path between any two vertices of 
a graph G, then G is connected^ otherwise disconnected. Alternatively, G 
is disconnected if we can partition its vertices into two nonempty sets, R 
and 5, say, such that no vertex in R is adjacent to a vertex in S. in this 
case we say that G is the disjoint union of two subgraphs. 

A cycle is a connected graph where every vertex has exactly two neigh- 
bours. The smallest cycle is the complete graph K^,. An acyclic graph is 
a graph with no cycles. A connected acyclic graph is called a tree, and an 
acyclic graph is called a forest., since each component is a tree, a spanning 
subgraph with no cycles is called a spanning tree. A graph has a spanning 
tree if and only if it is connected. A star Ki^ni which consists of a single 
vertex with n neighbours. 

An isomorphism from a graph G to itself is called an automorphism of G. 
An automorphism is therefore a permutation of the vertices of G that maps 
edges to edges and nonedges to nonedges. The set of all automorphisms 
of G forms a group, which is called the automorphism group of G and 
denoted by Aut{G). The symmetric group Sym{V) is the group of all 
permutations of a set F , and so the automorphism group of G is a subgroup 
of Sym{V{G)). 

The valency of a vertex v is the number of neighbours of and the 
maximum and minimum valency of a graph G are the maximum and min- 
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imum values of the valencies of any vertex of G. A graph in which every 
vertex has equal valency k is called regular of valency k or /c-regular. The 
distance dciv^u) between two vertices v and w in a graph G is the length 
of the shortest path from v to u. A leaf is a vertex of degree (valency) 1. 
An isolated vertex has degree 0. 

The complement G of a graph G has the same vertex set as G, where 
vertices v and u are adjacent in G if and only if they are not adjacent in 
G. A mapping / from V{G) to V{H) is a homomorphism if f{v) and f{u) 
are adjacent in H whenever v and u are adjacent in G. A graph G is called 
bipartite if its vertex set can be partitioned into two parts Vi and V2 such 
that every edge has one end in Vi and one inV2. A proper colouring of a 
graph G is a map from V{G) into some finite set of colours such that no 
two adjacent vertices are assigned the same colour. If G can be properly 
coloured with a set of k colours, then we say that G can be properly k- 
coloured. The least value of k for which G can be properly /c-coloured is 
the chromatic number of G, and is denoted by x{G)- The set of vertices 
with a particular colour is called a colour class of the colouring, and is an 
independent set. 

A homomorphism from a graph G to itself is called an endomorphism^ 
and the set of all endomorphisms of G is the endomorphism monoid of G. 
A monoid is a set that has an associative binary multiplication defined on 
it and an identity element. 

A line graph of a graph G is graph L{G) with the edges of G as its 
vertices, and where two edges of G are adjacent in L{G) if and only if 
they are incident in G. The star Ki^^i: consists of a single vertex with n 
neighbours. For more details see |128l I129| . 

The Cartesian product GHH of two graphs G and H is defined on the 
Cartesian product V{G) x V{H) of the vertex sets of the factors. The 
edge set E{G\3H) is the set of all pairs {{u, v), {x, y)} of vertices for which 
either u = x and {-?;, y} G E{H) or {w, x} G E{G) and v = y. 

The direct product (tensor product) G0H is defined also on the Carte- 
sian product V{G) X V{H) of the vertex sets of the factors. Two vertices 
(^1,^2), {vi,V2) are adjacent when uiVi G E{G) and U2V2 G E{H). 

The adjacency matrix of a graph on n vertices G is an n x n matrix. 
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denoted by M{G), having rows and columns labeled by the vertices of G, 
and ij — th entry defined as follows: 

L^^A^jJM-| if{^,,^j^E(G). (1.1) 

Note that M{G H) = M{G) (g) M{H). 

Two distinct vertices Vi and vj are said to be adjacent if {-u^, -Uj} G E{G). 
The degree of a vertex Vi G ^(C), denoted by 'Ocivi), is the number of edges 
adjacent to Vi. Two adjacent vertices are also said to be neighbours. The 
degree-sum of G is defined and denoted by do = Yl^i=i^G{vi)- Note that 
Dg = The degree matrix of G is an n x n matrix, denoted by 

A(G), having ij — th entry defined as follows: 



[A(G)] 



^G{v^) if i = j; 

ifi^j. [1.2] 



The combinatorial Laplacian matrix of a graph G (for short, Laplacian) 
is the matrix L{G) = A(G) — M{G). Notice that L{G) does not change 
if we add or delete loops from G. According to our definition of graph, 
L(G) ^ 0. Note that L{G ® H) ^ L{G) ® L{H). 

In Standard Quantum Mechanics (that is the Hilbert space formulation 
of Quantum Mechanics), the state of a quantum mechanical system associ- 
ated to the n-dimensional Hilbert space = C"^ is identified with an n x n 
positive semidefinite, trace-one, Hermitian matrix, called a density matrix. 
It is easy to observe that the Laplacian of a graph is symmetric and pos- 
itive semidefinite. The Laplacian of a graph G, scaled by the degree-sum 
of G, has trace one and it is then a density matrix. This observation leads 
to the following definition [77| . 



The density matrix of a graph G is the matrix 



cj{G) = ^L{G). (1.3) 
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1.4 Bloch Representation 

In this section, we discuss the geometric approach to a density matrix via 
its Bloch representation. The determination of a state on the basis of the 
actual measurement (experimental data) is important both for experimen- 
talists and theoreticians. In classical physics, it is trivial because there is 
a one-to-one correspondence between the state and the actual measure- 
ment. On the other hand, in quantum mechanics, where a density matrix 
is used to describe the state, it is generally nontrivial to connect them 
11301 11311 11321 11331 [5] . the Bloch representation of the density matrix can 
be constructed experimently giving the required connection between the 
density matrix and experiments. 

A^-level quantum states are described by density operators, i.e. unit 
trace Hermitian positive semidefinite linear operators, which act on the 
Hilbert space H, ~ C^. The Hermitian operators acting on H, constitute a 
Hilbert space themselves, the so-called Hilbert-Schmidt space denoted by 
HS{H), with inner product {p,(j)hs = Tr{p^a). Accordingly, the density 
operators can be expanded by any basis of this space. In particular, we can 
choose to expand p in terms of the identity operator 1^ and the traceless 
Hermitian generators of SU (N) Aj = 1, 2, • • • , A^^ — 1), 

The generators of SU{N) satisfy the orthogonality relation 

{X^, Xj)hs = Tr{\\j) = 26^j, (1.5) 

and they are characterized by the structure constants of the corresponding 
Lie algebra, fijk and gijk, which are, respectively, completely antisymmetric 
and completely symmetric, 

2 

\Xj = j;^^ij^N + ifijk^k + gijk^k- (1-6) 
The generators can be easily constructed from any orthonormal basis 



(1.4) 
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in Ti fl34 ] . The (orthogonal) generators are given by 

= {ujk, Vjk, wi}, (1.7) 



when i = 1, • • • , N — 1 



\ = wi = sjj^^YS\3){3\ - l\l + !)(/ + 1 < / < iV - 1, 

while for z = iV, • • • , (iV + 2)(iV - 1) /2 

= Ujk = \3){k\ + 
and for i = N{N + l)/2, ■ ■ ■ , - 1 

K = v,k = -i{\j){k\-\k){j\), 

I < j < k < N. 

The orthogonality relation (1.5) implies that the coefficients in (1.4) are 
given by 

n = —Tr{p\i). 

Notice that the coefficient of In is fixed due to the unit trace condition. 
The vector r = {rir2 ■ ■ ■ Tiv2_2^)* G which completely characterizes 

the density operator, is called Bloch vector or coherence vector. The rep- 
resentation (1.4) was introduced by Bloch fl35j in the N = 2 case and 



generalized to arbitrary dimensions in |134j . Any density matrix in two- 



level systems turns out to be characterized uniquely by a three-dimensional 
real vector where the length satisfies 

|A| = V^<1- (1-8) 

Therefore, if we define the Bloch- vector space 5(M^) as a ball with radius 
1: 

5(M3) = {A=(Ai,A2,A3)gM':|A|<1}, 

its element gives an equivalent description of the density matrix with the 
following bijection (one-to-one and onto) map from 5(M'^) to the set of 
density matrices. 

1 1 

A — > p = -I2 + -\(Ji 
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B(M^) is called the Bloch ball, its surface the Bloch sphere and its element 
the Bloch vector. The equality in Eq.(1.8) (i.e., |A| = 1), the surface of 
the ball (the Bloch sphere) which constitutes the set of extreme points of 
Bloch ball, corresponds to the set of pure states, the points interior to the 
Bloch ball correspond to mixed states. It has an interesting appeal from 
the experimentalist point of view, since in this way it becomes clear how 
the density operator can be constructed from the expectation values of the 
operators Aj, 

{\)=Tt{p\) = ^t,. (1.9) 

As we have seen, every density operator admits a representation as in 
Eq.(1.4); however, the converse is not true. A matrix of the form (1.4) is of 
unit trace and Hermitian, but it might not be positive semidefinite, so to 
guarantee this property further restrictions must be added to the coherence 
vector. The set of all the Bloch vectors that constitute a density operator 
is known as the Bloch- vector space 5(M^ ~^). from above discussion it is 
known that in the case N = 2 this space equals the unit ball in and 
pure states are represented by vectors on the unit sphere. The problem of 
determining 5(M^ ~^) when > 3 is still open and a subject of current 
research |136j . However, many of its properties are known. For instance, 



for pure states (p = p) it must hold 



||r||2 = Y ^\ r,rjg,jk = {N -2)rk, (1.10) 

where ||.||2 is the Euclidean norm on In the case of mixed states, 

the conditions that the coherence vector must satisfy in order to repre- 
sent a density operator have been recently provided in [ jl37l 1138] . Re- 



gretfully, their mathematical expression is rather cumbersome. It is also 
known [BEIIID] that B{R^'-^) is a subset of the ball Dr{R^'~^) of radius 



R = \l ^^^2 which is the minimum ball containing it, and that the ball 
D^(IR^'~i) of radius r = ^^J^ is included in 5(M^'-i). that is, 

DM^"-^) C BiM.^"-^) c DRiM.^"-^). (l.li; 
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In the case of bipartite quantum systems of dimensions M x N {H '^^ 
<S)C^) composed of subsystems A and B, we can analogously represent 
the density operators as 

P = -J^i^M ®In + Y^ n\ ®In + ^ SjIm ® Xj + ^ A, (g) Xj), (1.12) 

i j -ij 

where (A-,) are the generators of SU{M) {SU{N)). Notice that r G 
and s G are the coherence vectors of the subsystems, so that 

they can be determined locally, 



Pa = TtbP = J^i^M + X! ^^^^)^ PB = T^Ap = J^^^N + ^ S^\). (1.13) 

i i 

The coefficients tij^ responsible for the possible correlations, form the real 
matrix T G M(M^ — l)x [N'^ — 1), and, as before, they can be easily ob- 
tained by Uj = ^Tr{p\ ® Xj) = ^(A, Xj). 
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In this section, we give some basic of multilinear algebra and Higher-order 
tensors. Higher-order tensor decompositions are in frequent use today in 
a variety of fields including psychometrics [ jl41l 11431 1144] . chemometrics 



1145] , image analysis fl46l 11471 1148] , graph analysis fl49l 1150] , signal pro 



cessing jl51l 1152] and we will use it in separability problem, the two most 
commonly used decompositions of tensor are Tucker [141] and Kruskal 
11421 11431 1144] , which can be thought of as higher-order generalizations of 



the matrix singular value decomposition. 

We start by defining a product of matrices which are useful to us. The 
Khatri-Rao product fl53l 11541 1155] is the columnwise Kronecker product 
(tensor product). The Khatri-Rao product of matrices A G M^^^ and 
B G R"^^^ is denoted by AQ B and its {IJ) x K result is defined by 



AQ B = [a-A 6:1 0:2 <^ 6;2 " " " CL.K ® ^.k], 



(1.14) 
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where a-j, b-j are the jth columns vectors of the matricies A and B. Ob- 
serve that the matrices in a Khatri-Rao product all have the same number 
of columns. 

As an example, let A e R^^'^ and B e R^^^ be as follows 



A = 



1 4 7 10 

2 5 8 11 

3 6 9 12 



B = 



12 3 8 
4 5 6 10 



The Khatri-Rao product of A and B is AQ B ^ 



x4 



AqB = 



AqB = 




10 24 
25 48 
12 27 



3 106:4 
3 116:4 
3 126:4 

80 \ 
100 
88 
110 
96 



V 12 30 54 120 / 



1.5.1 Tensors 

Let A* be an /i X /2 X • • • X tensor over R. The order of A' is N. 
The nth dimension of X is /„. An element of X is specified as 
where ij G {1,2,- •• for j = 1,2, • • • , N . The set of all tensors 
of size Ii X I2 X • • • X In is denoted by §(/i, • • • , In)- Let A", 3^ G 
S(/i,-- - ,In)- The inner product of A* and y is defined as A* • 3^ = 

Correspondingly, the norm of A* , ||A||, is defined as ||A||^ = A*. A* = 

ii=l Z-^Z2=1 ' ' ' 2^iN=l iii2---iN' 

We say A is a unit tensor if 1 1 A 1 1 = 1 . 
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A decomposed tensor is a tensor U G • • • , In) that can be written 

U = u^'Ku^'K...ou^^\ (1.15) 

where o denotes the outer product and each u^^^ G W\ for j = 1, • • • , A^. 
The vectors w^-^^ are caUed the components of U. Ui^i^...i^ = u^^^u^^^ ■ ■ ■ ul^\ 
A decomposed tensor is a tensor of rank one. The set of aU decomposed 
tensors of size Ii x I2 x ■ ■ ■ x 1^ is denoted by T>{Ii, • • • , 1^). 
Lemma 1.5.1: jl56j Let U,V E V where U is defined as in Eq.(1.15) and 



V is defined by 

y = ^(l)oi;(2)o...oi;(^). (1.16) 

Then (a) U.V = Uf^^u^^\v^^\ (b) \\U\\ = n^^i||w(^)||2, and, {c)U + V eV 
if and only if all but at most one of the components of U and V are equal 
(within a scalar multiple). 

Let U, V G V he defined as Eq.(1.15) and Eq.(1.16) respectively with 
11^711 = ||-\/| I = 1. We say that U and V are orthogonal {U ± V) if U.V = 
IljLiU^^\v^^^ = 0. We say that U and V are completely orthogonal {U _Lc V) 
if for every j = 1, 2, • • • , A; u{j) _L v{j). we say that U and V are strongly 
orthogonal {U ±5 F) if [/ ± V and for every j = 1, • • • , A w(j) = ±t'(j) 
or u{j) ± v{j). 

Let A' G S be a tensor 

r 

:\: = Y,^^U^: (1-17) 

where cr > for i = 1, •• • , r and each Ui eV and ||?7|| = 1 for ^ = 1, • • • , r. 



• The rank of X ^ denoted rank(A'), is defined to be the minimal r such 
that X can be expressed as in Eq.(1.17), The decomposition is called 
the rank decomposition. 

• The orthogonal rank of X, denoted rank±{X), is defined to be the 
minimal r such that X can be expressed as in Eq.(1.17) and Ui _L 
Uj for all i j- The decomposition is called the orthogonal rank 
decomposition. 

• The complete orthogonal rank of X, denoted rankj_^{X), is defined 
to be the minimal r such that X can be expressed as in Eq. (1.17) 
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and Ui _Lc Uj for all i ^ j. The decomposition is called the complete 
orthogonal decomposition. 

• The strong orthogonal rank of A:', denoted rankj__,{^), is defined to 
be the minimal r such that (A:') can be expressed as in Eq.(1.17) 
and Ui _Ls Uj for all i j. The decomposition is called the strong 



orthogonal rank decomposition [ jl56] . 

Assume an Nth-order tensor A G (^(7)^x^2 x - x/w^ rpj^^ matrix unfold- 
ing A^n) G (C')^nx(^«+i/„+2 -/w/i/2 -/„-i) contains the element a^^i^..^^ at the 
position with row number in and column number equal to 

{in+l — l)^n+2^n+3 • • • ^nhh ■ ■ ■ ^n-l + " l)^n+3^n+4 • • • ^Nhh ■ ■ ■ ^n-l 
H H(iAr-l)/i/2 . . . 4_l + (^l-l)/2^3 • • • 4-l + (^2-l)^3^4 • • • Hn-1- 

For n = 1, we take the last term in-i = io = iN- 

As an example jl57] . define a tensor T^^) G M^^^^^, by tm = tiu = 





hn — —hu — 1, ^213 — ^311 — ^313 — ^121 — ^122 — ^221 — "^222 — 2, 

^223 = ^321 = t323 = 4, ^113 = ^312 = ^123 = ^22 = 0. The matrix unfolding 
Tj-^^ is given by 



r(3) 



If we refer in general to the vectors of an Nth-order tensor G (C^i^^^x - x/at 
as its "n-mode vectors", defined as the /^-dimensional vectors obtained 
from A' by varying the index and keeping the other indices fixed, then we 
have the following definition. The n-rank of A:', denoted by Rn = rankn{X), 
is the dimension of the vector space spanned by the n-mode vectors |157| . 
The n-mode vectors of A' are the column vectors of the matrix unfolding 

A(^n) and rankn{^) = rank{A(^n))- 

The n-mode product of a tensor 3^ G M'^ix-'^2x - x Jat ^ matrix A G 
]^/xj„^ is denoted by 3^ x„ ^, is an ( Ji x J2 x • • • J„_i x / x Jj+i • • • x 
JAr)-tensor of which the entries are given by (3^ x„ ^)jij2-- j«-iu«+i --jjv — 
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EJn 
j„ = 1 yjlj2 ■ ■ -in- ljnjn+1 ' ' 'JN ^ijn 



T5S 



As an example [ jl58j . let 3^ be the following 3x4x2 tensor: 



1 4 7 10 

2 5 8 11 

3 6 9 12 



/ 13 16 19 22 

Y,,2 = 14 17 20 23 

\ 15 18 21 24 

Let A be the following 2x3 matrix: 



A = 



1 2 3 \ 
^ 4 5 6 J ' 

Note that the number of columns in A is equal to the size of mode 1 of 3^. 
Thus we can compute 3^ Xi ^4, which is of size 2x4x2 and 



iyxiA),i = 
{yxiA),2 = 



' 22 49 76 103 ^ 
^ 28 64 100 136 y 

130 157 184 211 



^ 172 208 244 280 J 

(a) Given matrices A G B G R-^"^-^" 

y x„, A^nB = {y A) B = {y B) x,„ a {m^n) 

(b) liAeW^^\ 5gM^^^ 

yx^,Ax^B = yx^{BA). 

(c) If ^ G M^'''^" with full rank, then 

X = y x^A^y = X x^A^. 

The matrix unfolding and n-mode product are related via 
Proposition 1.5.2: |I5B] Let 3^ G R'^iX'^2x-x 
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(a) UAe W^-. Then 

X = y X^A^ X^^)= AY(^r^)■ 

(b) Consequently, if A^"-^ G j^^x j„ for all n G 1, 2, • • • , we have 

A' = 3^Xi^« X2^(2)--- X^^(^)^ 

Proposition 1.5.3: [158] Let A:', 3^ G M^ix^2x-x/a._ 

(a) Then ||A:'|| = , for n G {1, 2, • • • , A^}, where ||.||f is Frobe- 
nius norm of a matrix. 

(b) If A:' = a(i) o a^^) o • • • o a(^) and = b^^^ o 6(2) q . . . o b^^l Then 

(c) Let Q be a J X orthonormal matrix. Then ||A:'|| = x Q\\. 

In this thesis, we study detection and quantification of entanglement in 
multipartite quantum systems. The chapters are arranged as follows : 
Chapter 2 : In this chapter, we give a method to associate a graph with 
an arbitrary density matrix referred to a standard orthonormal basis in the 
Hilbert space of a finite dimensional quantum system. We study related 
issues such as classification of pure and mixed states. Von Neumann en- 
tropy, separability of multipartite quantum states and quantum operations 
in terms of the graphs associated with quantum states. In order to address 
the separability and entanglement questions using graphs, we introduce 
a modified tensor product of weighted graphs, and establish its algebraic 
properties. In particular, we show that Werner's definition (Werner [25] of 
a separable state can be written in terms of graphs, for the states in a real 
or complex Hilbert space. We generalize the separability criterion (degree 
criterion) due to Braunstein et al. [76] to a class of weighted graphs with 
real weights. We have given some criteria for the Laplacian associated with 
a weighted graph to be positive semidefinite. 

Chapter 3 : In this chapter, we settle the so-called degree conjecture 
for the separability of multipartite quantum states, which are normalized 
graph Laplacians, first given by Braunstein et al. [76]. The conjecture 
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states that a multipartite quantum state is separable if and only if the 
degree matrix of the graph associated with the state is equal to the degree 
matrix of the partial transpose of this graph. We call this statement to 
be the strong form of the conjecture. In its weak version, the conjecture 
requires only the necessity, that is, if the state is separable, the correspond- 
ing degree matrices match. We prove the strong form of the conjecture for 
pure multipartite quantum states, using the modified tensor product of 
graphs defined by Ah S. M. Hassan and P. S. Joag pTSQj . as both nec- 
essary and sufficient condition for separability. Based on this proof, we 
give a polynomial-time algorithm for completely factorizing any pure mul- 
tipartite quantum state. By polynomial-time algorithm we mean that the 
execution time of this algorithm increases as a polynomial in m, where m 
is the number of parts of the quantum system. We give a counter-example 
to show that the conjecture fails, in general, even in its weak form, for mul- 
tipartite mixed states. Finally, we prove this conjecture, in its weak form, 
for a class of multipartite mixed states, giving only a necessary condition 
for separability. 

Chapter 4 : In this chapter, we give a new separability criterion, a nec- 
essary condition for separability of A/^-partite quantum states. The crite- 
rion is based on the Bloch representation of a A-partite quantum state 
and makes use of multilinear algebra, in particular, the matrization of 
tensors. Our criterion applies to arbitrary A-partite quantum states in 
7i = H^^ ^ H^^ • • • H^^ . The criterion can test whether a A-partite 
state is entangled and can be applied to different partitions of the A-partite 
system. We provide examples that show the ability of this criterion to de- 
tect entanglement. We show that this criterion can detect bound entangled 
states. We prove a sufficiency condition for separability of a 3-partite state, 
straightforwardly generalizable to the case A > 3, under certain condition. 
We also give a necessary and sufficient condition for separability of a class 
of A-qubit states which includes A-qubit PPT states, 
chapter 5 : In this chapter we present a multipartite entanglement mea- 
sure for A-qubit pure states, using the norm of the correlation tensor which 
occurs in the Bloch representation of the state. We compute this measure 
for several important classes of A-qubit pure states such as GHZ states, W 
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states and their superpositions. We compute this measure for interesting 
applications like one dimensional Heisenberg antiferromagnet. We use this 
measure to follow the entanglement dynamics of Grover's algorithm. We 
prove that this measure possesses almost all the properties expected of a 
good entanglement measure, including monotonicity. Finally, we extend 
this measure to A^-qubit mixed states via convex roof construction and 
establish its various properties, including its monotonicity. We also intro- 
duce a related measure which has all properties of the above measure and 
is also additive. 

Chapter 6 : In this chapter, we summarize the work presented in this 
thesis and give the possible ways in which this work may be developed 
further. 



Chapter 2 



A combinatorial approach to 
multipartite quantum systems: basic 
formulation 



Before I came here I was confused about this subject. 
Having listened to your lecture I am still confused. But on a higher level. 

Enrico Fermi 



2.1 Introduction 



In chapter 1, we saw that quantum information is a rapidly expanding 
field of research because of its theoretical advances in fast algorithms, su- 
perdence quantum coding, quantum error correction, teleportation, cryp- 
tography and so forth [231 H^Oj E] • Most of these applications are based on 
entanglement in quantum states (see chapter 1). Although entanglement 
in pure state systems is relatively well understood, its understanding in 
the so called mixed quantum states [5], which are statistical mixtures of 
pure quantum states, is at a primitive level. Recently, Braunstein, Ghosh 
and Severini [761 EZj, have initiated a new approach towards the mixed 
state entanglement by associating graphs with density matrices and under- 
standing their classification using these graphs. Hildebrand, Mancini and 
Severini [TS] testified that the degree condition (see chapter 3) is equivalent 
to the PPT-criterion. They also considered the concurrence (see chapter 
1) of density matrices of graphs and pointed out that there are examples 
on four vertices whose concurrence is a rational number. In this chapter. 
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we generalize the work of these authors and give a method to associate a 
graph with the density matrix (real or complex), of an arbitrary density 
operator, and also to associate a graph with the matrix representing Her- 
mitian operator (observable) of the quantum system, both with respect to 
a standard orthonormal basis in Hilbert space. We define a modified tensor 
product of graphs and use it to give Werner's definition for the separabil- 
ity of a m-partite quantum system, in W'^ • • • W'"', as well as 
(g) . . . (g) QQm terms of graphs. We also deal with classification 
of pure and mixed states and related concepts like von-Neumann entropy 
in terms of graphs. 

The chapter is organized as follows. In Section 2.2, we define weighted 
graphs and their generalized Laplacians which correspond to density ma- 
trices, and discuss the permutation invariance of this association. We also 
deal with pure and mixed states in terms of graphs. Section 2.3 deals with 
von-Neumann entropy. Section 2.4 is concerned with separability issues 
as mentioned above. In Section 2.5, we deal with graph operations which 
correspond to quantum operations jT], 11611 1162j . In Section 2.6, we present 
a method to associate a graph with a general Hermitian matrix, having 
complex off-diagonal elements. We define the modified tensor product for 
complex weighted graphs and express the separability of mixed quantum 
states in a complex Hilbert space in terms of graphs, using Werner's defi- 
nition. In section 2.7, we present some graphical criteria for the associated 
Laplacian to be positive semidefinite. Finally, we close with a summary 
and some general comments. Sections 2.2 to 2.5 deal with graphs with 
real weights, that is, quantum states living in real Hilbert space. Graphs 
with complex weights, corresponding to density operators with complex off 
diagonal elements are treated in section 2.6. However, a large part of the 
results obtained for real Hilbert space in sections 2.2 to 2.5, go over to the 
case of complex Hilbert space (see section 2.8 (ix)). 
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2.2 Density matrix of a weighted graph 
2.2.1 Definitions 

A graph G = (V, E) is a pair of a nonempty and finite set called vertex 
set V{G), whose elements are called vertices and a set E{G) C V'^{G) of 
unordered pairs of vertices called edges. A loop is an edge of the form 
{vi^Vi] for some vertex Vi. A graph G is on n vertices if 1^(^)1 = n. We 
call the graph as defined above a simple graph. = m + s, where m 

is the number of edges joining vertices, s is the number of loops in G jl63l . 



A weighted graph (G, a) is a graph together with a weight function jl64 

a : V{G) X V{G) R 

which associates a real number (weight) a{{u^v}) to each pair {w, "u} of 
vertices. The function a satisfies the following properties: 

(i) a({w, v]) 7^ if {w, v) G E{G, a) and a({w, v]) = if {w, v) ^ E{G, a). 

(ii) a({w, v}) = a{{v^ u]) 

(iii) a{v^v) 7^ if {v^v} G E{G,a) and is zero otherwise. 

If e = {u^v} is an edge in E{G,a), property (ii) allows us to write a(e) or 
ttuv for a({w, v}). A simple graph can be viewed as a weighted graph with 
all nonzero weights equal to 1. 

In the case of simple graphs the degree dy of a vertex v G V{G) is defined 
as the number of edges in E{G) incident on v. For a weighted graph we 
set 



t'(G,a)(^) = f« = ^ CLuv (2.1) 
u€V{G,a) 

The adjacency matrix of a weighted graph with n vertices M{G, a) = 
[0'uv]u,v€V{G,a) is a n X n matrix whose rows and columns are indexed by 
vertices in V{G,a) and whose uv-th element is auv Obviously the adja- 
cency matrix M{G, a) is a real symmetric matrix with diagonal element vv 
equal to the weight of the loops on vertex v (i.e ayy). 
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The degree matrix for the weighted graph A(G,a) is a n x n diagonal 
matrix, whose rows and columns are labeled by vertices in V{G, a) and 
whose diagonal elements are the degrees of the corresponding vertices. 

A{G, a) = diag[dy; v e V{G, a)]. (2.2) 
The combinatorial Laplacian of a weighted graph is defined to be 

L(G, a) = A(G, a) - M{G, a). (2.3) 
The degree sum of (G, a) is defined as 

^{G,a)= ^v = TrA{G,a). (2.4) 

v&V{G,a) 

The Laplacian defined by Eq. (2.3) has no record of loops in the graph. 
Therefore we define the generalized Laplacian of a graph {G^a), which 
includes loops, as 

Q(G, a) = A(G, a) - M{G, a) + Ao(G, a) (2.5) 

where Ao(G', a) is a n x n diagonal matrix with diagonal elements equal to 
the weights of the loops on the corresponding vertices 

[Ao{G,a)]vy = ayy. (2.6) 

We call Ao(G', a) the loop matrix of the graph (G, a) . 

For a given weighted graph {G,a), the generalized Laplacian, defined 
by Eq. (2.5), is not necessarily a positive semidefinite matrix. When, for a 
given graph {G,a), the generalized Laplacian Q{G,a) is positive semidef- 
inite, we can define the density matrix of the corresponding graph (G, a) 
as 

a{G, a) = -^Q{G, a) = -^[L{G, a) + Ao(G, a)] (2.7) 

where Tr{a{G, a)) = 1. Note that, this definition of the density matrix of 
a weighted graph (G, a) reduces to that of the density matrix for a simple 
graph without loops [77] . 
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Whenever we can define the density matrix for a graph (G, a) we say 
that the graph (G, a) has density matrix. 

For any density matrix cr, we can obtain the corresponding graph as 
follows: 



(i) Determine the number of vertices of the graph from the size (n x n) of 
the density matrix. The number of vertices = n. Label the vertices 
from 1 to n. 

(ii) If the ij-th element of a is not zero draw an edge between vertices Vi 
and Vj with weight — Cij. 

(iii) Ensure that '0v^ = era by adding loop of appropriate weight to Vi if 
necessary. 

Example (2.1): For the following three matrices, we find the corresponding 
graphs. 



1 1 
1 1 



m 



Figure 2.1 
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Figure 2.3 



2.2.2 Invar iance under isomorphism 

Two weighted graphs (G, a) and (G', a') are isomorphic if there is a bijective 
map [1221 

(p : V{G,a) I — > V{G',a') 

such that 

e E{G'a') iff {v,, Vj} G E{G, a), j = 1, 2, • • • , n 

and 

^'(l){v^)(i){vj) = ^WiWj ^5 j = 1, 2, • • • , n. 
We denote isomorphism of (G, a) and {G',a') by (G, a) = {G'a'). 

Equivalently, two graphs (G, a) and (G', a') are isomorphic if there exists 
a permutation matrix P such that 

P'^M{G,a)P = M{G',a'). 

Note that, 

P^A(G, a)P = A(G', a'); ^^Ao(G, a)P = Ao(G', a') 
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Therefore we have 

P^Q(G,a)P = g(G',a'). (2.8) 

This means that Q{G^ a) and Q{G' ^ a') are similar and have the same eigen- 
values. Therefore, if Q{G,a) is positive semidefinite then so is Q{G',a'). 
Therefore, if (G, a) has the density matrix so does {G', a'). We have proved 

Theorem 2.2.1 : The set of all weighted graphs having density matrix is 
closed under isomorphism. ■ 

Since isomorphism is an equivalence relation, this set is partitioned by 
it, mutually isomorphic graphs forming the partition. 



2.2.3 Correspondence with quantum mechanics 

Henceforth, we consider only the graphs having density matrix unless 
stated otherwise. The basic correspondence with QM is defined by the 
density matrix of the graph. For a graph with n vertices the dimension of 
the Hilbert space of the corresponding quantum system is n. To establish 
the required correspondence we fix an orthonormal basis in the Hilbert 
space M*^^ <S) • • • <S) M*^"' of the system, which we call the standard basis 
and denote it by {\ijki ■ ■ i ■ ■ ■ = 1, 2, • • • , n = qiq2 ■ ■ ■ Qm^ or by 

{bi}}' i = ■ ■ ■ = qiq2 ■ ■ ■ Qm- We label n vertices of the graph (G, a) 
corresponding to the given density matrix by the n basis vectors in the 
standard basis. We say that the graph (G, a) corresponds to the quantum 
state (density operator) whose matrix in the standard basis is the given 
density matrix. Finally, we set up a procedure, by associating appropri- 
ate projection operators with edges and loops of (G, a) to reconstruct this 
quantum state from the graph {G,a). (See Theorem 2.2.7). In view of 
Theorem 2.2.1 , if {G^a) has density matrix a and {G,a) = [G'^a') with 
the corresponding permutation matrix P,then (G', a') has the density ma- 
trix P^aP. All of this paragraph applies to the complex weighted graph 
(section 2.6). 
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2.2.4 Pure and mixed states 

A density matrix p is said to be pure if Tr{p^) = 1 and mixed other- 
wise. Theorem 2.2.3 gives a necessary and sufficient condition on a graph 
(G, a) for cr(G, a) to be pure. For a graph [G^a) having k components 
(Gi, ai), (^2, 02), • • • , {Gk, cLk) we write (G, a) = (Gi, ai) W (^2, 02) W • • • W 
(Gfc, afc) where a^, i = 1, • • • , /c are the restrictions of the weight function 
of the graph (G, a) to the components. We can order the vertices such 
that M{G,a) = (Bi=iM{Gi,ai). When k = 1, {G,a) is said to be con- 
nected. From now on we denote by Xi{A), X2{A), • • • , Xk{A) the k different 
eigenvalues of the Hermitian matrix A in the nondecreasing order. The 
set of the eigenvalues of A together with their multiplicities is called the 
spectrum of A PE [IS51 [ffl] • 



Lemma 2.2.2 : The density matrix of a graph {G,a) without loops 
has zero eigenvalue with multiplicity greater than or equal to the num- 
ber of components of (G, a) with equality applying when weight function 
a =constant > 0, 

Proof : Let (G, a) be a graph with n vertices and m edges. Since Q(G, a) 
is positive semidefinite, for x G iR" we must have |164] 



X^Q{G, a)x — ^ CLikjki^ik ~ ^jk)"^ + ^ (^itk^lt — ^■ 

k=i t=i 
For the graph without loops the above inequality becomes 

x^Q{G,a)x = - x^J^ > 0. (2.9) 

For x^ = (1 1 • • • 1) we can see x^Qx = 0. This means that x^ = 
(1 1 1 • • • 1) is an unnormalized eigenvector belonging to the eigenvalue 
|129j . If there are two components (Gi,a) and (G2,a) of (G,a), with 



ni,mi and 712, m2 being the numbers of vertices and edges in (Gi,a) and 
(G2,a) respectively, we can decompose the sum in Eq.(2.9) as 



mi TO2 



x^Q(G, a)x = ^ a,^^j,^{x,^^ - Xj,f + ^ a,^^j,^{x,^^ - XjJ\ (2.10) 

ki = l k2=l 
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For = (1 1 1 • • • 1) both the terms in Eq. (2.10) vanish. Now consider 
two vectors = (0 • • • 01 1 • • • 1) with first rii components zero and last 
712 components 1 and = (1 1 • • • 1 • • • 0) with first ni components 
1 and last 77.2 components zero, (ni + 77-2 = n). Obviously the RHS of Eq. 
(2.10) vanishes for both xi and X2. This implies xi and X2 are two orthog- 
onal eigenvectors with eigenvalue zero. This means that the multiplicity 
of zero eigen value is at least 2 (number of components in (G, a)). 

The equality condition for =constant> V {w, v} G E{G, a) is 
proved in [77] . ■ 



Theorem 2.2.3 : The necessary and sufficient condition for the state 
given by a graph (G, a) to be pure is 

n m 

E + 2 E4« = (2-11) 

i=l fc=l 

where is the degree of the vertex vi, ai^j^ is the weight of the edge 

n 

{^*fe5 7^ and t)(G,a) is the degree sum t)(G,a) = ^i- 

i=l 

Proof : Equation (2.11) is just the restatement of the requirement Tr{a'^{G, a)) 
= 1, which is the necessary and sufficient condition for the state a{G,a) 
to be pure. ■ 
Lemma 2.2.4 : The graph {G,a) for a pure state a{G,a) has the form 
{Ki, h) 1+1 v^j^i 1+1 viJ^2 l±l • • • 1+1 1;^ for some 1 < f < n. 
Proof : Since the state is pure, it has the form 



= X]c^Ji;^,),l <ik<n. 



k=l 



We can permute the basis vectors to transform this sum to = qI-u, 

That is, the ^ basis kets contributing to the sum in the above equation 
become the vectors l-ui), |'U2), • • • , 1"^^) under this permutation. The result- 
ing density matrix has a block of first ^ rows and first ^ columns 
all of whose elements are nonzero, while all the other elements of density 
matrix are zero. The graph corresponding to this density matrix is just 
the required graph. ■ 
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Example (2.2) : We now give important cases of pure state graphs in 
which we use later. 



Figure 2.4 
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" 1 -1 " 


2ai2 


-CLl2 


ai2 
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-1 1 



(i) aiK2,a) = ^^ 
the corresponding graph is as shown in Figure 2.4. 



Figure 2.5 
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graph is as shown in Figure 2.5. 



= P[|i;i)], the corresponding 



2a 



0-^ 



2a 



Figure 2.6 



(iii) ai2 > 0, (j{K2, -a) = ^ 



2ai2 



au ai2 
ai2 ai2 



|i'2))], a > 0. The corresponding graph is as shown in Figure 2.6 



1 1 
1 1 



11 



12 



m 
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Figure 2.7 
(iv) 



(t{G, a) = j 



1 

-1 
-1 
1 



-1 1 

1 -1 

1 -1 

-1 1 



= ^[(|-)|-))]: 
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where |— ) = -^(ll) ~ 1 2)), the corresponding graph is as shown in Figure 
2.7. 




4 



Figure 2.8 
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= ^[(l+)l+))], 



where |+) = ^^^dl) + |2)) in the corresponding graph is as shown 

in Figure 2.8. 

It may be seen that in each of the cases in example (2.2), same density 
matrix on the standard basis corresponds to infinite family of graphs as 
the nonzero weight on each edge or loop is multiplied by a constant. But 
this is a false alarm because any weight a 7^ 1 only changes the length 
of the corresponding state vector in the Hilbert space (i.e. state becomes 
unnormalized) which does not have any physical significance. Another 
example pertaining to this situation is the random mixture (see lemma 
(2.3.1 )). 



(j(G', a) = — 
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an 



a 











a 
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n 



n 







2. A combinatorial approach to multipartite quantum systems: basic 
44 formulation 



However, this does not lead to any contradiction because of the uniqueness 
of the random mixture [5j . 

All the density matrices in (i), (ii), (iii), (iv), (v) above represent pure 
states. 

Remark 2.2.5 : Any graph with the weight function a = constant > 
has the same density matrix for all a > 0. This infinite family of graphs 
corresponds to the same quantum state (density operator). 
Definition 2.2.6 : A graph {H, b) is said to be a factor of graph (G, a) if 
V{H, b) = V{G, a) and there exists a graph {H', b') such that V{H', b') = 
V{G,a) and M{G,a) = M{H,b) + M{H\b'). Thus a factor is only a 
spanning subgraph. Note that 



by^y^ if {v,,Vj} e E{H, b) 



Now let (G, a) be a graph on n vertices "Ui, • • • ^Vn having m edges 
{^^i.^jjr-- and s loops {v^^,v,^} ■ ■ ■ {v,^,v^} where 1 < ii, ji, 

■■■,imjm<n, I < 11,12,- ■ ■ ,^s < n. 

Let {Hij^jj,, ai^.j^.) be the factor of {G,a) such that 



[M{H,,j„a,,j,)]u,y,- I otherwise 



a^^j^ if u = ik and w = jt or u = jk,w = ik 

(2.12: 



Let {Hi^^i^, a^i^) be a factor of {G,a) such that 

lAj/u M f ^^ttt when u = it = w , . 

[M(i?„,,o,,)]„„ = |j, ^^^^^.^^ (2.13) 

Theorem 2.2.7 : The density matrix of a graph {G,a) as defined above 
with factors given by equations (2.12) and (2.13) can be decomposed as 

^ m 1 * 

a{G, a) = V 2a,,j,a{H,^j„ a^.^J + V a^^^^a{H,^^^, a^^^J 

or 

1 ™ 1 1 
a{G,a) = V2a,,^,P[— =(|z;,,) - \vj,))] + V a,,,,P[|i;,, 

f (G,a) V2 f (G,a) ^ 
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Proof : From equations (2.12), (2.13) and theorem 2.2.3 and lemma 2.2.4, 
the density matrix 



2a 



ikjk 



is a pure state. Also, 

cr(i^Mt,aMt) = - — [^o(^zt,it,aMj] 
is a pure state. Now 

m s 

A(G,a) = ^A(i/,,,, 

5 ^itit) 



k=l 



t=l 



5 ^itit)- 



k=l 



t=l 



Therefore 



(7(G, a) = 



m 



^i^ikjk^ ^ikjk) - X] ^{Hikjk^ (^ikjk) 



.k=l 



k=l 



'(G,a) 



= /S^jHikjk, CLikjk) - MiHikjk^ (^ikjk)] + Z /_] ^o(^Mt, a^^ij 

^(G,a) ^ ^{G,a) ^ 

In terms of the standard basis, the wit;-th element of matrices a{Hij^jj^,ai^j^] 
and a{Hi^^^, a^^i^) are given by {vu\ct{H^^j^, , a^^jj\vyj) and {vu\a{H^^^^, a^^^^)\vyJ) 
respectively. In this basis 
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Therefore equation (2.14) becomes 

m s 

a(G,a) = Y.^a,,,,P[-^{\v,,)-\v,,))+- J] a,,,,P[|^,J] (2.15) 

■ 

Remark 2.2.8 : If all weights ai^^j^^ > then equations (2.14), (2.15) give 
a{G,a) as a mixture of pure states. However, in the next subsection we 
show that any graph (G, a) having a density matrix can be decomposed 
into graphs (spanning subgraphs) corresponding to pure states. 

2.2.5 Convex combination of density matrices 

Consider two graphs (Gi,ai) and {02,0,2) each on the same n vertices, 
having (j{Gi,ai) and g{G2,cl2) as their density matrices, respectively. We 
give an algorithm to construct the graph (G, a) whose density matrix is 

(7(G, a) = A(7(Gi, ai) + (1 — \)a[G2-, (12) 

with 0<A< 1,A = q;//?, q;,/?>0 being real. 

We use the symbol U to denote the union of the edge sets of two graphs 
(Gi, ai) and (G2, 02) on the same set of vertices to give (G, a). If {-Ui, Vj] G 
E(Gi,ai)and {vi,Vj} G E{G2,a2) then a{{v^,Vj}) = ai{{v^,Vj})+a2{{v^,Vj}). 
We write (G, a) = (Gi,ai) U (G2,a2)- If E(Gi,ai) and E(G2,a2) are dis- 
joint sets, then we call the resulting graph (G, a) the disjoint edge union 
of (Gi, ai) and (G2, 02), we write (G, a) = (Gi, ai) + (G2, a2)- 

The algorithm is as follows : 
Algorithm 2.2.9 : 

(i) Put X = a/f3 so that (1 — A) = where a > 0, (3 > are real, 
(h) Write (7(G, a) = ■^(q;(7(Gi, ai) + (/? — q;)(7(G2, 02))- 

(iii) Modify the weight functions of the two graphs (Gi,ai) and (G2,a2) 
to get a'l = aai and 02 = (/3 — a)a2. 
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(iv) The graph (G, a) corresponding to a in step (ii) is 

(G',a) = (G'i,a;)U(G'2,ay (2.16) 

such that 

aviVj = {ai)v,v, + {a2)v^vj (2.16a) 

aviVi = MviVi + ia2)viVi (2.165) 

where we take ia[2)vivj = = (^'12)^;^ if {vi.vj}, {vi.Vi} ^ E{Gi,ai) 

or E{G2, a2) 



We can apply this algorithm to any convex combination of more than 

k 

two density matrices a{G,a) = ^pia{Gi,ai), = 1, by writing pi = 

i=l i 

ai/P, ai, P > and real, i = 1, - ■ ■ ,k. 
Example (2. 3) : consider the density matrices 



1111 
1111 
1111 
1111 
whose graph is shown in Figure 2.9 



(i) a(Gi,ai) = ! + +)(+ + 




Figure 2.9 



CT(G2,a2) = 



l|ll)(ll| + l|^+)(t^+| 



1 
4 



2 

110 

110 
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whose graph is shown in Figure 2.10 




21 22 



Figure 2.10 



The graph corresponding to 



1 2 

a{G,a) = -cr(Gi,ai) + -cr(G2,a2) 



1 

12 



5 111 
13 3 1 
13 3 1 
1111 



is given in Figure 2.11 




Figure 2.11 



Lemma 2.2.10 : Let (Gi,ai), (02,02) and (G, a) satisfy 



(G, a) = (d, ai) U (6*2, 02) 



or, 



Then 



(G, a) = {Gi, ai) + (G2, 02)- 



g(G,a) = Q(Gi,ai) + g(G2,a2) 
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and 



a) = ^^^(7(Gi, ai) + 5^^(7(G2, ^2). 



Proof : For two factors of {G,a), (Gi,ai) and {G2,cl2) we have 

M{G,a) = M{Gi,ai) + M{G2,a2) 

A{G,a) = A{Gi,ai) + A{G2,a2) 

Ao(G,a) = Ao(Gi,ai) + Ao(G2,a2) 

L{G,a) = A{G,a)-M(G,a) 

Q{G,a) = L{G,a) + Ao{G,a) 

Substitute M{G,a),A{G,a),Ao{G,a) and L{G,a) in Q{G,a) as above to 
get 

g(G,a) = g(Gi,ai) + g(G2,a2) 

and also 

^(G, a) = ^^l^aiGu ai) + ^^a(G2, a2). 

0(G,a) 0(G,a) 

■ 

Remark 2.2.11 : Obviously, the operation U is associative. We can apply 
lemma 2.2.10 for more than two graphs, 

(G, a) = Ui{Gi, ai) Q{G, a) = ^ QiGi, ai) 

i 

and 



(G, a) = ^ d{Gi, ai)(T{Gi, a^). 



a 

Theorem 2.2.12 : Every graph {G,a) having a density matrix a{G,a) 
can be decomposed as {G,a) = L\i{Gi,ai) where a{Gi,ai) is a pure state. 
Proof : Every density matrix can be written as the convex combination 

k 

of pure states (j{G,a) = ^Pi\ipi){ipi\- 

i=i 

By applying algorithm 2.2.9, lemma 2.2.10 and remark 2.2.11, we get 
the result. ■ 



2. A combinatorial approach to multipcirtite quantum systems: basic 
50 formulation 



2.2.6 Tracing out a part 

Consider a bipartite system with dimension pq. Let (t{G^ a) be a state of 
the system with graph (G, a) having jgg vertices labeled by (ij), i = 1, • • • ,j9 
and J = 1, • • • , If we trace out the second part with dimension we 
get the state of the first part which is p x p reduced density matrix of 
cr(G, a). The corresponding graph {G\a') has p vertices indexed by (z) 
and its weight function a' is given by 

fc=l 

and 

q 

k=l l&V{G',a') 

where Z)^^ is the degree of vertex {ik) in original graph. 

Example (2.2.4) '• Consider a graph [G^a) as shown in Figure 2.12a in 
M^. The corresponding density matrix is 

9 -1-1 1 " 
-1 3 -1 -1 
-1 -13-1' 
1 -1-1 1 _ 

After tracing out the second particle the graph {G', a') on two vertices 
becomes as in Figure 2.12b with corresponding density matrix 



a^{G',a') = 



1 


"12 -2 " 


1 


"6 -1 ■ 


16 


-2 4 


~ 8 


-1 2 



1 



which is the same as the reduced density matrix a of a 
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(b) 



2.3 Von Neumann entropy 

The Von Neumann entropy of the n x n density matrix a is 

n 
i=l 

It is conventional to define OlogO = 0. The Von Neumann entropy is a 
measure of mixedness of the density matrix. For a pure state a, S{(j) = 0. 

2.3.1 Maximum and minimum 

Let 

(G,a) = Wti(A'l,a,) (2.17) 

where a^) is the graph on i-th vertex with a loop having weight > 0. 
Lemma 2.3.1 : Let (G, a) be given by Eq.(2.17) with the additional 
constraint that = c = = 1,2, ••• ,n. The density matrix of the 
graph (C, a) is the random mixture of pure states with cr(G, a) = ^In- 
Proof : For the graph (G, a), the first term in Eq.(2.14) vanishes. Then 

1 " 

a{G, a) = V Ao{Hi^i^, a) 

f (G,a) f-f 

where AQ{Hi^i^,a) is the n x n matrix with all elements zero except the 
{it, it)th element which is equal to a. This means 

o" 



1 




(t{G, a) = 



a 



\G,a) 
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because d(^G.o.) = i^^^- ' 
Theorem 2.3.2 : Let (G, a) be a graph on n vertices. Then 

(i) max^G,a)S{a{G,a)) = log2n 

(ii) min(Gf a) S{(7{G, a)) = 0, and this value is attained if (7{G, a) is pure. 
Proof : 

(i) By lemma 2.3.1 a{G,a) defined in the lemma has eigenvalue 1/n 
with multiplicity n. The corresponding Von Neumann entropy is log2 n. 
Since {G^a) is on n vertices, the support of a{G^a) has dimension < n. 
Any matrix having dimension of support < n cannot have Von Neumann 
entropy > log2 n. 

(ii) For a pure state S{a{G,a)) = and S{a{G,a)) ^0. ■ 

2.4 Separability 

In this section we primarily deal with the graphs representing a bipartite 
quantum system with Hilbert space of dimension pq. Obviously, 
the corresponding graph has n = pq vertices. We label the vertices using 
standard (product) basis {\vi) = \us+i) <S) \wt)}, < s < p — 1, 1 < t < 
q,i = sq + t. 

2.4.1 Tensor product of weighted graphs 

The tensor product of two graphs (G, a) and {H, b) denoted (G, a) <S> {H, b) 
is defined as follows. 

The vertex set of (G,a) (g) {H,b) is V{G,a) x V{H,b). Two vertices 
{ui^vi) and {u2,V2) are adjacent if {^1,^2} £ E{G,a) and {vi,V2} G 
E{H,b). The weight of the edge {{ui,vi), (^2,^2)} given by a{u,,u2}b{vi,v2} 
and is denoted by c({(wi,'Ui), {u2,V2)}). Note that either ui and U2 or vi 
and V2 or both can be identical, to include loops. 

The adjacency, degree and the loops matrices of (G, a) <S> {H, b) are given 

by 

M((G, a) {H, b)) = M{G, a) M{H, b) (2.18a) 
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A((G, a) (g) {H, b)) = A{G, a) ^ A{H, b) (2.186) 
Ao((G, a) {H, b)) = Ao{G, a) ® Ao{H, b) (2.18c) 

Note that 

L((G, a) ® {H, b)) ^ L{G, a) ® L{H, b) 

Q{{G, a) {H, b)) ^ Q{G, a) Q{H, b). 

In fact, in general, the tensor product of two graphs having a density matrix 
may not have a density matrix. 

For two simple graphs G and H we know that |167l , [77] 



This result is also satisfied by the tensor product of the weighted graphs. 

^'(G,a)®(i/,6) = '^{G,a) ' ^{H,b)- (2.19) 

2.4.2 Modified tensor product 

We modify the tensor product of graphs in order to preserve positivity of 
the generalized Laplacian of the resulting graph. 
Given a graph (G, a) we define {G^, a) by 

V{G''',a) = V{G,a) 

E[{G^, a)] = E{G, a) \ {{v,, v,] : [v,, v,} G E{G, a)} 

That is, (G'^, a) is obtained from (G, a) by removing all loops. 
Given a graph (G, a) we define (G, a) by 

V{G,a) = V{G, a) 

E{G, a) = E{G, a) \ {{v„ v^} : i ^ j, {v„ v^} G E{G, a)}. 

That is, (G, a) is obtained by removing all edges connecting neighbors and 
keeping loops. 

Note that in both (G*^, a) and (G, a), the weight function a remains the 
same, only its support is restricted. 
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Given a graph (G, a) we define (— G, a) = (G, —a). Given a graph (G, a) 
we define (G^, a') 

y(G^,a') = y(G,a) 

(G# a')=wr(^i,a;), 

where Ki is the graph consisting of ith vertex with a loop and a'^ is the 
weight of the loop on the ith vertex. If = then there is no loop on the 
ith vertex. a^,i = 1,2, ••• n are given by 

VkGV{G,a) 

Note that the term Vk = Vi is also included in the sum. 
We now define the graph operators on the set of graphs 



(ii) M{C{G,a)) = M{G,a) - Ao{G,a) 
A{C{G, a)) = A(G, a) - Ao(G, a) 

Ao(/:(G,a)) = [0] 
MG,a) = Tr{A{G, a)) - Tr{Ao{G, a)) = d^o^a) 

(iii) M{J\f{G,a)) = A(G,a) 
A{Af{G,a)) = A{G,a) 
Ao{M{G,a)) = A{G,a) 
dj^^G,a) = TriAiG,a)) = i>^G,a) 



(2.206) 



(i) r] : (G, a) — ^ (— G, a) = (G, —a) 

(ii) /::(G,a)^(G'^,a) 

(iii) TV : (G, a) ^ (G#, a') 

(iv) 1^ : (G, a) (G, a) 

Some properties of the graph operators defined in (2.20b) are, 

(i) M{r]{G,a)) = -M{G,a) 

A{r]{G,a)) = -A{G,a) (2.21) 
Ao(77(G,a)) = -Ao(G,a) 



(2.22) 



(2.23) 
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(b) 
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(c) 



(e) 



Figure 2.13 



(iv) M{n{G,a)) = Ao(G,a) 
A{n{G,a)) = Ao{G,a) 
Ao(0(G,a)) = Ao(G',a) 
fo(G,a)=Tr(Ao(G,a)) 



(2.24) 



Example (2.5) : Given a graph [G^a) as shown in figure 2.13a, if we act 
by ?7, jC ^Af and Q on {G,a), we get the graphs r]{G^a), jC{G,a), J\f{G^a) 
and a), as shown in figures 2.13(b), (c), (d) and (e), respectively. 
Definition 2.4.1 : Let {G, a) and {H, b) be two graphs with j9 and q (> p) 
vertices, respectively. Then their modified tensor product is defined by 



iG,a)B{H,b) = {C{G,a)0Cr]{H,b)} + {jC{G,a)0M{H,b)} 

+{Ar(G, a) (g) C{H, b)} + {^7(G, a) n{H, 6)}, (2.25) 



whose cardinality is pq. 

E{{G, a) □(if, b)} = disjoint union of the edge set of each term in (2.25). 
Lemma 2.4.2 : 



Proof : Consider the degree matrix of the modified tensor product. We 



and 



V{{G, a) □ {H, b)} = V{G, a) x V{H, b) 



(i) A((G, a) □ {H, b)) = A{G, a) A{H, b). 

(ii) Ao((G, a) □ (F, 6)) = Ao(G, a) ® Ao{H, b). 
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have 

A((G, a) = A(/:(G,a)0£r/(//,6)) + A(£(G,a)(g)Ar(F,5)) 

+A(A/'(G, a) C{H, b)) + A(Q(G, a) 6)). 

This follows from lemma 2.2.10. Using equation (2.18b) and equations 
(2.21) to (2.24) to the terms on the RHS of the above equation we get 

A((G, a) □ {H, b)) = A(G, a) (g) A{H, b). 

Equation (ii) can be proved similarly. ■ 
Corollary 2.4.3 : ^{G,a)B{H,b)ivhV2) = f (G,a)(^i) • ^ (f,6)(^2) 
Proof : This follows directly from equation (i) in lemma 2.4.2. ■ 
Remark 2.4.4 : From corollary we get t)(G,a)H(i?,&) = ^(Ca) ■ ^{H,b) 
Theorem 2.4.5 : Consider a bipartite system in MP (8) in the state a. 
Then cr = cti (8) (72 if and only if a is the density matrix of the graph (G, a) □ 
{H, b) , where {G, a) and (H, b) are the graphs having density matrices ai 
and (72, respectively. 

Proof : If part : Given (G, a), (H, b) we want to prove 
(t{{G, a) □ {H, b)) = (Ti{G, a) (g) (T2{H, b). 
From the definition of the modified tensor product we can write 

(7((G, a) □ [H, b)) = {Q[CiG, a) (g) CriiH, b) 

^{G,a)E]{H,b) 

+£(G, a) X{H, b) + M{G, a) C{H, b) + ^(G, a) Q(F, 6)]} 
Using lemma 2.2.10, remark 2.2.11 and remark 2.4.4 we get 

a((G,a)H(i/,6)) = \ [Q[C{G,a)®Cr^{H,b)) 

+Q{jC{G, a) M{H, b)) + Q{M{G, a) 

^C{H, b)) + g(Q(G, a) Q{H, b))]. (2.26) 
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We can calculate every term in (2.26) using (2.21) - (2.24) and substitute 
in (2.26) to get 

a{{G, a) □ (if, 6)) = a{G, a) ® a{H, b). 

Only if part : Given a = ai<S> cr2, consider the graphs (G, a) and (i7, b) for 
(7i and (72 respectively. Then the graph of a has the generalized Laplacian 

[L(G, a) + Ao(G', a)] (g> [L{H, b) + Ao{H, b)] 

= L{G, a) (g) L{H, b) + L{G, a) Ao(G, a) 

+Ao{G, a) (g) L{H, b) + Ao(G, a) (g) Ao(^, 6). 

Now it is straightforward to check that the graphs corresponding to each 
term are given by the corresponding terms in the definition of (G, a) □ 
{H,b). ■ 
Remark 2.4.6 : Note that the proof of Theorem 2.4.5 does not depend in 
any way on the positivity or the Hermiticity of the associated generalized 
Laplacians. Therefore we have 

Q((G, a) □(//,&)) = Q{G,a)®Q{H,b) 

for any two graphs (G, a) and (if, b) 

Corollary 2.4.7 : The modified tensor product is associative and dis- 
tributive with respect to the disjoint edge union +. 

Proof : Let (Gi,ai), (^2,02) and (G3,a3) be any graphs. Using theorem 
2.4.5 and remark 2.4.6, we can write 

Q(((Gi,ai)H(G2,a2))H(G3,a3)) 

= g((Gi, ai) □ (G2, a2)) g(G3, 03) 
= (Q(Gi,ai)(g)g(G2,a2))(g)Q(G3,a3) 

= g(Gi,ai)(8)(Q(G2,a2)(g)Q(G3,a3)) 
= g(Gi,ai)®Q((G2,a2)H(G3,a3)) 
= g((Gi,ai)H((G2,a2)H(G3,a3)) 



Therefore, 

((Gi,ai)H(G2,a2))H(G3,a3) = (Gi, ai) □ ((G2, 02) □ (G3, 03)). 
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Similarly, using lemma 2.2.10 and distributive property of the matrix 
tensor product we get 

Q{{Gi,ai)B{{G2,a2)) + {Gs,as))) 
= g((Gi, ai) □ (G2, a2)) + ((d, ai) □ (G3, as)), 

which gives 

(Gi,ai)H((G2,a2)) + (G3,a3)) = (Gi, ai) □ (G2, 02) + (Gi, ai) □ (^3, 03). 

■ 

Definition 2.4.8 : The Cartesian product of two weighted graphs {G, a) 
and {H,b) is denoted {G,a)n{H,b) with weight function c defined as fol- 
lows: 

V{G,a) X V{H,b). 

E{{G,a)\J{H,b)) = {{{%v),{x,y)}\ u = x and {v,y} G E{H,b), v ^ 
y, c({(w, v), (w, y)}) = du ■ b{{v, y]) 01 v = y and {u, x] G E{G, a), u ^ 
X, c{{{u, v), {x, v)}) = dy ■ a{{u, x}), 

where X^u and Z)„ are the degrees of the vertices u G V{G, a) and v G 
V{H,b), respectively. It is straightforward to check that 

{G,a)n{H,b) = C{G,a) ^ Af{H,b) ^ Af{G,a) ^ C{H,b), 

which can be taken to be the definition of the Cartesian product of 
graphs in terms of the operators C and J\f . We also note that 

(G, a)B{H, b) = {CiG, a)®CT]{H, 6)}+{(G, a)U{H, 6)}+{^^(G, a)m{H, b)}. 

Note that the isolated vertices in (G, a) or (i^, b) do not contribute to 
(G, a)n(i7, b) as their degree is zero. 

Example(2.6) : Consider {G,a),{H,b) where V{G,a) = {1,2}, E{G, a) = 
{{1,2}} and ViH,b) = {1, 2, 3, 4}, 6) = {{1, 2}, {2, 3}, {3, 4}} with 
weight functions a = 6 = 1, as shown in figure 2.14a, b. The modified 
tensor product of these graphs is given by figures 2.15a, b, c, d, for each 
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term in (2.25), and the resulting graph is as shown in figure 2.15e. The 
corresponding density matrix of the graph (G, a) □ {H, b) is 



aiiG,a)BiH,b)) = ^ 



1-10 0-11 

-12-10 1-2 1 

0- 12-10 1-21 
-1 1 1 -1 

-110 1-10 

1- 2 1 0-12-10 
1-2 1 0-12-1 
1 -1 -1 1 



which is the same as (j{G, a) (8) cr{H, b). 



(a) 



4 ' 3 
(b) 



Figure 2.14 



Corollary 2.4.9 : The density matrix of the modified tensor product of 
two graphs is separable. 

Proof : From theorem 2.4.5, we see that cr((G, a) □ iH,b)) is actually a 
product state. ■ 
Corollary 2.4.10 : a = ai<S)(J2'S)- ■ •'SxJk for a /c-partite system if and only 
if the graph of a is the modified tensor product of the graphs of di, • • • , dfe. 
Proof : Apply theorem 2.4.2 successively to (ai (12), {{ai (12) <S) (T3) • • • 
and then use the associativity of the modified tensor product corollary 
2.4.7. ■ 
Corollary 2.4.11 : A state cr of a /c-partite system is separable if and 
only if the graph (G,a) for a has the form 



(G,a) = U, 0^^,(0], a]) 
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Proof : Let a be separable, i.e., 

i i 

By algorithm 2.2.9 and corollary 2.4.10 the graph of a has the form 

(G,a) = U, H^^i (Glai). 

Now let the graph of a /c-partite state be 

(G,a) = U, (Gi^ai). 

Then by lemma 2.2.10, remark 2.2.11 and the above corollary to theorem 
2.4.5 

a{G, a)=Y^ w,a[^^ ^ • • • ^. 

i 

■ 

Corolary 2.4.11 says that Werner's definition [25] of a separable state in 
the W^^ (g) . . . (g)M* system can be expressed using corresponding 

graphs. 
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Lemma 2.4.12 : For any n = pq the density matrix ai^Kn-, a) is separable 
in MP <S) if the weight function is constant > 0. 

Proof : The proof is same as that given for corollary 2.4.3 in [77j, for 
simple graph. ■ 
Example (2.7) : Consider the graph [K^^a). The vertices of {K^^a) are 
denoted by 1, 2, 3, 4, where weight function is constant, say , a = 3 > and 
has loops in vertices 1,2. We associate to these vertices the orthonormal 
basis {|1) = |1)|1), |2) = |1)|2), |3) = |2)|1), |4) = |2)|2)}. In terms of this 
basis (j(K4, a) can be written as 



(j(K4, a) = 



1 

42 



12 -3 
-3 12 
-3 -3 
-3 -3 



-3 -3 

-3 -3 

9 -3 

-3 9 



1 

14 



4 
-1 
-1 
-1 



-1 -1 
4 -1 
-1 3 
-1 -1 



-1 
-1 
-1 
3 



and from equation (2.15) we can write a{K4,a) as 



<^(A-4.«) = ^{6P[|1)-J=(|1> - |2»] + 6P[^(|1) - |2))|1 

+6P[-J=(|1> - |2»|2>] + 6P[|2>-J=(|1> - |2))] + 6P[;^(|11) - |22»] 
+6P[-^(|12) - |21»] + 3P[|11)] + 3P[|12)]} 



<x(A'4,a) = ip[|l)-L(|i)_|2))] + ip[-L(|l)_|2))|l 

+ip[i=(|l) - |2»|2>] + ^^^[12)^(11) - |2»] 
+^{5^175(111) - |22»] + 5^175(112) - |21»]} 

+j-iPl\n)] + ^^Pl\i2)]. 

Each of the first four terms in the above expression is a projector on a 
product state, and also the last two terms are projectors, while the fifth 
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and sixth terms give rise to the separable density matrix ^p[\ —) \ +)] + 

where |±) =^ ^(|1) ± |2)) [77|. Thus a{K4,a), a constant, is 
separable in M^. 

Note that there exists a graph which is complete with a real weight 
function, which is entangled as the following graph shows in figure 2.16 



11 J 12 




21 ^ 22 



Figure 2.16 

Remark 2.4.13 : The separability of a{Kn, a) with constant weight func- 
tion > does not depend upon the labeling of V{Kji^ a) provided every 
vertex has a loop or there are no loops. Given a graph, an isomorphism 
from (G, a) i — > (G, a) is called automorphism. Under composition of 
maps, the set of automorphisms of (G, a) form a group denoted Aut{G, a). 
If a{Kn,a) is separable, and if the Aut{Kn,a) = Sn, [G^a) = {Kn^a) is 
also separable. Note that Aut{Kn, a) = Sn provided all weights are equal 
and either every vertex has a loop or there are no loops. 
Lemma 2.4.14 : The complete graph {Kn,a) on n > 2 vertices corre- 
sponding to a separable state with weight function constant > is not a 
modified tensor product of two graphs. 

Proof : It is clear that if n is prime then {Kn: a) is not a modified tensor 
product of graphs. We then assume that n is not a prime. Suppose that 
there exist graphs (G, b) and {H, c) respectively on p and s vertices such 
that {Kps,a) = (G, 6) □ {H,c), where b and c are constants. From the 
definition of the modified tensor product 

a({(wi, vi), (w2, V2)}) = b{{ui, U2}) ■ c{{vi, V2}), 

the degree sum is 

^iG,b)= Y b^v = 2b\E{G,b)\. 

ueV{G,b) ueV{G,b) weV{G,b) 
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We know that d(^G,b) < Kpi.P ~ 1)) also f (/f,c) = 2c\E{H, b)\ < cs{s — 1) 
and 

f(G,6) • ^{H,c) < bcps{p - l)(s - 1) = bcps{ps-p- s + 1). (2.27) 

Now observe that V{{G, b) □ {H, c)) = ps and, 

^'(G,6)H(ir,c) = a]3sbs - 1) (2.28) 

because (C, 6) □ (i7, c) = {Kps, a). 
We know that 

^{G,b)B{H,c) = ^{G,b) ■ ^{H,c)- (2.29) 

Substituting from Eq.(2.27) and Eq.(2.28) we see that Eq.(2.29) is satisfied 
only when p = 1 = s , i.e. n = ps = 1. ■ 

Lemma 2.4.12, lemma 2.4.14 and theorem 2.4.5 together imply that a 
complete graph (K^, a) on n > 2 vertices with a = constant > is a 
separable state but not a product state. 

Definition 2.4.15 : Consider a graph {G,a), without loops, pertaining 
to a bipartite system of dimension pq . The partial transpose of {G^a)^ 
denoted {G^^,a'), is a graph defined as V{G^^,a') = V{G,a), {il,kj} G 
E{G^^,a') ^ {ij, kl) e E{G, a) and a'{{il, kj}) = a{{ij, kl}). 
Lemma 2.4.16 : Consider a bipartite separable state a{G,a) with the 
associated graph (G, a) without loops. Then A(G, a) = A(G^^,a'), where 
{G^^,a') is the partial transpose of (G, a). 

Proof : Let Q{G, a) be the Laplacian of a graph (G, a) with real weights 
without loops, on n vertices. Let D be any n x n real diagonal matrix in 
the standard orthonormal basis {If^)}; i = 1, 2, . . . , n, such that D and 
Tr{D) = 0. This means that there is at least one negative entry in the 
diagonal of D. Denote this element by Da = 6^. Let iV^o) = X^j 1"^^^) and 
10) = J2jXj\vj) where 

r ifj^z 

\keRifj = i 
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Let \x) = 1^) + 10) = Elfi^Xj)\vj). Then 

(x|Q(G,a) + D|x) = (V^o|Q(G',a)|^o) + (V^o|Q(G,a)|0) + (0|Q(G,a)|V^o) + 

(0|Q(G, a)|0) + {iJo\D\i/jo) + (V^o|i^|0) + (0|i?|^o) + 

Since l^/^o) is (unnormalized) vector having all components equal unity, 
from equation (2.9) it follows that {ilJo\Q{G,a)\ilJo) = 0. Also (?/^o|^|V^o) = 
Tr{D) = 0. We have 

{cJ)\Q{G,aM = k\Q{G,a)),, = k% 

(^o|Q(G,a)|0) = (0|Q(G,a)|V^o)=O. 
Finally, the remaining terms in the above equation are given by 

{MD\ci)) = b,k = imiJo). 

Thus 

{X\Q{G, a) + D\x) = k\h + d,) + 2kh. 
So we can then always choose a positive k , such that 

{x\Q{G,a) + D\x) 

It then follows Q(G, a) + L> ^ 0. 

This expression is identical with that obtained in fTG]. For any graph 
G on n = pq vertices 

consider the degree condition A(G') = A(G^^). Now 

{L{G)f^ = (A(G) - A(Gr^)) + L(Gr^). 

Let 

D = A{G)-A{G^^). 
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Then D is an n x n real diagonal matrix with respect to the orthonormal 
basis 

\Vl) = \Ul) ® \wi), . . . , \Vpq) = \Up) (g) \Wq). 

Also 

Tr{D) = Tr{A{G)) - Tr{A{G^^)) = 0. 

We have two possible cases : D or D = 0. If D 0, that is the 
degree condition is not satisfied {i.e.A(G) ^ A(G^'^)), we have seen that 
L{G) + D ^ 0. As a consequence, L{G^^) + D ^ and then ^ 0. 

Hence (j{G) is entangled. ■ 
Lemma 2.4.17 : A graph (G, a) for a bipartite state corresponds to a 
separable state if {ij,kl} [i ^ k,j 7^ /) G E{G,a) =^ {il,kj} e E{G,a) 
and Qij^ki = CLii,kj- 

Proof : Suppose aij^ki = (iii,kj = CL^i k^j I. The contribution of the 
corresponding two edges is 

4P[-J=(iy) - \ki))] + pi-^(w - mm. 

which is a separable state. Thus all such pairs contribute separable states. 
Any other edge {ij, kl} with i = k or j = I has the contribution a^j /j;P[|z) 
— |/)))] which is separable. Loops contribute the product states 

..... ■ 

The reverse implication is not true in general. The counter-example is 

the graph (figure 2.12a) in example (2.4) which is separable. 
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A graph operation is a map that takes a graph to another graph |168] . We 



deal with four cases, namely deleting and adding an edge and deleting and 
adding a vertex. 

Deleting an edge {{vi^Vj}, ay^y.) from a graph [G^a] results in a graph 

def 

(G, a) - {{v^, v^), Gy^y^) = {V{G, a), E{G, a) \ {v^, Vj}) with ay^y^ = 0. Note 
the possibility Vi = Vj corresponding to the edge being a loop. Addition 

of an edge {{ui, vj}, aij) maps (G, a) to the graph (G, a) + ({-yj, -y^}, aij) 
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[V{G, a), E{G, a) U {-u^, Vj}] with a^.y. = aij. Deletion of a vertex Vi maps 

(G,a) to (G,a) - {v,} =^ [V{G,a) \ {v,},E{G,a) \ E,] where E, is the 
set of all edges incident to Vi (including the loop on Vi) with the weight 
function zero for the edges in Ei. Adding a vertex Vi to (G, a) maps (G, a) 

to (G, a) + {z;,} =^ {V{G, a) U {z;,}, E{G, a)). 

A very important point is that, in general, the set of graphs having a 
density matrix is not closed under these operations. Addition of an edge 
with positive weight and deletion of an edge with negative weight preserves 
the positivity of the generalized Laplacian resulting in the graph having 
a density matrix. However, addition (deletion) of an edge with negative 
(positive) weight may lead to a graph which does not have a density matrix. 
In the next section, we give a method for addition and deletion of an 
edge which preserves the positivity of the generalized Laplacian. Deletion 
and addition of vertices always preserves the positivity of the generalized 
Laplacian. 

Let B{7i^) be the space of all bounded linear operators on 7i^. A hnear 
map A : B{l-C^) B{l-C'^) is said to be hermiticity preserving if for every 
Hermitian operator O G B{H'^), A{0) is an Hermitian operator in B{l-C^). 
A hermiticity preserving map A : B{H^) B{H^) is said to be positive if 
for any positive operator O G B{H^), A(0) is a positive operator in B{7i^). 
A positive map A : B{W^) B{H^) is said to be completely positive if 
for each positive integer /c, (A ® /p) : B{H^ B{H^ H'^) is again 
a positive map. A completely positive map A : B{H^) B{H^) is said to 
be trace preserving if rr(A(0)) = Tr{0), for all O G B{7i'^). A quantum 
operation is a trace preserving completely positive map (for short, TPCP) 
11611 [7]. In standard quantum mechanics, any physical transformation of a 



quantum mechanical system is described by a quantum operation pj. We 
are going to use the following result: 

(Kraus representation Theorem) jl62f : Given a quantum operation 
A : B{H^) B{H^), there exist m x n matrices Ai, such that A{p) = 
'^AipA], where p is any density matrix acting on and ^AjAi = 

i i 

(the converse is true). The matrices Ai's are called Kraus operators. 
A projective measurement A4 = {Pi;i = 1,2,- •• ,n}, on a quantum 
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mechanical system S whose state is consists of pairwise orthogonal pro- 

n 

jectors Pi : Tis ^ T^s, such that = Idim{ns)- The i-ih outcome of the 

measurement occurs with probability Tr{Pip) and the post-measurement 
state of S is -^f^^- Whenever the i-th outcome of the measurement occurs, 
we say that Pi clicks. Last two paragraphs apply to complex Hilbert space 
and so also to real Hilbert space. 



2.5.1 Deletion and addition of an edge for a weighted graph 
with all weights > 

Here we describe how to delete or add an edge by means of TPCP. Our 
method of deleting an edge from a weighted graph with all positive weights 
is a simple generalization of the method in |77]. Let (G, a) be a graph 
on n vertices vi,--- ,Vn and m edges {vi.Vj^} ■ ■ ■ {v^^Vj^J,ik jk,k = 
I,-- - ,m and s loops {vi^Vi^} ■ ■ ■ {vi_^Vi^} . Our purpose is to delete the 
edge {v.^Vj^}, ik jk- Then we have 

^ m ^ s 



{G,a) 



and 



V2 



a{{G,a) - [v^^Vj^]) 



t=i 



'(G,a) 



- 2a 



m ^ 

i=l 



\v 



t=i 



A measurement in the basis M. = {-^(1^1^) ± |%))? \vi) ' i ^ 'ik,jk and 
z = 1, 2, • • • , n} is performed on the system prepared in the state ^(G, a). 
The probability that P+ = ^^[^(1^^^^) + |%))] chcks is 



Tr[P+a{G,a)] = E(^,|P+(7(G', a)| 
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(2.30) 

The state after the measurement is ^[-^(1^^^^) + |%)]- Let t/^ and t/^ 
be n X n unitary matrices such that [/+[--i=(|i;jJ + \vj^))] = -^i\vie) - \vje)) 
for £ = 1, • • • , /c - 1, /c + 1, • • • , m and t/^+J^d^) + 1%))] = ^ 
I ■ ■ ■ ■ ,s. Now, with probabiUty 2aj^j^/(t)(G',a) - '^CLim) we apply t/^J on 
+ for each £ = !,••• — + - ,m, and with prob- 

abihty ai^ij {Vi^G,a) - '^a^jk) we apply t/^ on ^[^(1^;^,) + |%>] for each 
t = 1, • • • , s. Finally we obtain cr((G', a) — {t'^t'^i;}) with probability given 
by Eq.(2.30). The probability that P[^{\v,^) - \vj^))\ clicks is 



2t) 



{ 7. ^j^iJ^jfei^ ~ ^^fei^ ~ ^ife*^ + ^ifejJ ~^z.^ititi^itik~^itjk) } (2-31) 

(G,a) frf 



£ = 1 *=1 



the state after measurement is P[-^i\vi,^) — \vjk))]- Let C/^^ and ?7^^ be n x n 
unitary matrices such that 

for £ = 1, • • • ,k — l,k + 1, ■ ■ ■ ,171 and 

Ukt-^iK) - 1%» = K) 

for t = 1, -,8. With probability 2a^fjJ(d^G,a) - ^^a^^jk) we apply t/^^ 

on P[^(\vi^} — 1%))] for each i = 1, • • • , A; — 1, A; + 1, • • • , m and with 
probability (G,a) - 2a^,^-J we apply [/^"^ on P[^(|i;,J - 1^;^,))] for 

each t = 1, 2, • • • , s. Finally we obtain (7((G, a) — {vi,^Vji,}) with probability 
given by Eq.(2.31). 
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The probability that -P[|fi)] where i ^ ik^jk and i = 1, • • • , n clicks is 



\G,a) 



t=i 



(2.32) 



and the state after measurement is Pflz;^)] . Let Ua and Uu henxn unitary 
matrices such that = -j^i\vu)~\vje)] for£ = 1, • • ■ k — l,k-\-l, • • • ,m 

and Uit[\vi)] = \vi^) for t = 1, • • • , s. With probability 2ai^jJ{d(G,a)-'^aikjk) 
we apply Uu on for each i = 1, • • • ,k — l,k -\- 1, • • • ,m and with 

probability ai^^J {d^G,a)-'^atkjk) we apply Uu on P[\vi)] for each t = 1, • • • , s. 

We obtain (7((G, a) — with probability given by Eq.(2.32). 

This completes the process. 

The set of Kraus operators that realizes the TPCP for deleting the edge 
{^i.^^jj is then 



U{, 



20,/ 



Ut 



A; — 1, A; + 1, • • • , m} 



2a 



ikjk 



t/+P[-^(K) + K))]:t=l,---,s} 



2n - ■ 1 

Ui<iP[^(\vi:) - 1%))] : ^ = 1, • ■ • , fc - 1, fc + 1, • • ■ m} 



'^{G,a) - 2a,: 



IkJk 



1 



'(G,a) 



2a. 



t/,-p[-^(K)-K))]:i=l,---,s} 



— "^^na^ t/i^P[|M,)] : i = 1, • • • ,n,ij^ik,jk-J = '^r-- 

k + 1, - ■■ ,m} 



a. 



-UitP[\vi)] : i = !,■■■ ,n,iy^ik,3k;t =!,■■■ ,s} 



f (G,a) - '^di^jk 

The set of Kraus operators that realizes TPCP for adding back edge 
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to iG,a) - {vi^Vj^} is. 



U 



U 



u 



u 



u 



i)(G,a) + 2a,,,, ^V2 



t'(G,a) + 2a,,,, ^V2 



"""" -V,-p[^{K)-\v,,))]:t = lr--.s 



i'(G,a) + 2a„„ ^^2 



2ai ■ 

-VuP[\vi)] ■.i = l,2,--- ,n,ij^ikJkJ = '^,'^,--- ,rn 



^{G,a) + 2aj,j, 



a 



-yifP[|t;i>] : i = 1, 2, • • • , n, i ^ ifc, jfc, ^ = 1, 2, 



/ f (G,a) + 2ai,„ 

where V^, V^^, V^, Vu, Vu are n x n unitary matrix defined as follows: 

1 1 

+ 1%)) = - l^j,)), for £=l,2,---,m 

yn^ihk) + 1%)) = K)^ for t = 1, • • • , s, 

- bjj) = - l^i.)). for ^ = l,2,---,m 

- 1%)) = K))^ for ^ = 1, • • • , s 

Vi^l^Ji) = -^(li;^,) - l^;,,)), for i = 1, 2, • • • , n, i 7^ ifc, j^, £ = 1, • • • , m 
= bzt), for i = 1, • • • 7^ = 1, • • • ,s. 



For deleting a loop {vi^,,Vi^,} a measurement in the basis {\vi),i = 
1, • • • , n} is performed on the system prepared in the state (7(G, a). Then 
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the probability that Pil'^i)] clicks for z = 1, • • • n is 



m 



\G,a) 



t = 1 



1=1 



(2.33) 



The state after the measurement is P[|'Ui)]. Let Uu be nxn unitary matrices 
such that ?7i4|vi>] = ^jd^i^) - \'^h))- For i = 1, • • • ,m and Uit[\vi)] = \vi^), 
for t = 1, • • • , - 1, + 1, • • • , s. With probability 2ai^jJ (G,a) - a^t/,**/) we 
apply ?7if on P[|vj)] for each £ = !,••• , m and with probability a^ij (f (G,a) — 
Oij/v) we apply Uu on -P[|?^i)] for each t = 1, ■ ■ ■ ,t' — l,t' + 1, ■ ■ ■ s. We 
obtain (7((G, a) — {vi^,,Vi^,}) with probability given by Eq.(2.33). 

The set of Kraus operators that realizes the TPCP for deleting the loop 



f (G,a) - av,v 



U,^P[\v^)\ z = l,--- ,m, £ = 1, 



m 



U 



f (G,a) - av,v 

The set of Kraus operators that realizes the TPCP for adding the loop 



2a 



H3(. 



f (G,a) + ^VV 



-yifP[|i;i)] : i = 1, • • • , n, t = 1, • • • , s L 



where V^^, V^^ are nxn unitary matrices define as follows : 
Vu\vi) = -^i\'^H) - 1%))^ for f = 1, • • • , m, z = 1, 
Viil^^i) = for t = 1, • • • , s,i = 1, • • • ,n. 



n 
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2.5.2 Deletion and addition of an edge with real weight, which 
preserves the posit ivity of the generalized Laplacian 

Let {G, a) be a graph with real weights on its edges not necessarily positive. 
We are basically concerned here with the deletion of {vi, vj} with a^.y. > 
and the addition of {vi,Vj} with a^^y. < 0, because in other cases the 
positivity of the Laplacian is preserved. We define the sets 

E+ = {{vi,Vj} G EiG,a),ay,y. > 0}, (2.34) 

E- = {{vi, Vj} e E{G, a),ay^y^ < 0} (2.35) 

and E = E+UE'. 

We define a graph operator S as 

E[E] = EU {{vi, Vi}, {vj, Vj} : a^.y. = ay.y. = 2\ay.y. \ and {vi, Vj} G E~} 

(2.36) 

Suppose we wish to delete a positive weighted edge {vi^^Vj,^} e E'"^, then 
we define the resulting graph as 

EC{{G, a) - {v,^,vjj) 

where the graph operator C is defined in (2.20b). 

For adding a negative weighted edge between Vi and Vj^i ^ j, we act 
on E(G, a) by the appropriate element of the set of operators {G^}, i,j = 
1 , • • • ,n,i j defined as 

Eij [E] = EU {{Vi,Vj}, {Vi,Vi}, {Vj.Vj} : ay.y. < O.ay.y. = 2\ay.y. \ = tty^y^} 

(2.37) 

To obtain the set of the corresponding TPCP operators we decompose 
the resulting graph, {G',a') given by EjC{{G,a) - {v,.^,Vj^}){ay^^y^^ > 0) 
(Eq. (2.36) ) or by G,, ((G,a) + {v,,Vj}){ay^y^ < 0) (Eq. (2.37)) or by 
((G,a) - {ui,,VjJ){ay^^y^^ < 0) or by (G, a) + {vi,Vj}){ay^y^ > 0) into 
spanning subgraphs determined by the sets E^ and E~ and treat the 
spanning subgraph corresponding to E~ replace the weights au,^y^ of edges 
{vi:,Vj} G E~ by —ay.y., so that both the spanning subgraphs have only 
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positive weights. For getting the Kraus operators we go through the fol- 
lowing steps. 

(a) First we determine the degree sums for the resulting graphs (C, a') 
in four cases. 

(i) Deletion of a positive weighted edge {w^^, Vj^} 

i {ui,Vj}eE- 

(ii) Addition of a positive weighted edge {vi, Vj}. 

^{G',a') = ^{G,a) + '^CLy^vj- (2.39) 

(iii) Deletion of a negative weighted edge {vi^, Vj,J\ 

^{G',a')='0{G,a)-'^a^i,v,^. (2.40) 

(iv) Addition of a negative weighted edge {vi, vj} 

^{G',a') = f(G,a) + 2a^;,,^. + 4|a^.,^.|. (2.41) 

(b) We construct the Kraus operators separately for and G~ for 
deleting the same edge {vi, vj} from G^\J{vi, vj} or adding the edge {vi, vj} 
to G^, using the method given in section 2.5.1. However, the probabilities 
of applying various unitary operator U^^ and [/^, Ua and Ua are determined 
using f (G'^a') as in step (a) above. 

(c) Let {Ai} and {Bi} denote the sets of Kraus operators for the graph 
operations on G'^ and G~ as described in (b). Then 

(t{G', a') = A,a{G, a)4 - ^ Bja{G, a)B] (2.42) 

i j 

and 

J2a]A-J2b]B, = I (2.43) 

i j 

which can be justified by construction. 
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We comment here that it is possible to modify the graph, after deleting 
a positive edge or adding a negative edge, which can preserve positivity 
in different ways, leading to different sets of Kraus operators. The basic 
idea is to add new loops. In our method we try to minimize the addition 
of loops. Further, in our method we cannot reverse the graph operation 
for deleting a positive edge or adding a negative edge. But this is not a 
problem since the quantum operations given by super operators are, in 
general, irreversible. 

2.5.3 Deleting Vertices 

In order to delete a vertex Vi from a graph (G,a), 

(i) Delete edges, including loops, on Vi, one by one, by successively apply- 
ing the procedure in section 2.5.2. The resulting graph (C, a') has a 
density matrix with the i-th row and i-th column containing all zeros. 

(ii) We now perform, on a{G', a'), the projective measurement M = {/„ — 
P[|t'j)], P[|i;i)]}. Since P^'^i)] is the matrix with all elements zero 
except the i-th diagonal element, while a{G', a') as all zeros in the i-th 
row and column, the probability that Pfl'^i)] clicks = Tr{aP[\vi)]) = 
0. Thus when M is performed on a{G' , a')' , In — P[\vi)] clicks with 
probability 1 and the state after measurement is (t{G', a') — {vi}) and 
is the same as cr(G', a') without i-th row and i-th column. 

2.5.4 Adding a Vertex 

Let (G, a) be a graph on n vertices vi, • • • and m edges {vi^^, t'j^}, k = 

,m,ik ^ jk and s loops {vi^.Vi^}, t = I,-- - ,s; 1 < ik,Jk,H < n. 
Consider the following density operator 



where |w2)} form an orthonormal basis of C^. We associate vertices 

Ui,i = 1, 2 to the state \ui). Consider the graph 

H = ({wi, U2}, {{ui, Ui}, {W2, U2}}) 
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p = a{{H,b)B{G,a)) 



a{H, b) a{G,a). 



Thus {H, 6)Q(G, a) is the graph on 2n vertices labeled by uiVi, • • • , ui'Un, W2'i^i, 
• • • ,U2Vn and with 2m edges and 2s loops (see section 2.4.2) 

{uiVi^,UiVj^} • • • {uiVi^,UiVjJ{u2Vi^,U2VjJ • • • {u2Vi^, U2VjJ 

and loops {uiVi^^uiVi^}^{u2Vi^,U2Vi^],t = I,-- - ,s. So {H^b) □ (G, a) = 
(ifi, ai) l±) (if25 ^2) where 

(//l, ai) = • • • UiVn}: {{uiVi^.UiVj^} ■ ■ ■ {uiv^^, UlVj}}) 

{H2, 02) = {{U2VI ■ ■ ■ U2Vn}, {{u2Vi^, U2Vj^} ■ ■ ■ {w2^z™, U2Vj^}}). 

We first delete all edges and loops of {H, b) □ (G, a) which are incident 
to the vertex U2V1 G V{H2,a2) as in section 2.5.2. Now we perform the 
following projective measurement on (j{{H, b) □ (G, a)), 



The probability that l2n — XI P[^2Vi)] is 1 and the state after the measure- 

ment is a{{Hi^ ai) + {u2Vi}). 

Example (2.8) : Consider the graph as given in the figure 2.16 , we want 
to delete the edge {1,2} with positive weight by means of TPCP. Calculate 
the Kraus operators for G^ and G as in section 2.5.2 , where Ai for 

i = 1, 24 and Bi for G~ ,i = 1, ,4 and substitute in the following 

equation 



n n 



i=2 1=2 



n 



a 



{G\ a') = ^ AiuiG, a)A\ - ^ BjuiG, a)B] 
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where 



11 , 12 




21 ' 22 



Figure 2.16 



a{G,a) = i 



1-1-11 

-1 3 -1 -1 

-1 -1 3 -1 

1 -1-1 1 



we get 



^(G", o!) = - 



2 0-11 

2-1-1 

-1 -1 3 -1 
1-1-13 



and 



24 



Y,A\A,-Y,B]B, = I 



1=1 



2.6 Representation of a general Hermit ian operator 
by a graph 

In this section, we generahze sections 2.2 - 2.4, to quantum states in a 
complex Hilbert space, that is, to the density matrices with complex off- 
diagonal elements. We have also given rules to associate a graph to a 
general Hermitian operator. We believe that any further advance in the 
theory reported in this chapter will prominently involve graph operators 
and graphs associated with operators. 
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2.6.1 Representation of a general density matrix with complex 
off diagonal elements 

Consider an n x n density matrix with complex off-diagonal elements. We 
associate with this density matrix an oriented graph (G, a) on n vertices, 
m edges and s loops with weight function 

a : V{G) X ^ C. 

The weight function a has the following properties: 

(i) a{{u,v}) 7^ if {u,v} G E{G,a) and otherwise. 

(ii) a{{u, v}) = a^d-u, u}) 

we write a{{u, v}) = |a({M, e*^"", 0^;^; = 0. 

Note that when = lir,! = 0, 1, • • • , i.e. a{{u, v}) is real, positive 
when I is even and real negative when I is odd. 
The degree Xiy of vertex v is given by 

^{G,a)iv) = ^v= + «({^'^}) (2-44) 

uGV{G,a),uj^v 

f(G,a) = Y 

vGV{G,a) 

The adjacency matrix M(G, a) of a complex weighted graph with n 
vertices is an n x n matrix whose rows and columns are indexed by vertices 
in V{G,a): 

Muv = a{{u,v}) = a%{v,u}) = {M^uT- 

The degree matrix A(G, a) of the complex weighted graph is an n x n 
real diagonal matrix, whose rows and columns are labeled by vertices in 
a) and whose diagonal elements are the degrees of the corresponding 
vertices. 

A(G, a) = diag\Pv] v G V{G, a)] 
where dy is given by equation (2.44). 
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The loop matrix Ao((j', a) of a graph (G, a) is an n x n real diagonal 
matrix with diagonal elements equal to the weights of the loops on the 
corresponding vertices 

The generalized Laplacian of a graph (G, a), which includes loops, is 

Q{G,a) = A{G,a) + M{G,a) - Ao(G,a) (2.45) 

Note that Q{G, a) is a Hermitian matrix . If the generalized Laplacian 
Q{G, a) is positive semidefinite, we can define the density matrix of the 
corresponding graph (G, a) as 

a) = Q{G, a) (2.46) 

where Tr{a{G,a)) = 1. 

For any n x n density matrix a with complex off diagonal elements we 
can obtain the corresponding graph as follows: 
Algorithm 2.6.1 : 

(i) Label the n vertices of the graph by the kets from the standard or- 
thonormal basis. 

(ii) For every nonzero ijth element with j > i given by a{{i,j}) draw an 
edge between vertices labeled \vi) and \vj), with weight a{{i,j}). 

(iii) Ensure that '0v^ = era by adding loop of appropriate weight to Vi if 
necessary. 

Example (2.9) : (1) 



1 


" 1 -z " 


1 


1 e-'''/'^ 


2 


i 1 


~ 2 


g-i7r/2 I 



where 1?/,+) = "^(|1) + ^|2)) and the corresponding graph is as shown in 
figure 2.17a 

(2) 



1 


" 1 i' 


1 


1 e''^/2 


2 


-i 1 


~ 2 
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i -i 1 



(a) 



I i 1 




(b) 

Figure 2.17 

where \y, —) = -^(ll) — ^|2)) and the corresponding graph is as shown in 
figure 2.17b 
(3) 



1 

% 
% 

-1 



1 



The corresponding graph is as shown in figure 2.17c. 

Note that remark 2.2.1 is valid also for complex weighted graphs. 
Remark 2.6.2 : Theorem 2.2.3 applies to complex weighted graphs with 
equation (2.11) changed to 



(2.47) 



i=\ k=l 

also, lemma 2.2.4 applies to complex weighted graphs. 
Definition 2.6.3 : A graph (H, b) is said to be a factor of graph {G, a) if 
V{H, b) = V{G, a) and there exists a graph (H', b') such that V{H', b') = 
V{G,a) and M(G,a) = M{H,b) + M{H',b'). Thus a factor is only a 
spanning subgraph. Note that 



by.y. if {vi.Vj} e E{H,b) 
Kv, if {vi,Vj}eEiH',b') 



Now let (G, a) be a graph on n vertices f i, • • • having m edges 
{^^l,^Jl}, • • • and s loops {vi^.Vi^} ■ ■ ■ {vi^.Vi^} where 1 < nji, 



^l^2 



is < n. 
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Let {Hii^j^, ciikjk) be the factor of (G, a) such that 

\M(H n \] - / ^"'^'^ if w = ifc and w = jk or a*^^.^if u = jk,w = ik 
[M[iii^j,,ai,jj\u,nj - Q otherwise 

(2.48) 

[A{H,,,,^a,,j,)U = l M .^ = ^^ = ^ ^ = i. = «^ (2.49) 
L V «fcjfc' «Ufeyjw,«^ I otherwise ^ ^ 

Let {Hi^^i^^a^i^) be a factor of {G^a) such that 

a,..)l„.„ = [A(/f..... «..,.)].,„ = I "^'iZ^Zs^ = = " (2.50) 

Theorem 2.6.4 : The density matrix of a graph {G,a) as defined above 
with factors given by equations (2.48), (2.49) and (2.50) can be decomposed 
as 

^ m 1 * 

a{G,a) = - V2|a({ifc,ifc})|(7(/fi,j,,ai,jJ + y^^ith^i^hit^^Hit) 

^{G,a) ^(G,a) ^ 



or 

TO 



1 1 

a{G,a) = - ^2|a({z,,j,})|P[-=(|^,,) - e 

1 ^ 



%))] 



Where ^^^^j^ = tt for any edge {ik, jk} with real positive weight and = 
for any real negative weight. 

Proof : From equations (2.48), (2.49), (2.50) and remark 2.6.2, the density 
matrix 



'^\^ikjk\ 



is a pure state. Also, 
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is a pure state. Now 

m s 

5 ^itU) 



k=l 



t=l 



t=l 



k=l 

Therefore, from Eq. (2.46) 

1 



(7(G, a) = 



m 



m 



.A;=l 



1 



'(G,a) 



^ m 1 

^(G,a) ^ 0(G,a) ^ 

In terms of the standard basis, the ww;-th element of matrices a{Hi^,j^,, o^ikjk) 
and a{Hi^,^,a,^,J are given by {vu\(7{H^^j^, , ai,^jj\vuj) and (i;y|cr(//ij^^ai^iji;^i;), 
respectively. In this basis 

Therefore equation (2.51) becomes 

777. 

a{G,a) = J_^2|a({4,i,})|P[i.(|i;,J -e^^^^^H%» 

(2.52) 
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Where 04^-^, = tt for any edge {zfc, jk} with real positive weight and ^^^j^ = 
for any real negative weight. ■ 
Example (2.10) : 

(i) For a graph given in figure 2.17b, the density matrix is 

<7{G,a) = i{2P[-i=(|l>-e"/^|2»]} = P[-J=(ll)-»|2))l 

(ii) For a graph given in figure 2.17c, the density matrix is 
a(G,a) = i{2P[-i=(|ll)-e-^-/^|12))] + 2P[-i=(|ll)-e-^^^^ 

+2P[-i.(|ll>-|22))] + 2P[-i.(|12) + |21))] 

+2P[-i.(|12) - 6-^^^/2122))] + 2P[-J=(|21) - e-/2|22))] 
-2P[|11)] - 2P[|22)] - 2P[|12)] - 2P[|21)]} 



a{G, a) = ^ 



1 -I -I -1 

I I I -I 

I I I -i 

-I I I I 



Example (2.11) : Consider the state 

cy = \pu+)\y.+)] + \pu+)m 

where |y, +) = ^(|1) + i\2)) and = ^(|1) + iV2\2)) 

7 -(3 + 4v/2)i -7i 

J_ (3 + 4v/2)i 11 3 + 4v/2 

36 7i 3 + 4v/2 7 

_-(3 + 4\/2) Hi (3 + 4v/2)i 



(7 



-(3 + 4^2) 

-III 
-(3 + 4V2)z 
11 



1 

36 



7 (3 + 4x/2)e-^'^/2 7e-^'^/2 -(3 + 4^2) 

(3 + 4^2)6^^/2 11 7 lle-*''/^ 

7eW2 3 + 4^2 7 (3 + 4^2)6-^^/2 

-(3 + 4^2) ll6^''/2 (3 + 4^2)6^'^/^ 11 
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-(6+8K2) ,^ 




Figure 2.18 



(6+8VT) 



-(6+81/2) 



21 -(3+4l/^)i 22 



The corresponding graph is as shown in figure 2.18, and using the equation 
(2.45) to get the matrix from graph in figure 2.18, 

<T(G,a) = l{2(3 + 4V^)P[-J=(lll>-e-"/^|12»] 

+2 X 7Pl^{\n) - e-"/^|21»] + (3 + 4v/2)P[i=(|ll> - |22»] 

+2(3 + 4V2)P[^(|12) + |21))] + 2 X 11P[^(|12) - e-"'''|22))] 
V 2 V 2 

+2(3 + 4V2)P[^(|21) - e-'^/^\22))] - (6 + 8V2)P[|11)] 
v2 

-(6 + 8v^)P[|22>] - (6 + 8v^)P[|12)] - (6 + 8v^)P[|21)]}. 

We can check that 

7 -(3 + 4V2)i -7i -(3 + 4\/2) 

(3 + 4\/2)i 11 3 + 4^/2 -Hi 

7i 3 + 4v^ 7 -(3 + 4\/2)i 

_ -(3 + 4^/2) Hi (3 + 4V2)i 11 

We can also check that this state is not pure by applying remark 2.6.2 on 
the graph. 

Remark 2.6.5: Lemma 2.2.10 and remark 2.2.11 are valid for the complex 
weighted graphs with disjoint edge union + 



(j(G, a) = — 
^ ' ^ 36 



2.6.2 Separability 

Remark 2.6.6: The definition of the tensor product {G,a) <S) {H,b) of 
two complex weighted graphs (G, a) and {H, b) is the same as given before. 
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However, note that {vi,V2} G E(G,a), {wi,W2} £ E(H,b) implies 
c{{{vi,wi),{v2,W2)}) = a{{vi,V2})b{{wi,W2}) 

and 

c{{vi,W2),iv2,wi)}) = a{{vi,V2})b{{w2,wi}) = a{{vi,V2})b%{wi,W2}). 

Remark 2.6.7 : Equations (2.18a) and (2.18c) are valid for the tensor 
product of complex weighted graphs. Also, Q{{G, a) <^ {H, b)) ^ Q{G^ a) ® 
Q{H,b). Equation (2.18b) holds good only for graphs without loops, for 
graphs with only loops or when one factor has no loops and other factor 
has only loops. For such graphs equation (2.18b) immediately gives 

2.6.3 Modified tensor product 

The modified tensor product of two complex weighted graphs requires the 
operator J\f to be redefined in the following way. We replace the equation 
(2.20a) by 

^ \a{{vi,Vk})\ + a{{vi,Vi}) (2.53) 

VkGV{G,a) 

The definitions of the operators rj, C and Q, remain the same. Equations 
(2.21) to (2.24) are satisfied by these operators on the complex weighted 
graphs. We further have 

(i) M{AfC{G, a)) = A{G, a) - Ao(G, a) 
A(AA/:(G, a)) = A(G, a) - Ao(G, a) 
AoiAfC{G, a)) = A(G, a) - Ao(G, a) ^ ^ ^ 

Q{NC{G, a)) = A(G, a) - Ao(G, a) 

The modified tensor product of two complex weighted graphs (G, a) and 
(iif, b) with p and q {> p) vertices respectively is 

(G, c) = (G, a) □ {H, b) = 
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a) (g) CiH, h) + a) (g) M{H, h) + a) (g) 6)+ 

a) (g) 6) U 2A/'/:(G, a) (g) J\fCr]{H, b)} (2.55) 

The weight function c of (G, a) □ 6) is obtained via the definition of 
tensor product and the disjoint edge union. 
Lemma 2.6.8 : A((G, a) □ {H, b)) = A(G, a) A{H, b). 
Proof : Since lemma 2.2.10 applies to disjoint edge union of complex 
weighted graphs, 

A((G, a) □ (if, 6)) = A{C{G, a) (g) C{H, b)) + A{C{G, a) (g) 7V(if, 6))+ 

A(Ar(G, a)®C{H, 6)) + A(^^(G, a)^n{H, b)) + A{2NC[G, a) ®J\fCr]{H, b)) 

The last two terms are justified because the graphs involved are real weighted 
graphs. Using remark 2.6.7 we get 

A((G, a) □ (i/, b)) = A(£(G, a)) (g) A(i:(if, 6)) + A(£(G, a)) (g) A(A/'(i/, 6))+ 

A(A/'(G, a)) A(/:(//, 6)) + A(Q(G', a)) A(Q(//, 6))+ 

2A(A/'£(G, a)) (g) A{AfCr]{H, b)) 

Using equations (2.21) to (2.24) and (2.54), we get, after some simplifi- 
cation, 

A((G, a) B {H, b)) = A((G, a)) A{{H, b)) 

Corollary 2.6.9 : d^G,a)B{H,b)iv,w) = V{G,a){v) ■ '0{H,b){w) 
and 

^{G,a)B{H,b) = f (G,a) " ^{H,b) 

Proof : The first result follows directly from lemma 2.6.8. For the second 
note that 

rr(A((G, a)B{H, b))) = Tr(A(G, a)®A{H, b)) = Tr(A(G, a))-Tr{A{H, b)) 
where Tr denotes the trace. ■ 
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Theorem 2.6.10 : Consider a bipartite syatem in C^^O in the state a. 
Then o" = o"i (T2 if and only if a is the density matrix of the graph (G, a) El 
b) where (C, a) and (if, b) are the graphs having density matrices ai 
and (72, respectively. 

Proof : If part : Given (G, a), (ii, b) we want to prove 
(7((G, a) □ (if, 6)) = (7i(G, a) (g) (72(ii, b). 
From the definition of the modified tensor product we can write 

a[{G, a) □ (if, b)) = {Q[C{G, a) C{H, b)+ 

0{G,a)B{H,b) 

C{G, a) Af{H, b) + A/'(G, a) O C{H, b) + {n{G, a) ^(if, 6))U 

Using remark 2.6.5 and corollary 2.6.9 we get 

a((G, a)H(ir, b)) = ^^^^^[QiCiG, a)®C{H, b))+Q{C{G, a)m{H, 6))+ 

g(A/'(G, a) (g) £(//, 6)) + g(f^(G', a) ® n{H, b) U 2MC{G, a) (g) MCr]{H, b))] 

(2.56) 

We can calculate every term in Eq.(2.56) using Eqs.(2.21) to (2.24) and 
Eq.(2.54) and substitute in Eq.(2.56) to get 

(j((G, a) B {H, b)) = a{G, a) a{H, b). 

Only if part : Given a = cti (72 consider the graphs (G, a) and {H, b) for 
(7i and (72 respectively. Then the graph of a has the generalized Laplacian 

Q{G,a)^QiH,b) = (A(G,a) + M(G,a)-Ao(G,a))0(A(/f,6) + 

M{H,b)-Ao{H,b)) 
= A{G,a)0A{H,b) + A{G,a)0{M{H,b)- 
Ao{H, b)) + (M(G, a) - Ao(G, a)) A(ii, 6) + 
(M(G, a) - Ao(G, a)) (M(^, 6) - Ao{H, b)) (2.57) 

Using equations (2.21) to (2.24) and (2.54) we see that RHS of equation 
(2.57) is the generalized Laplacian for (G, a) □ {H, b) ■ 
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Remark 2.6.11 : The proof that the modified tensor product is associative 
and distributive with respect to the disjoint edge union is the same as that 
for the case of real weighted graphs (corohary 2.4.7). 

Remark 2.6.12 : The definition of the Cartesian product of graphs is the 
same as given in definition 2.4.8. 

Remark 2.6.13 : Corollaries 2.4.9 and 2.4.10 apply to complex weighted 
graphs without any change. 

2.6.4 Convex combination of density matrices 

Consider two graphs (Gi,ai) and {02,0,2} each on the same n vertices, 
having a{Gi, ai) and (t{G2, 02) as their density matrices respectively, where 
ai and a2 are complex weight functions. Let (G, a) be the graph of the 
density matrix a{G,a) which is a convex combination of cr(Gi,ai) and 

cr(G2, 02), 



It is straightforward, using the definitions of the operators J\f, C and 77, 
to verify that 



(G, a) = [W(Gi, ai) U (1 - A)A/'(G2, 02)] U [A/:(Gi, ai) U (1 - A)£(G2, 02)] 



We can apply this equation to any convex combination of density matrices. 
Let 



a{G, a) = Xa{Gi, ai) + (1 - X)a{G2, 02), < A < 1. 



Ur]C[XjC{Guai) U (1 - X)C{Gi,a2)]. 



(2.58) 





where a and {ai} are complex weight functions, a{{vi, v^}) = ^'^{{vi, Vk}) 
and a{{vi, vi}) = J2i ^'iii'^h M) with a\ = pitti. 
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Lemma 2.6.14 : Let (Gi,ai), (^2,02) and (G, a) satisfy Eq.(2.58). Then 

a) = ai) + ^^^(7(^2, 02). 

^{G,a) f(G,a) 

Proof : Similar to that of lemma 2.2.10. ■ 
In general, if {G,a) satisfies Eq.(2.59) for some set of graphs {{Gi,ai)}, 
we have 

<^{G,a) = -^J2^^G.,aMG^.a^). (2.60) 

Theorem 2.6.15 : Every graph {G,a) having a density matrix a{G,a) 
can be decomposed as in Eq.(2.59), where a{Gi, ai) is a pure state. 
Proof : The same as that of theorem 2.2.12. ■ 
Corollary 2.6.16 : A state of a /c-partite system is separable if and only 
if the graph (G, a) for a has the form 

(G, a) = [U,A/- Hj^i (GJ, ai)] U B^, (Gj, aj)] 

UrjClU.CB'^.iGiai)]- (2-61) 



Proof : The same as of corollary 2.4.11, where we refer to theorem 2.6.4 
instead of theorem 2.4.5 and lemma 2.6.14 instead of lemma 2.2.10 and 
Eq.(2.60) instead of remark 2.2.11. ■ 
Corollary 2.6.16 says that Werner's definition [25j of a separable state 
in C^^ (g) • • • C^'= system can be expressed using corresponding 
graphs. 



2.6.5 Representation of a Hermitian operator (observable) by a 
graph 

In order to represent a general Hermitian matrix corresponding to a quan- 
tum observable A we lift the requirement that the Laplacian be positive 
semidefinite and rr[74] = 1. In other words we take the generalized Lapla- 
cian as the matrix for the graph. 
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Given a Hermitian matrix A, the algorithm 2.6.1 can be implemented 
to get its graph (G, a). The corresponding observable A of a graph (G, a) 
is 

A = Y^ 2a,,j,P[-^{K) - e^^'.^-. \vj,))] + ^hhP[K)] (2-62) 



k=l 



t=l 



Example (2.12) 


: Give the | 


^raph of 


(1) = 


"or 

1 






The corresponding graph of ax is 


(2) ay = 


-i 

1 




e 



The corresponding graph of ay is shown in figure 2.19b 



(a) 

Figure 2.19 



(b) 



Using Eq. (2.62) to get the operators from graphs, 



a, = -2P[^(|1> - |2»] + F[|l>] + F(|2>] = |1>(2| + |2>(1| = 

ay = 2P[-L(|l)-e-"/^|2»]-P[|l)]-Pl|2)] = -»|l>{2|+j|2>(l| = 

2.7 Some graphical criteria for the positive semidef- 
initeness of the associated Laplacian 

In this section we address the following question. Given a graph, can the 
positive semidefiniteness of the associated Laplacian be determined using 
the topological properties of the graph? A general answer to this question 
seems to be difficult because the theory of weighted graphs, with negative 
and complex weights is almost unavailable. Many results obtained for 
simple graphs do not apply to the weighted graphs with real or complex 



1 

1 
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weights. Nevertheless, we give here the above mentioned criteria in some 
special cases. 

Lemma 2.7.1 : Let (G, a) be a graph with real or complex weights, having 
one or more nonisolated vertices with degree zero. Then the Laplacian of 
(G, a) is not positive semidefinite. 

Proof : Such a graph (G, a) has one or more zeros along the diagonal of 
its Laplacian with nonzero entries in the corresponding rows. However, a 
Hermitian matrix with one or more zeros in its diagonal has at least one 
negative eigenvalue unless all the elements in the corresponding rows and 



columns are zero [ jl69j . ■ 
Lemma 2.7.2 : Let {G,a) be a n vertex graph with real weights, having 
at least one loop and let the weights on all the loops be negative. Then 
the Laplacian of (G, a) is not positive semidefinite. 
Proof : For the given (G, a) and some x in we have 



x^[Q(G,a)]x = ^az,j,(x,, - Xj,f - ^ 



,2 



where the first sum is over edges and the second sum is over loops. It is 
easy to check that x^[Q(G, a)]x < for a: = (1 1 1 ... 1)^. ■ 
Lemma 2.7.3 : Let (G, a) be a graph without loops satisfying a{u,v) = 
a^^e**^"^, {(j)uv 7^ 27rn). Then the associated Laplacian is positive semidefi- 
nite. 

Proof : This follows directly from theorem 2.6.4. ■ 
Observation 2.7.4 : Let (G,a) be a graph satisfying a{u,v) = a^^e*"^™ 
and a{{v, v}) > for all vertices in V{G, a). Then the associated Laplacian 
is positive semidefinite. 

Proof : The Laplacian is a Hermitian matrix which is diagonally domi- 
nant. Therefore, by Gersgoin circle criterion jl65l 11661 1170] it is positive 
semidefinite. ■ 
On a n vertex graph (G, a), we define a new graph operator Q{ui) which 
deletes the vertex Ui and rolls the edges incident on Ui into loops with same 
weights on the edges connecting neighbors of Ui as shown in figure 2.20. 
We call the resulting subgraph principal subgraph. The Laplacian of the 
principal subgraph obtained by operating Q{ui) on [G^a] is the principal 
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submatrix of the Laplacian of (G, a) obtained by deleting ith row and ith 
column. 




6(4) 



Cf— 



4 3 

Figure 2.20 



Lemma 2.7.5 : If one or more principal subgraphs of {G,a) are not 
positive semidefinite, then (G, a) is not positive semidefinite. 
Proof : This follows from the result that all the principal submatrices of 
a positive semidefinite matrix are positive semidefinite jl66] . ■ 



Lemma 2.7.6 : Let {G, a) be either a n vertex tree (n > 2) or a n vertex 
cycle (n > 4). We assume that there are no loops in (G, a) and that a(w, v) 
is real for all {u,v} G E{G,a). Then (G, a) has a positive semidefinite 
Laplacian if and only if a({w, v}) > for all (w, v) G E{G, a). 
Proof : Only if part : We prove that a{{u,v}) < for any one {u^v} G 
E'(G, a) =4> Q(G, a) ^ 0. Let (T, a) be a tree with -ui, . . . , vertices, and 
let {vi^Vi+i} be an edge in (T, a) with negative weight a{{vi,Vi+i}) < 0. 
We operate on (T, a) by Q{vi+i). There are two possibilities. If vi+i is 
a leaf, we get only one component with a negative weighted loop on Vi. 
By lemma 2.7.2, the Laplacian of this principal subgraph is not positive 
semidefinite and by lemma 2.7.5 the Laplacian of (T, a) is also not positive 
semidefinite. If "U^+i is not a leaf then Q{vi+i) will result in two or more 
principal subgraphs. The principal subgraph containing vertex Vi is a graph 
having one loop with negative weight. By lemma 2.7.2 the Laplacian of 
this principal subgraph is not positive semidefinite and from lemma 2.7.5 
Q{T, a) ^ 0. Let G„ be an n-cycle and let a{{vi, Vi+i}) = a < 0. We operate 
by Q{vi+i) which results in a n — 1 vertex path, say Pn-i with Vi having 
a negative loop and having a positive or negative loop. If both the 
loops are negative we can use lemma 2.7.2 and lemma 2.7.5 in succession 
to show that Q(G„,a) ^ 0. Suppose the loop on Vi+2 is positive. Then for 
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some X £ we have 

n-2 

x^Q{Pn-i, a)x = a,^j^{x,^ - xj^f + a{v^+l, v,+2)xl.^^ - \a\xl.. 

k=i 

It is straightforward to check that x^Q{Pn-i, a)x < for = (1 ... 1 . . . 1), 
that is a vector x with all components 1 except Vi+2th component which is 
zero. Thus Q{Pn-i,a') ^ 0. By lemma 2.7.5 Q{G,a) ^ 0. 
If part : Assume Q{G,a) ^ 0. This implies that there exists at least one 
X £ satisfying 

x^Q{G,a)x = Ya{ikJk)ixi^ -Xj^f < 0. 

k 

Since {xi^ ~ Xj^Y > for all k, the above inequality is satisfied only 
when a{ik,jk) < for some k. This proves the if part. ■ 

We observe that the proof of if part applies to all graphs as it should. 
Lemma 2.7.7 : Let all loops on a graph {G, a) have real positive weights. 
Let every edge {u, v} G E{G, a) having a{u, v) < he common to pair of 
C3. Let all such pairs of C3, each containing a negative edge be disjoint. 
Let all the edges in each pair of C3, other than the contained negative edge 
have positive weights satisfying a{u, v) greater than the absolute value of 
the weight on the negative edge. Then the Laplacian of (G, a) is positive 
semidefinite. 

Proof : Consider a negative edge common to two G3S as shown in the 
figure 2.21. 




a> 

Figure 2.21 



By hypothesis bi > a, i = 1,2,3,4. We can write 6^ = a + q,q > 
0, i = 1, 2, 3, 4. We can decompose this graph as the edge union as shown 
in figure 2.22. 
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The first graph on RHS has aU positive weights and hence has a positive 
semidefinite Laplacian. It is straighforword to check that the second and 
the third graphs on RHS correspond to the projectors P[-^{\j) ~ + 

and P[^{\j) ~ + \k))] respectively. Hence they have positive 
semidefinite Laplacians. The Lapalcian of the graph on LHS is the sum 
of the Laplacian of the graphs on RHS (lemma 2.2.10), each of which is 
positive semidefinite. But we Know that the sum of positive semidefinite 
matrices is a positive semidefinite matrix jl66j . Now the graph (G, a) 
can be written as edge union of the factors (spanning subgraph) as figure 
2.21 (possibly more than once) and the remaining factor which has all 
positive weights. The Laplacian of each factor is positive semidefinite and 
the Laplacian of the given graph, being the sum of positive semidefinite 
matrices, is positive semidefinite. ■ 
Lemma 2.7.8 : If all the negative edges of a real weighted graph {G,a) 
occur as in the following subgraph as shown in figure 2.23, where q > 
6; z = l,2,...,8 and 6 > a > 0; then the associated Laplacian is positive 
semidefinite. 

Proof : We can decompose the above graph into factors as shown in figure 
2.24 

From the graphical equation in figure 2.24 we see that the first factor on 
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b > a > 

Figure 2.23 



c2 k ^3 j i b _ 



u 




Figure 2.24 



RHS corresponds to )|— )], the second factor has a positive semidefinite 
Laplacian from lemma 2.7.7 and the third factor has a positive semidefinite 
Laplacian as it has all positive weights. Since this graph occurs (once or 
more) as disjoint subgraphs of (G, a) it can be written as edge union of 
one or more of these subgraphs and the remaining graph containing only 
positive or complex edges. Since each of these has a positive semidefinite 
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Laplacian, (G, a) also has a positive semidefinite Laplacian. ■ 
Lemma 2.7.9 : Let {G'^'\a) be a complete signed graph with weight 
function aij G { — 1,1} without loops on 2" vertices n > 1. Let Ei denote 
the set of edges incident on ith vertex {\Ei\ = 2" — 1) and let Ef^E^ 
denote the sets of edges incident on the iih vertex with weight +1 and —1 
respectively, [Ei = E^ + E^). Let (G^",a) satisfy the following condition 
(i), \E'[\ = 2"'~-^ — 1, i = 1, 2, . . . , 2", so that the degree of every vertex= 1. 
Then (G^", a) corresponds to a pure state in 2"^ dimensional Hilbert space. 
Proof : We need to prove that condition (i) in the statement of the lemma 
can be realized for all n and that the resulting signed graph corresponds 
to a pure state for all n. We use induction on n. It is clear that the 
assertion is true for n = 1 with the corresponding pure state given by 
P[-i=(|l) — 12))]. Now assume that assertion (that is condition (i) and purity 
of the corresponding state) is true for n = k. For the graph corresponding 
to n = /c + 1 with \V{G, a) \ = 2^+^, consider the modified tensor product 

(G^a)□(G'^a) = iG^'^\a) 

= {£(G'^ a) (g) Cr]{G''\ a)} + {C{G^ a) N'{G''\ a)} + 

{X{G^ a) (g) C{G^\ a)} + {Q{G^ a) Q{G^\ a)}(2.63) 

where C^rj^N and are graph operators defined in equation (2.20b) and 
C^, G^ are graphs with number of vertices 2 and 2^ respectively. Note that 
the last term corresponds to an empty graph as G does not have any loops. 
Since the modified tensor product of two complete graphs is also a complete 
graph, (G^'^^^a) is a complete graph. Therefore |£'i(G^''^\ a)| = 2^+^ — 
1,1 = 1,2,...,2^+^ To show that condition (i) is realized for [G'^^^^.a) 
given the induction hypothesis, we note that the first term in equation 
(2.63) contributes \E^{G'^^ ,a)\ negative edges to (l,z)th vertex in (G^*"^^) 
and the third term contributes \E~iG'^ , a) \ negative edges, while the other 
two terms have no contribution. Therefore 

\E^,{G''^\a)\ = \EriG'\a)\ + \E7{G'\a)\ = 2^-1 

Similarly, the first three terms contribute |£^^~(G^'', a)|, 1 and \E^(G'^'' , a)\ 
positive edges to (1, z)th vertex. Therefore, 

\E+{G^'^\a)\ = \E^{G^\ a) I + \E+{G^\ a)\ + l = 2^-^ - 1 + 2^-^ + 1 = 2^ 
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and similarly for £"2^. That {G^^^^ , a) corresponds to pure state follows from 
the fact that the state corresponding to the modified tensor product of two 
graphs is the tensor product of the states corresponding to the factors. 
Since the state for (G^ , a) is pure by induction hypothesis and {G'^^a) is 
pure, the preceding statement means that(G2'^\a) is a pure product state. 



2.8 Summary and Comments 

Following is a brief summary of the main features of the chapter 

(i) We have given rules to associate a graph with a quantum state 
and a quantum state to a graph, with a positive semidefinite 
generalized Laplacian, for states in real as well as complex 
Hilbert space sections (2.2 and 2.6). 

(ii) We have shown that projectors involving states in the stan- 
dard basis are associated with the edges of the graph (theo- 
rems 2.2.7 and 2.6.4) 

(iii) We have given graphical criteria for a state being pure. In 
particular, we have shown that a pure state must have a graph 
which is a clique plus isolated vertices (theorem 2.2.3, 2.2.4, 
remark 2.6.2) 

(iv) For states in a real Hilbert space, we have given an algorithm 
to construct graph corresponding to a convex combination of 
density matrices, in terms of the graphs of these matrices 
algorithm 2.2.9. 

(v) We have defined a modified tensor product of two graphs 
in terms of the graph operators C,r],J\f,Q and obtained the 
properties of these operators sections (2.4.2, 2.6.2). We have 
shown that this product is associative and distributive with 
respect to the disjoint edge union of graphs (corollary 2.4.7, 
remark 2.6.11). 
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(vi) We have proved that the density matrix of the modified ten- 
sor product of two graphs is the tensor product matrices of 
the factors, (theorem 2.4.5, 2.6.10 ). For density matrices, we 
show that a convex combination of the products of density 
matrices has a graph which is the edge union of the modified 
tensor products of the graphs for these matrices (corollary 
2.4.11, 2.6.6). Thus we can code werner's definition of sepa- 
rability in terms of graphs. 

(vii) We have generalized the separability criterion given by S. L. 
Braunstein, S. Ghosh, T. Mansour, S. Severini, R.C. Wilson 
[76] to the real density matrices having graphs without loops 
(lemma 2.4.16). 

(viii) We have found the quantum superoperators corresponding 
to the basic operations on graphs, namely addition and dele- 
tion of edges and vertices, it is straightforward to see that all 
quantum operations on states result in the addition / dele- 
tion of edges and/ or vertices, or redistribution of weights. 
However, addition / deletion of edges / vertices correspond to 
quantum operations which are irreversible, in general. Hence 
it seems to be difficult to encode a unitary operator, which 
has to be reversible, in terms of the operations on graphs. 
Further, graphs do not offer much advantage for quantum op- 
erations which only redistribute the weights, without changing 
the topology of the graph, as in this case the graph is nothing 
more than a clumsy way of writing the density matrix. 

(ix) In section 2.6, we generalize the results obtained in sections 
2.2 - 2.4, to quantum states in a complex Hilbert space, that 
is, to the density matrices with complex off-diagonal elements. 
In fact, all the results previous to section 2.5 go over to the 
complex case, except lemma 2.4.16. Many of these results 
have been explicitly dealt with (e.g. theorem 2.6.4, remark 
2.6.2, section 2.6.2 etc). We have also given rules to associate 
a graph to a general Hermitian operator. We believe that any 
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further advance in the theory reported in this chapter will 
prominently involve graph operators and graphs associated 
with operators. 

(x) Finally, we have given several graphical criteria for the pos- 
itive semidefiniteness of the generalized Laplacian associated 
with a graph. Note that by lemma 2.7.3 and observation 2.7.4 
all graphs with complex weights, either without loops or with 
positive weighted loops, have positive semidefinite generalized 
Laplacians. This characterizes a large class of graphs coding 
quantum states. 

This chapter is essentially a generalization of the work by Braunstein, 
Ghosh and Severini |[77| in which the idea of coding quantum mechanics of 
multipartite quantum systems in terms of graphs was implemented. The 
motivation in both Braunstein, Ghosh, Severini and this chapter is to ex- 
plore the possibility of facilitating the understanding of mulipartite and 
mixed state bipartite entanglement using graphs and various operations 
on them. Whether such a goal can be reached is too early to say. In order 
to code arbitrary quantum states and observables in terms of graphs, we 
have to deal with weighted graphs with real or complex weights. Unfor- 
tunately, a mathematical theory of such graphs is lacking. Many results 
pertaining to simple graphs are not available for such weighted graphs. We 
hope that, through the need of understanding entanglement and related is- 
sues the mathematical structure of weighted graphs gets richer and in turn 
gives a feedback to our understanding of entanglement. 



Chapter 3 



On the degree conjecture for 
separability of multipartite quantum 
states 

Models are to be used, not believed. 
H. Theil ( Principles of Econometrics) 

3.1 Introduction 

In chapter 2, we have given a formulation to use graphs to address the sep- 
arabihty and related problems pertaining to multipartite quantum states. 
In this chapter, we use this formulation to settle a conjecture due to Braun- 
stein et al. [73, [77] regarding the separability of quantum states. 

Braunstein et al. [76] have made a conjecture, called degree conjecture, 
for the separability of multipartite quantum states. The conjecture states 
that a multipartite quantum state is separable if and only if the degree 
matrix of the graph associated with the state is equal to the degree matrix 
of the partial transpose (with respect to a subsystem, see below) of this 
graph. We call this to be the strong form of the conjecture. In its weak 
version, it requires only the necessity, that is, if the state is separable the 
corresponding degree matrices match. We prove the strong version for a 
m-partite pure state (section 3.2) and give a polynomial-time algorithm 
for factorization of a m-partite pure state (section 3.3). We show that the 
conjecture fails, in general, for mixed states (section 3.4). Finally, we prove 
the weak version of the conjecture for a class of multipartite mixed states 
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(section 3.5). 

3.2 The Separability Criterion and its Proof 

For the definition of a weighted graph with real or complex weights, denoted 
{G, a), we refer the reader to chapter 2. We denote by V{G, a) and E{G, a) 
the vertex set and the edge set respectively of a weighted graph. The 
degree of a vertex v G V{G^a) is denoted by Dy and ^{G,a) = X^^^w is 
the degree sum of the graph. The degree matrix of a weighted graph is 
denoted by A(G, a) (see chapter 2). The combinatorial Laplacian of the 
weighted graph is denoted L(G', a) and the generalized Laplacian of (G, a) 
is denoted Q{G, a) (see chapter 2). If the generalized Laplacian of a graph 
(G, a) is positive semidefinite, we can define the density matrix of the graph 
(G, a) as a{G, a) = -^^^^Q{G, a). Conversely, given a density matrix, we can 

assign a graph to it (see chapter 2). A real weighted graph gets assigned 
to a density matrix with all real elements, while a complex weighted graph 
is assigned to a density matrix with complex off-diagonal elements. The 
vertices of the graph are labeled by the elements of the standard basis in 
the state space of the multipartite quantum system which is used to set 
up the density matrix. Also, in what follows, we use the definition and 
properties of the modified tensor product of graphs proved in chapter 2. 

Let (G, a) be the graph corresponding to a m-partite pure state in d = 
did2 • • • djn dimensional Hilbert space H = H'^'^ 'S)TC^^'^ • • • 07-^^", where di 
is the dimension of 7{^% 1 < i < m. Each vertex of (G, a) is labeled by an 
m tuple {viV2 ■ ■ ■ Vjn) where I < Vi < di, z = 1, 2, • • • , m. In other words, 
we set up (G, a) using the standard basis in 7{. Since (G, a) is the graph 
of a pure state, it must be a clique on some subset of V(G, a), all vertices 
not belonging to this subset being isolated (see chapter 2). We divide the 
m parts of the system in two nonempty disjoint subsets (partitions) whose 
union makes up the whole system. We call them s and where t is the 
complement of s in the set of all parts of the system. That is, s and t are the 
nonempty subsets of {1, 2, • • • , m}, sUt = {1, 2, • • • , m}, and sOt = 0. This 
corresponds to H = H^'^ ^n^^\ where H^'^ = ^ n^^z • • • W^^s ^ 
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and HS^^ = TC^^i (g) l-f^^^ ^ H'^js • • • (g) TC^^t^ with {n, ^2, • • • , ^s} = s, and 
{jiij2, • • • ^jt} = t- As a result of this division, we can divide the m tuple 
V = {viV2 ■ ■ ■ Vm) into the corresponding partitions (strings) which we call 
Vs, Vt. Thus, ^5 = Vi^Vi^ ' ' 'Vi^ and Vt = Vj^Vj^ ■ ■ ■ Vj^. We can label each 
vertex equivalently by {vg^vt). We call Vg The s part and Vt the t part of 
the vertex label. For example, consider a four partite system whose parts 
are labeled 1, 2, 3, 4. Let s = {1, 4}, t = {2, 3}. Then, a vertex label (1122) 
can be written as (12, 12). A vertex label {V1V2VSV4) becomes {viV4,V2Vs)- 

3.2.1 Partial transpose with partition s. 

This is a graph operator denoted Tg operating on E{G, a) which we define 
separately for the graphs with real and complex weights. 
Definition 3.2.1 : Let (G, a) be a graph with real weights. The operator 
Tg is defined as follows: 

Ts : (G,a) 1 — > {G^%a'), (G, a) 3 {{vs.Vt), {ws.Wt)} 

I — > {{ws,vt)jivsjWt)} e {G^%a), 

with 

a{{{ws, Vt), {vs, wt)}) = a{{{vs, Vt), (w,, Wt)}), 

that is, a'[Tse) = a(e). Note that, in general, E{G,a) is not closed under 
Tg. The operator Tf for the partition t giving the complement of s in the 
m-partite system is defined in the same way. Note that Tg = Tt. 
Definition 3.2.2 : Let (G, a) be a graph with complex weights. The 
operator Tg is defined as follows. 

Ts : (G, a) 1 — > (G^%a'), (G, a) 9 {{vs.Vt), (wg^wt)} 

I — > {{ws,vt),{vs,wt)} e {G^%a'), 

with 

\a'{{{ws,vt),{vs,wt)}) \ = \a{{{vs,vt),{ws,wt)})\, 

that is, |a'(r,e)| = |a(e)|. Note that, in general, E{G,a) is not closed under 
Tg. The operator Tt for the partition t giving the complement of s in the 
m-partite system is defined in the same way. Note that Tg = Tf. 
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Definition 3.2.3 : Partial transpose of (G, a) with real or complex weights 
with respect to partition s, denoted {G'^%a'), is the graph obtained by 
acting Ts on E{G,a). Note that V{G,a) = V{G^%a'), but in general 
E{G,a) E{G^%a'). 

An edge joining the vertices, whose labels have either the same s part 
or the same t part or both are fixed points of Tg. We have 

Ts{{Vs,Vt),{Vs,Wt)} = {{Vs,Vt),{Vs,Wt)}, 

with the condition on weights automatically satisfied. Similarly 

Ts{e) = Ts{{Vs,Vt),{Ws,Vt)} = {{Ws,Vt),{Vs,Vt)} = {{Vs,Vt),{Ws,Vt)}, 

with the condition on weights automatically satisfied. Note that in the case 
of complex weighted graphs the action of on {(f s, ft), (ws, Vt)} changes its 
orientation and hence a'{Tg{e)) = a*(e). However, by the definition of T,, 
this edge is still invariant under T,. Further, the definition of the operator 
Tg leaves the phase of a'(Ts{e)), e G E{G, a), as a free parameter. We shall 
use this freedom later in fixing the phases of the graphs corresponding to 
the factors in the tensor product decomposition of a density matrix. 

If both s and t parts of two vertices of e G E{G, a) are the same, then 
both vertices are identical and we have a loop, which is obviously preserved 
under Tg. Thus, Tg divides E{G, a) into two partitions, one containing all 
fixed points of T^, that is, edges with same s or t part and all loops, which 
we call J- set and the other containing the remaining edges which we call 
C set. In other words, (G, a) is the disjoint edge union of the two spanning 
subgraphs corresponding to the set and the C set. 

Note that is the identity operator (Tg = 1). Thus, Tg is its own 
inverse and is one to one and onto. 

Lemma 3.2.4 : Let (C, a) be a graph of a pure state in the Hilbert space 
n = W^' ^H'^^ ^ ■ ■ ■ ^ TC^"'. Let {G'^%a') be the partial transpose of {G, a) 
with respect to a partition s as defined above. Then, E{G, a) is closed under 
Ts, E{G,a) = E{G^%a'), if and only ?/A(G,a) = A{G^%a'). 

We emphasize that the closure of E{G, a) under means, for every 
e G E{G,a), Ts{e) = e' G E{G,a), and a(e') = a(e) or |a(e')| = |a(e)| 
as appropriate. 
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Proof : Only if part : We are given that E{G,a) is closed under T,. We 
divide E{G^ a) into two partitions C and J- as above. Note that if E{G, a) 
is closed under Tg then the sets C and are separately closed under Tg. 
Consider now the set of edges incident on a vertex in V{G,a). The edges 
in this set which belong to set are not shifted by Tg. Since E{G^a) is 
closed under and T^^ = 1, every incident edge belonging to C is the image 
of an edge in C with same weight ( or same absolute value for the weight) 
under the action of Tg on C. Thus, the degree of each vertex is preserved 
under the action of T, on E{G^a)^ for both, real and complex weighted 
{G,a) (see chapter 2), so that A(G, a) = A{G^%a'). 

If part : We are given A(G, a) = A{G'^% a'). The edges in E{G, a) belong- 
ing to J- remain in E{G, a) under the action of Tg. Now suppose that the 
set C is not closed under Tg. Since Tg is its own inverse, the edge e for 
which Tg{e) ^ E{G,a) cannot be the image of any other edge in E{G,a) 
under Tg. Therefore, the degree of the end vertices of e is changed under 
the action of Tg. This contradicts the assumption A(G, a) = A{G^% a'). ■ 

Lemma 3.2.5 : Let (G^a) be the graph of a pure state in the Hilbert 
space H, = TV^^ <S) 'H'^'^ • • • T^'^" of a m-partite quantum system. Then, 
(G, a) = (Gs, h) □ (Gt, c), where {Gg, b) is the graph of a pure state in the 
Hilbert space made up of s factors "H^*^ = TC^'^ <S) TC^'^ <S) H"^'^ <S) ■ ■ ■ <S) H'^" 
and {Gt^c) is the graph of a pure state in the Hilbert space made up of 
t = m — s factors H^^^ = TC^^^ ® H'^^^ TC^^^ ® ■ ■ ■ ® TC^^t , if and only if the 
edge set E{G, a) is closed under Tg, s = {zi, ^2, "^3, • ' ' 5 is}- 
Proof : Only if part : Case I : Graphs with real weights. We are given 

(G,a) = (G„6)H(Gi,c) 

= CiGg, b) (g) Ct){Gu c) + C{Gg, b) N{Gt, c) + N{Gg, b) ® C{Gu c) 

+n{Gs,b)®n{Gt,c) (3.1) 

Using the definition of the operators jC, r],J\f, Q, (see chapter 2) and that 
of the tensor product of graphs, we can make the following observations. 
The second term is a spanning subgraph of (G, a) each of whose edges 
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has a common t part and hence is fixed point of Tg. The third term is 
a spanning subgraph of (G, a) each of whose edges has a common s part, 
so that each edge is a fixed point of Tg. The fourth term is a spanning 
subgraph of (G, a) which contains only loops all of which are fixed points 
of Tj. Again, from the definition of the tensor product, we see that in 
the first term, any edge {vs^Wg} in C{Gs,b) with weight b{{vs,w,s}) and 
any {vt^wt} in jCr]{Gt,c) with weight c{{vt,wt}) gives us, under the tensor 
product, two edges {(fs, Vt), {wg, Wt)} and {{vg, Wt), (ws, Vt)} with the same 
weight 6({fs, Ws})c{{vt, w^t}), which are the images of each other under Tg. 
This proves that E{G, a) is closed under Tg. 

Case II : Graphs with complex weights . We are given 

(G, a) = (G„ 6) □ (Gi, c) 

= CiGg, b) (g) C{Gu c) + >C(G„ b) Af{Gu c) + A/'(G„ b) C{Gu c) 

^{n{Gs,b)0n{Gt,c)U2AfC{Gg,b) ®N'Cr]{Guc)} (3.2) 

Note that the first three terms are similar to those in Eq. (3.1) and the ar- 
guments corresponding to these terms in the paragraph following Eq.(3.1) 
apply, except that we require |a(rse)| = |a(e)|. Again, fourth term cor- 
respond to graph with loops (and no edges) which are fixed points of Tg. 
This proves that E{G,a) is closed under Tg. 

If part : We begin by noting that the graph (G, a) has the structure of a 
clique and isolated vertices, with |y(G, a)| = d = did2 ■ ■ ■ dm^ where di is 
the dimension of l-C^', 1 < i < m, because it is the graph of a pure state 
in 7i = Ti"^' Ti"^^ ■ ■ • TV^"". Let {K^, a) denote the chque in (G, a) 
and Vfc(G,a) be the set of vertices on the chque. Let \Vk{G,a)\ = n. We 
are given that E{G,a) is closed under Tg. Note that all loops are on the 
clique and no loops are on the isolated vertices. Consider a vertex [vg, Vt) 
on {Kni a). Let q denote the number of vertices in {Kn-, a) having the same 
s part as {vg,vt) and p denote the number of vertices on {Kn,a) with the 
same t part as {vg^vt). We note that p and q are the same for all vertices 
on {Kn, a), otherwise the set C is not closed under Tg. We draw {Kn, a) as 
a lattice of p rows and q columns, such that all vertices in one row have 
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common s part and all vertices in one column have common t part. Since 
{Kn, a) is a complete graph, from figure 3.1, we see that any vertex (v^^ Vf) 
has (l? — 1) + — 1) neighbors giving edges in the T set and {p —l){q— 1) 
neighbors giving edges in the C set. Since {Kn^a) is complete {ys-,vt) has 
n — 1 neighbors giving n = pq. 
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Figure 3.1: Every row contains vertices with same s part, and every column contains 
vertices with same t part. In a vertex label, first number stands for the s part and the 
second for the t part. For example, the edges between the vertex 34 and all the vertices 
in the 3rd row and 4th column are in the set JF, while the edges between vertex 34 and 
all the vertices in the four blocks obtained by deleting third row and fourth column are 

in the set C. 

Now consider C set on (Kn, a). From the definition of the tensor product 
of weighted graphs (see chapter 2), we can factorize each pair {{vg, Vt){ws, Wt)}, 
{{ws,Vt), {vs,Wt)} in the set C as the tensor product of two edges {vs,Ws} 
and {vt, Wt} with weights b' and c' satisfying 



a{{{vs, Vt), {ws, Wt)}) = b'{{vs, Ws}) ■ c{{vt, Wt}) = a{{{ws, Vt), {vg, Wt)}) 
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or 

\a{{{vs,vt),{ws,wt)})\ = \b'{{vs,Ws})\-\c'{{vt,Wt})\ = \a{{{ws,Vt), {vs,Wt)})\ 

Writing each pair {e, Tge} in C set in this way and taking the disjoint 
edge union (see chapter 2) in all these tensor products, we get {Gc^a) = 
{G'i,b') {G2,c') where {Gc,a) is the spanning subgraph on {Kn,a) cor- 
responding to set C. {G'l^b') and (G2,c') are graphs on p and q vertices, 
respectively. Again, from the definition of the tensor product of weighted 
graphs, we know that isolated vertices in the factors produce isolated ver- 
tices in the product. Therefore, we can add {dg — p) isolated vertices to 
{G'l^h') and {dt — q) isolated vertices to (G'2,c'), where ds = di^di^ ' ' ' di^ 
and dt = dj^dj^ ■ ■ ■ dj^. We call this graphs C{G'g, h') and C{G[, c') (the op- 
erator C defined in chapter 2 removes loops from a graph). The tensor 
product jC{G'g,b') jC{G[,c') gives the spanning subgraph of {G,a) corre- 
sponding to set C. Now consider a row in figure 3.1 containing vertices 
with common s part say Vg. This row generates q{q — l)/2 edges of the 
form {{vg, Vt), {vg, Wt)} all in the set J-. By the definition of the Cartesian 
product of the weighted graphs (see chapter 2), each of these edges is the 
Cartesian product of the vertex Vg in, say, (Gi, b) with the edge {-Ut, Wt}, in, 
say, (G'2,c), where a{{{vg,vt), {vg,wt)}) = dy^c{{vt,wt}). Thus, the graph 
(G2, c) is a graph of q vertices obtained by projecting each of the q{q — l)/2 
edges {{vg, Vf)^ (f^, Wf)} to {-u^, wt]^ with corresponding weight assignments, 
and thus is a complete graph on q vertices. As we vary Vg through its p 
possible values, one row corresponding to each value, the definition of the 
Cartesian product of weighted graphs generates the same (^2,0) possibly 
with different weights on edges. In exactly the same way, q columns in 
figure 3.1 generate the complete graph on p vertices (Gi, b) by employing 
the Cartesian product of weighted graphs. Noting that Cartesian product 
of an isolated vertex in (Ci, b) with an edge in {G2-,c) or vice versa gives 
an isolated vertex in the product graph, we have shown that the span- 
ning subgraph of (G, a) corresponding to the T set edges gets generated 
by (G5, 6)n(Gt, c) containing dg and dt vertices, respectively. Finally, note 
that there are no loops in (G^, 6)n(Gt, c) because loops contribute to the 
Cartesian product of weighted graphs only via the degrees X>v^ and D^^. In 
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the above analysis, question can be raised regarding the weight functions 
of the factors. Thus, more than one weight functions can generate the 
spanning subgraph corresponding to the C set, while the T set edges on 
different rows in figure 3.1 may be generated by different weight functions 
on the factors in the Cartesian product. These points will be addressed 
later in this proof. We denote by G\ and G'^ to be the graphs underlying 
and {G[,b') and similarly define G2, G2, G, G[ <S) G'2 and Gi\I]G2- 
We first note that \V{Gi)\ =p = \V{G[)\ and 1^(^2)1 = q= iViG'^)]. In 
fact, we can identify the set of vertices for Gi and G'l and for G2 and G2, 
respectively, that is, V{Gi) = V{G'i) and V{G2) = We now show 

that the identity map on V{Gi) is an automorphism taking Gi to G'l- In 
other words, Gi and G'l are identical. Consider a vertex (-Ug, Vt) in G and 
let NG{vs,Vt) denote its neighborhood in G. We denote by A^CiDGal^s? ^0 
and NG[^G'^{vs,Vt) the neighborhoods of {vs,Vt) in Gi\3G2 and G[ G2, 
respectively. In other words, NG[(g,G'2{vs,Vt) contains neighbors of {vs^vt) 
with edges in set C and A^G'inG2('^s? "^0 contains the neighbors of {vg^vt) 
with edges in set T. Clearly, NG[(g>G'^{vs,Vt) and iVciDGal^s, ^t) partition 
NG{vs,Vt). Let Vk{Gi) and Vk{G2) be the set of vertices on the chque in 
Gi and G2, respectively. Consider jl67j 

NG[®G'^iVs:Vt) = VxiG) - {{Vs:Vt)} - iVCiDGs (^^s, ^i) 

= Vk{G) - {v,} X {vt} - {{v,} X NgM U NgMs) X {vt}} 

= Vk{Gi) X Vk{G2) - {vs} X {vt} - {{v,} X NG,{vt) U NgM x {vt}} 

= {Ng,{v,) U {v,}} X {NG,{vt) U {vt}} - {v,} X {vt} - {{vs} X NG,{vt) U 
NgAvs) X {vt}} 

= Ng,{vs) X NG,{vt) U Ng,{vs) X {vt} U {v,} X NG,{vt) U {v,} X - 
{v,} X X NG,{vt) U iVGi(^.) x {vt}} 

= iVGi(^,s) X iVcaK) = NG,^G2{Vs,Vt). 

Therefore [W], iVcil^^.) x NG'^{vt) = Ng,{v,) x TVg^K) 



or iVG"i(^s) = A^Gi(^s) and NG'^{vt) = Ng^M. 

Therefore, Gi and G[ and G2 and G2 are identical. In particular, Gi 
and G2 are cliques. Gs = G'^ and Gt = are cliques plus isolated vertices. 
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Thus, we get 

(G, a) = CiGs, h') ® C{Gu d) + (G,, h)U{Gu c) + {G\ a), (3.3) 

where (C", a) is the graph obtained from (G, a) by removing all edges and 
keeping loops. 

The only remaining gap is to show that a consistent assignment of 
weights to the factors Gs and Gt is possible so as to express (G, a) as 
modified tensor product. To get the required weight assignments, we use 
the requirement that both the factors in the modified tensor product must 
correspond to pure states. Indeed, we know that both Gg and Gt have the 
form of clique plus isolated vertices as required for them to represent pure 
states. We know that the graph (G, a) corresponds to pure state. There- 
fore, its weight function a must satisfy (see chapter 2), assuming (G, a) to 
he a real weighted graphs 

E + 2 E = f(G,a)- (3-4) 

{Vs,Vt)&V{G,a) e£E{G,a)) 

Splitting £'(G, a) into C and sets and using the definitions of the 
tensor and Cartesian products of weighted graphs, we get [note that a 
paired label is for a vertex in V{G^ a) and single labels with suffix s and t 
are vertices in V{Gs, b) and V{Gt, c), respectively] 

Ks,vt)+'^ ia^i{ivs:Vt):iws,Wt)})Wi{ivs,Wt),iws,Vt)})) 

{Vs,Vt) {Vs,Vt) {Ws,Wt) 

{Vs,Vt) {Vs,Wt) {Vs,Vt) {Ws,Vt) 

Since E{G, a) is closed under T^, we can write 

{Vs,Vt) {Vs,Vt) {Ws,Wt) 
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{Vs,Vt) {Vs,Wt) {Vs,Vt) iWs,Vt) 

,2 



%,a)- (3.5) 

Using splitting of the weight functions in C set and set, we get 



+2E<E E «^(K»»}+2EalE E ft'iK, = 



Vt VJt Vt Vs Vis 

Wt^Vt Ws^Vs 



(3.6) 

Since the graphs (Gg, 6), (G^, c), (G^, 6'), and (G^, c') correspond to pure 
states, b',c',b and c must satisfy 

E< +2E E = ok,!.) (3-7) 

Vs Vs V)s 

Ws^Vs 

E < + 2 E E = "h.M) (3-8) 

Wsj^Vs 

E "S. + 2 E E = (3-9) 

Wt^Vt 

E + 2 E E = (3-10) 

Vt Vt w;* 

Wtf'Vt 

We see that Eqs. (3.7)-(3.10) are consistent with Eq. (3.6) provided 

(i) '0{vs,vt) = '^vs^vt foi" all ('^^sj'^^i) ^ V{G,a) and consequently 
t'(G,a) = ^{Gs,b)^iGt,c) and 

(ii) 6^({?;s, u;J) = b'^{{vs, Ws}) and c^d?;^, wt}) = c'^i{vu wt}). 
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We first fix a vertex {vs^vt) G V{G^a) and obtain its degree ^{v„vt)- 
Summing tiie edges witii tiie same s part, we get 

a{{{vs,vt),ivs,wt)}) =dy^ c{{vt,wt}). (3.11) 

Adding edges with the same t part, we have 

Y a[{{ys,vt),{ws,vt)]) = 'i)v, Y K{vs,Ws}). (3.12) 

Adding over edges in the C set, we get 

Y a{{{vs,vt),{ws,wt)}) = Y b'{{vs,Ws}) Y c{{vuwt}). (3.13) 

Adding these three terms and the weight of the loop on (fs,ft), we get 
The requirement that '0{y^^y^) = is satisfied provided 

(iii) h\{vs,Ws]) = b{{vs,Ws}) and c'{{vuWt}) = -c{{vuWt}) 
which leads to 

(iv) - a{{{vs, vt), {vs, vt)}) = dy^dy^ - b{{vs, Vs})c{{vu Vt}). 

This is satisfied provided a{{{vs,Vt), {vs.Vt)}) = b{{vs,Vs}) ■ c{{vt,Vt}). 
This requirement is consistent with {G", a) = Q,{Gs, b) <S) ^{Gt, c). 
We can write, therefore, 

(G, a) = CiGs, b) CiGt, -c) + (G„ 6)0(0^, c) + 0(G„ 6) (g) 0(Gi, c) 

= (a,6)H(G,,c). 

Now, let us deal with the case where (C, a) is a graph with complex 
weights (see chapter 2). In this case Eq. (3.4) is replaced by 

E + 2 E i«wi' = 'k«) (3-4') 

{vs,Vt)eV{G,a) eGE{G,a)) 
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and Eqs. (3.5)-(3.10) become 

{Vs,Vt) {Vs,Vt) {Ws,Wt) 

Y |a({(^s,^i),K^0})l'+2 5] Y Wi{ivs,vt),iws,vt)})\^ = 

^(G,a), (3.50 



Y <,v.)+^Y Y i^'(K,^J)i'E Y i^'(K^^})i^ 



(3.6') 



Vs Vt Vt Vs '^f 

Wt^Vt Wsi^Vs 



E< + 2E E IHK»J)P = ak*)' (3.7') 



Vs Vs Ws 

WsT^Vs 



E< + 2E E !'''({''.. = (3-8') 



Vs Vs 

Ws^^Vs 



+ 2E E l'=({"*."'«})P = »^o.,c). (3.9') 



Vt Vt Wt 

wtTVt 
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E< + 2E E I' = ^h..<^r (3-10') 

Vt Vt wt 

WtT^Vt 

We see that Eqs. (3.7'), (3.8'), (3.9'), and (3.10') are consistent with Eq. 
(3.6') provided 

(v) ^{vs,vt) = ^vs^vt for all {vs, Vt) e V{G, a) and consequently 
f (G,a) = ^{G,,b)^{Gt,c) and 

(vi) \b{{v,,w,})\' = \b'i{vs,w,})\'^nd\ci{vuwt})\' = \c'{{vt,wt})\\ 

Eqs. (3.11)-(3.13) become 



{Vs,Wt) 
Wt^Vt 



Wt 
Wt^Vt 



J2 \a{{{v„vt),{w,,vt)})\ = d,, J2 H{ys,w,})\, (3.12') 



{Vs,Wt) 
Ws^Vs 



Ws 



^ \a{{{vs,vt),{ws,wt)})\= ^ \b\{vs,Ws})\ ^ \c{{vuwt})\. 

Wt^Vt 

(3.13') 

Adding these three terms and the weight of the loop on (vg, Vt), we get 
?)(^^^^^). The requirement that '0(^jj^^vt) = '^Vs^vt is satisfied provided 

(vii) \b'{{vs,Ws})\ = \b{{vs,Ws})\ and \c'{{vuwt})\ = \c{{vuWt})\ 
and 



3.2 The Separability Criterion and its Proof 



113 



(viii) a{{{vs,vt),{vs,vt)}) = 

b{{vs,Vs})c{{vt,vt}) -2 \b{{vs,Ws})\ ^ \c{{vt,wt})\. 

Ws Wt 

Requirement (viii) is consistent with 

(G",a) = ^7(G„6)(g)0(G^,c)U2A/'i:(G„6){g)A/'>C77(G^,c). (3.14) 

We now choose phases of the weight functions b and c. Consider the 
edges e = {{vs.Vt), {wg^Wt)} and T^e = {{ws.Vt), {vs.wt)} in E{G,a). We 
know that |a(e)| = |a(T5e)| . Let e*^^ and e*^^ be the phases of ai = a{e) 
and a2 = a(Tse), respectively. If we require that these two edges in E{G, a) 
be produced by the tensor product of the edge {vs, Wg} in {Gs, b) with the 
edge {vt, Wt} in {Gt, c) (see figure 3.2) then the phases of weights b and c on 
the corresponding edges must be 0i = (^i + 62) /2 and 02 = (^1 — ^2)/2, 
respectively. 




Figure 3.2 



This completely fixes the weight functions b and c on (G^, b) and (Gt, c), 
respectively. We now have, for every edge e in E{G,a), the correspond- 
ing edges 61,62 in {Gs^b) and {Gt^c), respectively, such that 6 = 61 62 
and a{e) = 6(6i)c(62). Thus, we have, finally, from Eq. (3.3), (vii), and 
Eq.(3.14), 



(G,a) = £(G„6)(g)£(Gt,c) + £(G„6)(8)7V(Gt,c)+7V(G„6)(8)£(G't,c) 
+MGs, b) Q{Gt, c) U 2MC{Gs, b) AfCr]{Gt, c)} 
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= {Gs,b)B{Gt,c). (3.20 

■ 

Lemma 3.2.6 : Let a, os and at be density matrices for pure states. Then 
a = Gs ® ot if and only if {G, a) = {Gg, b) □ {Gt, c) where {G, a), {Gg, b), 
and {Gt, c) are the graphs for a, os and at, respectively. 
Proof : This lemma is identical with theorems 2.4.5, 2.6.10 in chapter 2. 

■ 

Theorem 3.2.7 : Let {G, a) be the graph of a m-partite pure state a in 
the Hilbert space W^^ ®rC^^ ® ■ ■ ■ ® W^^. Let H^'^ = H"^'^ ® H'^'^ ^ U'^h (g) 
• • • (8) Ti^'^ , {ii, ^2, • • • 5 h} = s correspond to an s set, HS^^ = TY^^^i ® Hf"^-^ ® 
H^is <S) • • • <S) H^^t, {ji,j2, • • • ,jt} = t, correspond to the t set which is the 
complement of s set in {1,2, ■■■ , m}. Then, a = Gs® ot, where os and at 
are pure states in HS^^ and 7Y^*^ with graphs {Gg, b) and {Gt, c), respectively, 
if and only if A{G, a) = A{G'^% a'). 
Proof : Using lemmas 3.2.4, 3.2.5 and 3.2.6, we have 

a = as <S) (Jt {G, a) = {Gg, b) □ {Gt, c) E{G, a) is closed un- 
der Ts <(=^ A{G,a) = A{G^%a'). A state a is entangled if A{G,a) ^ 
A(G-^% a') in every partition s and t of {1, 2, • • • , m}. ■ 



3.3 Algorithm 

While proving the if part of lemma 3.2.5, we have shown that the number 
of vertices in the cliques of the factors Gi and G2 defined there {p and 
q, respectively) are the factors of the number of vertices on the clique in 
{G,a), that is, n = pq. This means that the m-partite pure state I?/') 
corresponding to {G,a) has two factors I'i/'i) and \il^2), corresponding to Gi 
and G2, respectively, such that lives in a j9-dimensional subspace of 
'H'^h (g) 7-^^'2 (g) 7-^^'3 . . . (g) and |?/'2) lives in a g-dimensional subspace 
of TY^^^i ?^'^^2 (g) H'^h • • • (g) l-C^it, If the weighted versions of Gi and G2, 
namely, {Gg, b) and {Gt, c), can be further factorized, the dimensions of the 
corresponding subspaces will be the factors of p and q, respectively. This 
procedure can be iterated at most until the dimensions of the subspaces 
for the factors of |?/^) are the prime factors of n. Therefore, the dimension 
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of the subspaces containing the factors of are the prime factors of n or 
the products of such factors. This fact can be used to get a polynomial 
algorithm to find the full separability of in the following way. By 
full separability, we mean expressing as a product state whose further 
factorization is impossible. Denote hy pi > p2 ^ • • • > Pk the prime factors 
of n. Let di > d2 > d2, > ■ ■ ■ > djnhQ the dimensions of the Hilbert spaces 
of m parts arranged in a nonincreasing order. Let si > be the least 
integer satisfying did2 • • • dg^ ^ Pi- We implement our algorithm (theorem 
3.2.7) on partitions (s, t) with si < s < m — 1. The total number of times 
the algorithm has to run, in the worst case, is 

which is a polynomial of degree si in m. Thus, we have a polynomial 
algorithm to check separability of a m-partite system. Suppose we get the 
separability as liS^^ ® Then, the factor in liS^^ cannot be further fac- 
torized as it corresponds to the largest prime factor of n and Ti^^^ contains 
factors corresponding to j92 > J^a • • • ^ Pk- We repeat the above algorithm 
on 7t:W with P2 as the largest prime factor. Its worst case complexity is 
given by a polynomial of degree (m — s)'*% where S2 is defined like si above. 
We carry out these iterations until full separability is obtained. Thus, if 
we do not get any factorization in the first iteration, corresponding to the 
largest prime factor then the state is fully entangled, such as GHZ or 
W state. Unless the factorization carries up to m factors, the factors of 
the state contain one or more entangled states involving less than m parts. 
The total algorithm is polynomial in m. Note that, if n is prime, then all 
that is necessary is to look for some Vi common to the m tuples for all 
vertices on the clique. If, say, Vi is common, then 

with 10) G 7t:^i (g) • • • (g) Ti^'^-i Ti^'+i • • • TV^^ and \v^) e W^^ Otherwise, 
the given state |?/^) is entangled. 
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3.4 A CounterExample 

We note that theorem 3.2.7 may not apply to mixed states as the following 
example shows. Consider the bipartite separable state 

a = l/2\y,-)\y,-){y,-\{y,-\ + l/2\x,+)\x,+){x,+\{x,+\, 

where \y,—) = -^{\0)—i\l)) and |x, +) = -^(10)4- The corresponding 
density matrix in standard basis is 

2 1 + i 1 +i " 

1-i 2 2 1+i 

1-i 2 2 1+i ' 

1-i 1-i 2_ 

and the corresponding graph is shown in figure 3.3. We see that the C 
set contains only one edge for all possible partitions and hence cannot be 
closed under any Tg. We show, in section 3.5, that the degree conjecture 
applies to states with real weighted graphs without loops. Therefore, the 
above example shows that the degree conjecture does not apply to all mixed 
states. 




Figure 3.3 
A counterexample. 




3.5 Proof of Degree Criterion for A Class of Multipartite Mixed States 



117 



3.5 Proof of Degree Criterion for A Class of Multi- 
partite Mixed States 

In this section, we prove the degree criterion for the class of states whose 
graphs have no loops and have real weights. 

Theorem 3.5.1 : Let cr(G', a) be a density matrix acting on TY^^ TC^'^ 
■ • • 7i^"\ with a real weighted graph (G, a), on d = did2 ■ ■ - dm vertices, 
having no loops. If a) is separable in s, t cut where s, t is a partition 
of {1, 2, ... , m}, so that cr(G, a) = ^iPiO'i^^ ® (j'f \ where cr-*'' and af^ are 
density matrices acting on H^^^ and H^*^ with graphs {Gf \ b) and {Gf\c) 
respectively, then A{G,a) = A{G^%a'). 

Proof : Let Q{G, a) be the Laplacian of a graph (G, a) with real weights 
without loops on d vertices. For a graph without loops, a) = L{G, a) 
(see chapter 2). Let D be any d x d real diagonal matrix in the standard 
orthonormal basis i = 1,2, ... ,d, such that D and Tr{D) = 0, 

where Tr is the trace of D. This means that there is at least one negative 
entry in the diagonal of D. Denote this element by Da = bi < 0. Let 



Let \x) = l^o) + 10) = Eti(l + Then, 

{x\L{G,a)^D\x) = {ML{G,a)\iJo)^{HL{G,aM^{c^\L{G,aMo)^ 

{<j)\L{G,aM + {ilJo\D\iJo) + mD\cl)) + (0|/^|V^o) + 



Since L{G, a) is positive semidefinite, we must have, for the quadratic 
form associated with L{G, a) of a graph without loops (see chapter 2) and 



nents equal to unity, from Eq.(3.16) it follows that {i/jo\L{G, a)\i/jo) = 0. 




(3.15) 
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Also (^ol^l^o) = Tr{D) = 0. We have 

{(j)\L{G,a)\<j)) = k\L{G,a))u = k\, 

where denotes the degree of ith vertex. For a real weighted graph without 
loops, the sum of the elements in any row of its Laplacian is zero. This 
leads to 

{il;o\L{G,a)\cj)) = (0|L(G, af |V^o) = a)|^o) = 0. 

Finally, the remaining terms in Eq. (3.15) are given by 

(01^10) = b,e, 
(^o|i^|0) = b,k = imijo). 

Thus, 

{x\L{G, a) + D\x)= k\h + + 2A;6,. 
So we can then always choose a positive /c, such that 

{x\L{G,a) + D\x) <0. 

Then, it follows that 

L{G,a) + D ^ 0. (3.17) 

This expression is identical with that obtained in chapter 2 and 
For any graph (G, a) on <i = dsdt vertices, 

consider the degree condition A{G,a) = A{G^%a'). Since M^^iG^a) = 
M[G^% a'), where M is the adjacency matrix, we get, 

{L{G, a)f^ = {A{G, a) - A{G^% a')) + L{G^%a'). (3.18) 

Let 

D = A(G, a)-A{G^%a'). 

Then, D is an <i x d real diagonal matrix with respect to the orthonormal 
basis, 

\vi) = \ui) l^i), . . . , \vd) = \ud,) ^ \wd,). 
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As the degree sum of a graph is invariant under partial transpose, 

Tr{D) = Tr{A{G, a)) - Tr{A{G^% a')) = 0. 

We have two possible cases : D or D = 0. If D 0, that is, the 
degree condition is not satisfied [i.e.A{G,a) ^ A(G'-^%a')], we can write, 
via Eq.(3.17), L{(f-% a') + D ^ because {G^% a') is a real weighted graph 
without loops and D ^ with Tr{D) = 0. Using Eq.(3.18), this means 
[LiG^a))^" ^ 0. As cr(G',a) = ^-^L(G,a) (see chapter 2), a{G,a) is en- 

(G,a) 

tangled [ 155] . 



In order to test the separabifity of the m-partite state a{G, a), we have 
to apply the degree criterion to all the bipartite cuts (s, t) of m-partite 
system. This procedure may not detect the full separability of m-partite 
state jl71j . However, all the known separability criteria for multipartite 
states test separability only in bipartite cuts, and hence are not enough to 
guarantee full separability. 

In this chapter, we have settled the degree conjecture for the separability 
of multipartite quantum states. Recently, Hildebrand et al. have proved 
that the degree criterion is equivalent to the PPT criterion for separability 
[75] . However, the importance of degree conjecture ensues from the oppor- 



tunity it offers to test the strengths and limitations of a nascent approach 
to the separability problem. We see that this approach has contributed to 
test the separability of a class of multipartite mixed states (theorem 3.5.1) 
as well as to the efficient factorization of multipartite pure states. 
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Chapter 4 



Separability criterion for 
Multipartite Quantum States Based 
on The Bloch Representation of 
Density Matrices 

We must all hang together, or assuredly we shall all hang separately. 

B. Franklin 

4.1 Introduction 

In chapters 2 and 3, we presented the approach based on graphs to the 
problem of separabihty of quantum states. In this and the next chapter, 
we give another approach to the separabihty problem, based on the Bloch 
representation of a density operator. 

A A^-partite state acting on 7-^ = H'^^ TY'^^ • • • ® H^^ is separable [25] 
( or fully separable) if it can be written as a convex sum of tensor products 
of subsystem states 

N 

= = Pw>0; Y.Pro = l. (4.1) 

w w j—1 w 

A state is called k separable if we can write 
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where a^; z = 1, 2, . . . , /c are the disjoint subsets of {1,2,..., A^} and p'^^^^ 
acts on the tensor product space made up by the factors of Ti labeled by 
the members of a^. The understanding of multipartite entanglement has 
progressed by deahng with some special classes of states like the density 
operators supported on the symmetric subspace of 7i |116j . A lower bound 
on concurrence on the multipartite mixed states is obtained [ jl07j . K. Chen 
and L. Wu have given a generalized partial transposition and realignment 
criterion to detect entanglement of a multipartite quantum state jll2j . 

There are two definitions commonly used for the entanglement of mul- 
tipartite quantum states, the one from Ref. fl72] (ABLS) and the one 
introduced in flllj (DCT). In DCT, all possible partitions of N parties are 
considered and it is tested for each partition if the state is fully separable 
there or not. A state is called N partite entangled if it is not separable for 
any partition. If a state is separable for a bipartite partition, it is called 
biseparable. In ABLS, a state is called biseparable if it is a convex com- 
bination of biseparable states, possibly concerning different partitions. A 
A^-partite entangled state is one which is not biseparable. 

In this chapter we derive a necessary condition for the separability of 
multipartite quantum states for arbitrary finite dimensions of the sub- 
system Hilbert spaces and without any further restriction on them. The 
criterion is based on the Bloch representation of a multipartite quantum 
state, which has been used in previous works to characterize the separabil- 
ity of bipartite density matrix, in particular, our work is a generalization 
of de Vicente's work on bipartite systems pD]. We make use of the al- 
gebra of higher order tensors, in particular the matrization of a tensor 

[1571 mi ESI [IMl [111 El [1761 . 



The chapter is organized as follows. In section 4.2 we present the Bloch 
representation of a A/^-partite quantum state. In section 4.3 we obtain the 
main results on separability of a A^-partite quantum state. In section 4.4 
we give a sufficient condition for the separability of a 3-partite quantum 
state generalizable to the case A^ > 3. In section 4.5 we investigate our sep- 
arability criterion for mixed states, in particular, bound entangled states. 
Finally we summarize in section 4.6. 
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4.2 Bloch Representation of a A^-Partite Quantum 
State 



Bloch representation |135l 11361 11371 11381 1177j of a density operator acting 
on the Hilbert space of a (i-level quantum system is given by 



i 

Eq.(4.3) is the expansion of p in the Hilbert-Schmidt basis {Id,Xi]i = 
l,2,...,(i2 - 1} where are the traceless hermitian generators of SU{d) 
satisfying Tr{XiXj) = 26ij and are characterized by the structure constants 
of the corresponding Lie algebra, fijk, Qijk which are, respectively, com- 
pletely antisymmetric and completely symmetric. 

2 

XiXj = 2^ij^d + ifijkXk + QijkXk (4.4) 

s = (si, S2, . . . , SfP-i) in Eq.(4.3) are the vectors in constrained 
by the positive semidefiniteness of called Bloch vectors (138j . The set of 



all Bloch vectors that constitute a density operator is known as the Bloch 
vector space 5(M^^~^). The problem of determining 5(M'^^^^) where d > 3 
is still open jl36l I137j . However, for pure states {p = p^) the following 



relations hold. 



|s||2 = A/ ^^^7-^; s^Sjg^jk= {d-2)sk (4.5) 



where ||.||2 is the Euclidean norm in M'^ ^ 

It is known [HB] that B{R'^'-^) is a subset of the baU Dr{R'^'-^] 



of radius R = , which is the minimum ball containing it, and that 

the ball Dr{W^^~^) of radius r = ^J ^i^T) included in B{W^^-^). That is, 

D,(M^'-i) C 5(M^'-i) C Dr{R'^"-^) (4.6) 

In order to give the Bloch representation of a density operator acting 
on the Hilbert space C^^ C^^ . . . (g) ^f^w Qf A^-partite quantum system. 



4. Sepctrability criterion for Multipartite Quantum States Based on The 
124 Bloch Representation of Density Matrices 



we introduce following notation. We use /c, /c^ = 1, 2, • • • ) to denote 
a subsystem chosen from N subsystems, so that /c, /c^ = 1, 2, • • • ) take 
values in the set A/" = {1, 2, • • • , N}. The variables ak or for a given k 
or ki span the set of generators of SU{dk) or SU{dkJ group (Eqs.(4.3) and 
(4.4)) for the kth or kith subsystem, namely the set {Ai^, , A2j, , • • • , A^2 _i} 
for the kith subsystem. For two subsystems ki and k2 we define 

Ait^ = {Id, /d. • • • Xa,^ Id,^^, • • • /d,) 
^at^^at^ = (^^1 ^Id,^---^ K,, ® Id,,^, (S) • • • A,,^ (g) h,^^, h,) (4.7) 

where Aq,^^ and Aq,^^^ occur at the kith and /c2th places (corresponding 
to kith and /c2th subsystems respectively) in the tensor product and are 
the ak,th and a^^th generators of SU {d^^), SU{dkJ, otki = 1? 2, . . . , c?^ — 
1 and = ^^'^^ ■ ■ ■ — 1 respectively. Then we can write 



p = + E E + E E a£) Ag^) + ■ ■ ■ + 

V V t \{ki)xik2) ...xikM) ^...^ 

{A;i,A;2,-" ,A;m} Q:fciQ!fc2"""feM 

^ ^ ^aia2---aArA^^''A^2'^ • • • A^^-*}. (4-8) 
aia2---oiN 

where s*^^-' is a Bloch vector corresponding to A;th subsystem, s^^^ = 
[sa^^a^2i which is a tensor of order one defined by 

= yrr^Agl = |rr[p,A„.], (4.9a) 

where pk is the reduced density matrix for the kth subsystem. Here 
{All, A;2, • • • 2 < M < A^, is a subset of J\f and can be chosen in 

(m) ways, contributing (^) terms in the sum Y.{k,M.-M Eq.(4.8), 
each containing a tensor of order M . The total number of terms in the 
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Bloch representation of p is 2^. We denote the tensors occurring in the 

sum Y.{k,M,-,kMY (2 < M < iV) by T^^^^^^-'M = [t^^^^^^...,J which are 
defined by 

_ dk.dk^ ■ ■ ■ 4m ^ r X (fci) Ufc^) Uku)] 



*i*P^rr[p,,,,^^,,„(A„,, ® A„,^ ■ . . A„,^,)I (4.96) 



where Pkik2...kM the reduced density matrix for the subsystem {kik2 . . . k]\j}- 
We call The tensor in last term in Eq. (4.8) T^^\ 

4.3 Separability Conditions 

Before we obtain the main results we need following definition. Through- 
out the chapter, we use the bold letter for vector and normal letter for 
components of a vector, matrix and tensor elements. 

A rank-1 tensor is a tensor that consists of the outer product of a number 
of vectors. For Mth order tensor T^^'^^ and M vectors 
this means that ti^i2,,,ij^j = u^^^^u^^^ . . . uf^^^ for all values of the indices. This 



is concisely written as T^^^^ = u^^) o u^^) o • • • o u^^^) |I7Sl [T75] 



Also, given two tensors T^^^^ and S^^^ of order M and N respectively, 
with dimensions Ii x I2 x ■ ■ ■ x Im and Ji x J2 x ■ ■ ■ x Jj^ respectively, their 
outer product is defined as |175l 1157] 



(T*- ^ o tS*- ^)iii2-.iMiii2---iN ^ ^iii2--iM^iii2---iN ('^•■^0) 

Proposition 4.3.1 : A pure A/^-partite quantum state with Bloch repre- 
sentation Eq.(4.8) is fully separable (product state) if and only if 

q-ikiMr- ,kM} _ g(/ci) Q g(fc2) o . . . o S*-^^-^^ (4-11) 

for 2 < M < A. In particular T^^^ = s^^) o s^^) o • • • o s^^) holds. Here 
{/ci, A:2, • • • , ^Af} C {1,2,..., A}, and s^^^^ is the Bloch vector of kth sub- 
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system reduced density matrix. 
Proof : Notice that Eq.(4.8) can be rewritten as 

••• + •••+ y ^ [taia2-aN ~ ^aiSa2 ■ ■ ■ ^aN]^ai ^012 ' ' ' ^aj}- (4-12) 
aia2---aN 

For full separability, the sum of all the terms apart from the first 
term must vanish. Note that for every subsystem k = 1,2, ... ,N the set 
{Id, Xi;i = 1,2, ... ,dl — l} forms an orthonormal Hilbert-Schmidt basis for 
the kth subsystem. Hence A^; xi'^hi'f^ xi'^hi'f^ ■ ■ ■ X^J; • • • ; A^^A^ 

• • • A^^'' are the vectors belonging to the orthonormal product basis of the 
Hilbert-Schmidt space of the whole A/^-partite system. By orthonormality 
of the tensor product of A's occurring in different terms, the required sum 
will vanish if and only if coefficients of each term vanish separately, that is 
if and only if 

that is, 

q-ihM,- M = s(fci) o g(A;2) o . • . o s^^^) -,2 < M < N . ■ 

In fact, the condition (4.11) for all N parts is enough to decide the 
separability of pure A/'-partite quantum states, as the following proposition 
shows. 

Proposition 4.3.1a : A pure A/'-partite quantum state with Bloch rep- 
resentation (4.8) is fully separable (product state) if and only if 

r(^)=s(i)os(2)o...os^^), 



where s^^-* is the Bloch vector of A;th subsystem reduced density matrix. 
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Proof : Suppose p is a product state p = pi <^ p2 • • • 'S) Pn- Then, 

d\d'^ . . . dj^ ^ . . \ \i 

ta,a2...aN = Tr[{pi ^ p2 ^ ■ ■ ■ ^ pN)[K, ^ ^ ' ' ' ^ Kn)\ 

did2 • • • d]\[ . , . , . , 
= Tr[{piXa,) {P2K2) • • • (S) [PNKr,)] 



2 

did2 ■ ■ ■ d]\f 



2N 

Suppose the condition holds, that is, [s'^-^) o s^^) o • • • o s^^^]a^a2...aN = 
f Then 

did2 ■ ■ ■ d]s[ I 



2N 

d\d2 ■ ■ ■ d]\[ 



-Tr[{piXa,) (p2Aa2) (8) • • • (pnKn)] 



2N 



Tr[{pi (g) /?2 • • • PAf)(Aai Aq2 (g) • • • Xaj^)] 



did2 ■ ■ ■ djSf r /\ \ AM 

= ta,a2...aN = ^ Tr[p[Xa, A^.^ • • • (g) Aa^jJ. 

The equahty 

Tr [(yOi y92 • • • Pn) {Xa, A^, • • • A«^)] = Tr[p{Xa, A^, • • • A^^)] 

is satisfied for all elements in the orthonormal basis {0^iAQ;fe}, < < 
dl — 1, {ak = for J^) where {Xa^} are the d\ — l generators of SU{dk)- 
This means that the joint probabilities obtained from the ensemble of mea- 
surements of {Xai • • • Aajv) foi" states p and p = pi® p2® • • • ® Pn are equal. 
This implies 

p = pi® p2® ■ ■ ■ ® pN- 
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Note that this criterion is easily amenable with experiments. In order 
to check it for an element of T^^'^ we have to measure the corresponding 
generators on each subsystem and then check whether the product of the 
averages equals the average of the products. 

Thus in order to check whether a given pure state is a product state 
we have to check whether T*^^^ = s*^^) o s*^^) o • • • o s^^^^ where the Bloch 
vectors s^^^, s^^^, . . . , s^^^ can be constructed from the reduced density ma- 
trices /9i,/92, ••• ,PN for subsystems 1,2, ••• , iV (s^^ = ^Tr{pkXak): k G 
{1,2,- •• ,iV}, see Eq.(4.9a)). 

In the case of mixed states we can characterize separability from the 
Bloch representation point of view as follows. 

A N -partite quantum state with Bloch representation Eq.(4-8) is fully 
separable if and only if there exist vectors Uw^ G M.^^-^ satisfying Eq.(4-5), 
and weights satisfying < Pw < 1 and Pw = \ such that 

r<'" = X:P»0£i4t', sW = 2^p„4t' (4.13a) 

w w 

and 

R 

^{fci,fc2,-,fcM} = J^Py: 0^=l uii'^ (4.136) 
w 

for2<M<N; for all subsets {h, k2, . . . , ku] C {1, 2, . . . , N}, 

where s^^-* is the Bloch vector of the mixed state density matrix for kth. 

(k) 

subsystem and Uw represent the Bloch vector of the pure state of the A;th 
subsystem contributing to the wth term in Eq. (4.1). 

This follows from proposition 4.3.1 and Eq. (4.1). However, in view 
of proposition 4.3.1a, the necessary and sufficient condition is given by 
Eq.(4.13a), so that Eq.(4.13b) can be dropped. The above result can not be 
used directly, as it amounts to rewriting Werner's definition of separability 
in a different way. However, it allows us to derive a necessary condition for 
separability for A^-partite quantum states. 
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We need some concepts in multilinear algebra. Consider a tensor T^^^ G 
^hxhx-xiN^ where Ik = dl - 1. The nth matrix unfolding of T^^) (n = 
1, 2, • • • , iV) [Wl is a matrix T^^^^ G ^/„x(/„+i/„+2.../Ar/i/2.../„-i), contains 
the element ti^i^_,,ij^ at the position with row index in = 1, 2, • • • , /„) and 
column index 



{in+l — l)^n+2^n+3 • • • ^Nhh ■ ■ ■ ^n-l + " l)^n+3^n+4 • • • ^Nhh ■ ■ ■ ^n-l 
H h(^iV-l)/l/2 . . . In-l^{il-l)l2h ■ ■ ■ In-l^{l2-'^)hh ■ ■ ■ Hn-1- 

For n = 1, we take the last term = io = i^. This ordering is called 
backward cyclic fl75j . To facilitate understanding, put points on a circle 



and label them successively by ii, ^2, • ' ' ? ^iv- The consecutive terms in the 
expression for the column index in T^^jJ'* corresponding to ti^^i^^... become 
quite apparent using this circle, for more detials see section 1.5. 

For T^^) G M^i^^2x/3 matrix unfolding T^}^ contains the elements 
ti^i^i^ (ifc = 1, 2, • • • , Ik] /c = 1, 2, 3) at the position with row number ii and 
column number equal to {12 — l)h + ^3, contains ^^1^2*3 position 
with row number 12 and column number equal to (^3 — + ii and T}. 



(3) 

contains ^^1^2*3 position with row number ^3 and column number equal 

to (^1 — l)/2 + 12- Example see section 1.5 

Note that there are possible matrix unfoldings of T^^\ The matrix 
unfolding is called the matrization of the tensor jl57l 1176] . We can now 
define the Ky Fan norm of the tensor T*^^^ (of order N) over N matrix 
unfoldings of a tensor, as 

||r(^)|U^ = max{||7;^^)|U^}, n = l,...,iV; (4.14) 
where is the Ky Fan norm of matrix r/j^^ defined as the sum of 



singular values of T^^^^ |179j . It is straightforward to check that 



defined in Eq.(4.14) satisfies all the conditions of a norm and is also uni- 
tarily invariant [801 11791 . 
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The tensors in Eq.(4.13a) are called Kruskal tensors with the restriction 
< i?w < 1, ^wPw = 1 fl42l I175j . We are interested in finding the matrix 



unfoldings and Ky Fan norms of T^^^ occurring in Eq.(4.13a). The kth 
matrix unfolding for Kruskal tensor is |176 



tI^^^ = t/(^)E(t/(^) • • • • • • U^'Y. (4.15) 



Here U^''^ = [uf ^u^^^ . . . uj^] G M^^^^; A; = 1, 2, . . . iV and is the rank 
of Kruskal tensor [ jl42l 11581 1176] . i.e. the number of terms in Eq.(4.13a). 



uf ) is a vector in and is the ith column vector in the matrix E is 
the R X R diagonal matrix, E =diag[j9i . . -Pr]. The symbol denotes the 
Khatri-Rao product of matrices [ jl76j U G and V G W'^^-^ defined as 



?7 y = [ui Vi U2 V2 ... -vji] G W"^^^ where Uj and v^, i = 
1,2, ... R are column vectors of matrices U and V respectively. Eq.(4.15) 
can be rewritten as 

tI^^^ = t/(^)E[vf)vf) . ..vfY = ^(')sy(^)" (4.16) 



where v- ;i = 1,2, . . . , R are the column vectors of the matrix 
vCkf e RiNi.-,i.-2...h^,h^,...hxR and vf ) = uf ) uf • • • 

uf^'^^ • • • uf\ Using Eq.(4.16) we can write xj-^^^ as 

= Y.P-^^^^'^S^^-^ k = l,2,...,N. (4.17) 

w—l 

Theorem 4.3.2 : If a A/^-partite quantum state of did2 . . .dN dimension 
with Bloch representation Eq.(4.8) is fully separable, then 



ll-rf'^'llA-F < y^njLA(rf,-i). (4.18) 

Proof : If the state p is separable then T^^^ has to admit a decomposition 

of the form Eq.(4.13) with ||uL^^||2 = /c = l,2,...,iV. From 

definition of KF norm of tensors, Eq.(4.14), 
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\r^''^\KF = max{\\Tl^,\\KF} ; /c = 1, . . . , iV. 



From Eq.(4.17), 

I |r(^) I = max{\ I Pw^i'H'^' I \kf} ■,k = l,...,N 



w 



< max{J2p^\\u^^^v^^^^\\KF} = max{J2p^J^U^dk{dk - \\kf} 



2^ 

w w 



where Uw^ , are unit vectors in M^fe ^ and M'^^ ^ W^^-i ^ • • • 
L^'+i-i^M'^i-i-i®. • respectively, so that \ \ui!;hl!;^^\\KF = 1 for all 

k = l,2,...,N. Using E^l^- = 1 we get ||r(^)|U^ < ^^U^AK^. 



For a subsystem we get, 

Corollary 4.3.3 : If the reduced density matrix of a subsystem consisting 
of M out of iV parts is separable then I I lifi. < ^^U^l^dk{dk - 1). 

The negation of the above condition, that is. 



is a sufficient condition of entanglement of A^-partite quantum state. This 
leads to a hierarchy of inseparability conditions which test entanglement 
in all the subsystems. 

For iV = 2 the condition \\T^^^\\kf < \J ^di{di - 1)^2(^2 - 1) has 
been shown in Ref. jHD], to be a sufficient condition for entanglement 



associated with any bipartite density matrix. Note that for A^-qubits, 
di = 2, i = 1, 2, . . . , A/^, the above criterion becomes, for a separable state, 

||r(^)|kF<i. 

Consider a qudit system Hs = ®k^i^k ^ state p, supported in the 
symmetric subspace of Tis- It is straightforward to see that all the tensors 
in the Bloch representation of p are supersymmetric, that is (see Eqs.(4.8) 
and (4.9)), ta,^a,^...au^^ = tp{au,)P{au^)-P{au,,). 2 < M < N, whcrc P is 
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any permutation over indices {a^^, ak^, • • • , afcjv/}- We have, neglecting the 
constant multipHers, 



= Tr[pk,k,-kM{Xp(a^^) (g) Ap(^,^) • • • ^PK^,))] 

where P is the appropriate permutation matrix permuting the A matri- 
ces within the tensor product |179j , P^ being the transpose of P satisfying 
P^ = P^^. In particular T^^^ is supersymmetric. All matrix unfoldings 
of a supersymmetric tensor have the same set of singular values fl57] and 



hence the same KF norm. Thus, for a A^-qudit system in a state supported 
in the symmetric subspace, it is enough to calculate the KF norm for any 
one of the N matrix unfoldings to get max{||r^^^''||i^i?}. 

4.4 A Sufficient Condition for Separability of a 3- 
Partite Quantum State 

Consider the Bloch representation of a tripartite state p acting on 7i^^ 
W^^^W^^, di<d2< ds. 



did2d^ 

Qi a2 as a\a2 



Q!2 



+ Vt A^A^^) , y-^ A(2)a(3)+ V i A^i)a(2);^(3) Mi9a) 



a\a-i OL2OLS axa^OLs 



where r, s and q are the Bloch vectors of three subsystems respec- 
tively , T^^'^^ = [^a^a^] the correlation matrix between the subsystems 
/i, z/; {/i, z/} C {1,2,3} and T^^^^ = [tQ;^a2a3] the correlation tensor among 
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three subsystems. Before stating proposition 4.4.1, we need the foUowing 
definition and result. 

Kruskal decomposition of a tensor T^^^ 



R 



n-(N) t (1) (2) {N) 

— ^ o ^ o • • • o 

is caUed completely orthogonal if {ut\ uf^) = 5ku ^ = 1, 2, • • • ,N] k,l = 



1, 2, • • • ,R jl56j . where (, ) denotes the scalar product of two vectors. If 



T^^-* has completely orthogonal Kruskal decomposition, then it is straight- 
forward to show that 

l|r(^^lk^ = E^., (4.20) 
i=i 

where R is the rank of T*^^) and ^j, j = 1,2, •• • ,R are the coefficients 
occurring in the completely orthogonal Kruskal decomposition of T^^\ 
In the proof of proposition 4.4.1, we assume that completely orthogonal 
Kruskal decomposition of T^^^ /c > 2 is available. A completely orthogonal 
Kruskal decomposition may not be available for an arbitrary tensor |156l . 



The general conditions under which the completely orthogonal Kruskal de- 
composition is possible is an open problem. We conjecture that completely 
orthogonal kruskal decomposition is available for all tensors in the Bloch 
representation of a quantum state, but we do not have a proof. As it 
stands, this issue has to be settled case by case. 

Proposition 4.4.1 : If a tripartite state p acting on 7i^^ <S) H'^'^ H'^^, 
di < d2 < d^, with Bloch representation Eq.(4.19a), where T'^^) has the 
completely orthogonal Kruskal decomposition, satisfies 



'2 (di-l) || 11 ^ / 2(d2-l) || II ^ rMck-l) 
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r T V ^M^^ V did2d3 

(4.21a) 



dnd,^ V did2d^ 



then p is separable. 
Proof : The idea of the proof is as follows. 

(i) We first decompose all the tensors in the Bloch representation of p 
as the completely orthogonal Kruskal decomposition in terms of the outer 
products of the vectors in the Bloch spaces of the subsystems (coherence 
vectors). 

(ii) We prove that we can decompose p using the Kruskal decomposi- 
tions described in (i) above, as the linear combination of separable density 
matrices, which is a convex combination if the coefficient of identity is pos- 
itive. This condition is the same as the condition stated in the proposition. 

Let T^^'"^; {p,iy} C {1,2,3} in Eq.(4.19a) have singular value decom- 
position r{/^'^> = ; with ||aj^)||2 = ||af)||2 = 1 , for 
{/i, z/} C {1,2,3} and let T in Eq. (4.19a) have the completely orthog- 
onal Kruskal decomposition T = ^jUj o v^- o w^- |176l 11421 [T80] with 
ll^jlb = 11^^112 = ||wj||2 = 1- We define 



= V . Me {1,2,3} 



so that we can rewrite 



y,,.., ^ U(d, -m. - 1) (4 22a) 



Similarly, we define 



can write 



8{di-l){d2-l){ds-l) ^^ , . . 
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If we substitute Eqs. (4.22a) and (4.22b) in p Eq.(4.19a), we get 

^ ^ ai a2 Q!3 



+ 



E E /^'^ " ^..(a!^)).,(ar))..Aff Afc) 



+ 



V aia2Q!3 i 

(4.196) 

The coherence vectors a|^^ occur in ^^(R^^"^), al^"* occur in i:),(]R^'-i), 
Uj occur in Dr(M^i"^), Vj occur in Dr-(]^'^^~"^) and wj occur in ^^■(IR^^""'') 
(see Eq.(4.6)), so that they correspond to Bloch vectors. 

We can decompose p Eq.(4.19b) as the following convex combination of 
the density matrices 

Pj 1 p'j , p"j ) p'j ; Qi , q'i , n , r[ , T^i , Ti[ ] pr , Ps , pq and ^^^-^didads; 



^ V — dj, — 2 + ft) + L y — — r 



4(d2-l)(d3-l) <, , 2(rfi-l) 



2 (rf2-l) „ 11 ^ / 2(rf3-l) || II / 2(rfi-l) || II 
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'8{di-l){d2-l){d3-l)..^.. JdM 



did2d3 did2d. 
where pj in Bloch representation is 



\\r\\KF)iTr. (4.23) 



a\ OLi a3 

Q:iQ!2 ttlttS Q!2a3 

+ E (fi.)<..(vi)«(w,)„,AWAgAg,'' 



Q!lQ!2Q!3 



ai a.2 as 

(4.24) 

Note that ||T||^^ in Eq.(4.23) is defined via Eq.(4.20), which is based 
on completely orthogonal Kruskal decomposition of T. 

The Bloch vectors, correlation matrices and correlation tensors of the 
density matrices pj , p^- , fl\ , fl'^ ] Qi , q[ , ti , , T^i , ; pr , Ps , pq 
are: 

For pj, 

— , — Vj , — Wj , Jj- — u^v^- , — u^w^- , — 



For p'-, 

= Uj , = -Yj , q^- = -Wj , T/ = -u^-vj , T/ = -u^-w^- 



For p';, 

^ J — 5 s J — Vj , q J — -Wj , 1- — -UjV_^- , 1- — UjWj- 

T-,//{2,3} ~ ~ y ^„ _ _ _ 
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For p7, 

r J — Uj , s J — Vj , q j — Wj , ij — UjVj , ij- — 

~ ~ T rrT'"{2,3} ~ ~ T 'Till ~ ~ ~ 



For 



For ^' 



I' 



Q 



;;(!) 



^'{1,3} 



= T 



^>'{2,3} 



^(2) J _ n W{1,2} _ 5(1)5(2)7 

, — u , i • — a^ a^ 



= , = 



For Ti, 



T, 



t{2,3} 



= 0, q 

= , 7:- = 0. 



~(3) 



t{1,3} _ 5(1)5(3)7 



For r/, 

„r' _ 5(1) ^r' _ A „r' _ 5(3) r^r' {1,2} _ . ^r'{l,3} _ 5(1)5(3)7 

— , _ u , q^ — -a^ , — u , — a • a^ 

^r'{2,3} _ Q ^ ^r' _ Q_ 



For TT, 

rr = 0, sf = af\ q- = af , 7^-{l,2} 

^7r{2,3} _ ~(2)-(3)T 



= , T, 



7r{l,3} 



= 



T'' = 



For tt', 

rf = , sf = -af ^ , qf = -af , Tf {I, 2} = , 7^"'^^''^ = 

^7r'{2,3} _ 5(2)~(3)T _ ^ 



For pr, 
{1,2,3}, 



di 



2{di 

% = 



I" 2 



= , q, = , r/'^'"^ = ; V{/i,z/} C 
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For ps, 

rs = , ss = V^I^A , q. = , T^^^ = ; V{//, C 
{1,2,3}, Z = 0. 

For pq, 

r, = , s, = , q, = , Ti''''^ = ; V{/., C 

{1,2,3}, % = 0. 

If we write all matrices p'-, p'-'; Qi, g'„ n, t-, tt,, tt-; pr, ps, pq 
(as we have done for pj in Eq.(4.24)) in the Bloch representation and 
substitute them in Eq.(4.23) we get p as in Eq.(4.19b). 

To understand this let us see how the first term in Eq.(4.23) adds up to 
give the last term in Eq.(4.19b). The definition of pj, pl-^ p'-^ p'-' (denoting 
the Bloch vectors by si, S2, S3, S4, ••••) can be summarized in the tabular 
form 



Table 4.1 Correspondence between the first term in Eq.(4.23) and the 

last term in Eq. (4.19b). 





Si 


S2 


S3 


S1S2 


S1S3 


S2S3 


S1S2S3 


Pl 


Uj 




Wj 


UjVj 


UjWj 


VjWj 


UjVjWj 


p', 


Uj 


-Vj 


— Wj 


— UjVj 


-UjWj 


VjWj 


UjVjWj 


p" 


— Uj 




-Wj 


— UjVj 


UjWj 


— VjWj 


UjVjWj 




-Uj 




Wj 


UjVj 


-UjWj 


-VjWj 


UjVjWj 



The contribution of each column to pj + p'j + p'- + p'-' is zero except the 
last column which reproduces the last term in Eq.(4.19b). We can get the 
contributions of each term in pj^ p'j^ p'-^ p'-' to their sum by just keeping 
track of their signs. Thus we only need the following table (dropping j) 
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Table 4.2 Contributions of various terms in fl" to their sum. 





Si 


S2 


S3 




S1S3 


S2S3 


S1S2S3 


p 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


p' 


+ 










+ 


+ 


p" 




+ 






+ 




+ 


p'" 






+ 


+ 






+ 



In the same way, the contributions of the terms involving ^, r, tt are 
obtained by using the table corresponding to table 4.2 for the bipartite 
ED]- ^, T, TT which contain tensors of order two correspond to three 



case 



2-partite subsystems 12,13 and 23 . The corresponding tables are 



Table 4.3 Contributions to q' 







S2 


S3 


S1S2 


S1S3 


S2S3 


S1S2S3 


Q 


+ 


+ 





+ 











d 









+ 












Table 4.4 Contributions to r + r' 





Si 


S2 


S3 


S1S2 


S1S3 


S2S3 


S1S2S3 


r 


+ 





+ 





+ 








t' 












+ 








Table 4.5 Contributions to tt + tt' 




Si 


S2 


S3 


S1S2 


S1S3 


S2S3 


S1S2S3 


TT 





+ 


+ 








+ 





tt' 















+ 






Tables 4.2, 4.3, 4.4, 4.5 encode the procedure to construct the possible 
separable state given in Eq.(4.23). 
We now note the following points 
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(i) If the condition (4.21a) holds, then the coefficient of the matrix 
Idid2d3 in Eq.(4.23) is positive which ensures that the decom- 
position (4.23) of p is positive semidefinite. 

(ii) By virtue of Eq.(4.6), all the coherence vectors occurring in 

p'j , p'j , p'j' ; Qi : Qi ^ n , tI , TTi , ttI ; pr , ps , pq bclong to 
the corresponding Bloch spaces. 

By (i) and ( ii) we conclude that Pj , Pj , p'j' 'i Qi •, q'i , t'^ , tt^ , tt^ 5 
Pr , Ps , Pq constitute density matrices. Further, all these matrices satisfy 
condition (4.11) so that, via proposition 4.3.1, all these matrices correspond 
to pure separable states, equal to the tensor products of their reductions. 
Therefore, they constitute density matrices and they are separable and so 
must be p. ■ 

We can generalize proposition 4.4.1 to the N-paitite case by construct- 
ing the tables successively for = 4, 5, 6, • • • . First note that the number 
of p s in the first term of Eq.(4.23) lifted to the A/'-partite case is 2^^"^. For 
N = 4 we have eight. The corresponding table is 



Table 4.6 Generalization of Table 4.1 to = 4. 





Si 


S2 


S3 


S4 




S1S3 




S2S3 


S2S4 


S3S4 


S1S2S3 


S1S2S4 




+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


p(2) 


+ 


+ 






+ 










+ 






p(3) 


+ 




+ 






+ 






+ 






+ 




+ 






+ 






+ 


+ 






+ 




p(5) 




+ 


+ 








+ 


+ 








+ 


p(6) 




+ 




+ 




+ 






+ 




+ 




pi7) 






+ 


+ 


+ 










+ 


+ 


+ 












+ 


+ 


+ 


+ 


+ 


+ 
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(Table 4-6. Continued) 



S1S3S4 


S2S3S4 


S1S2S3S4 


+ 


+ 


+ 


+ 


+ 


+ 




+ 


+ 




+ 


+ 


+ 




+ 


+ 




+ 






+ 






+ 



We see that the contribution of each column to the sum Y2i p'^^^ ^^^^ 
except the last one corresponding to the Kruskal decomposition of T*^^) 
occurring in the Bloch representation of the given state p. For general case 
of A^-partite state we construct the table for p^^\ i = 1,2, ••• , 2^~^ as 
follows. First column consists of 2^^^ plus signs followed by 2^"^ minus 
signs. Second column comprises alternating 2^^^ plus and minus signs. 
Continuing in this way upto 2^~-^ = 1 we get alternating plus and minus 
signs in the (A^ — l)th column. We set the A^th column to ensure that 
there are zero or even number of minus signs in each row. Rest of the 
columns can be constructed by appropriate multiplications. This procedure 
can be checked on table 4.6. We denote the sequence of such tables for 
TV = 2,3,4, ••• as T,, z = 2,3,4, 

The tables corresponding to (A^ — 1), (A^ — 2), 2 partite subsystems 
giving rise to the remaining terms in the Eq. (4.23), lifted to A^-partite 
obtained from r/v-i,r/v-2, ■■■,T^,T2, exactly as described in the 
proof of proposition 4.4.1. In this way we can lift Eq.(4.23) to the A^-partite 
case, with the total numbers of terms X^^q^ (^)2-^~-^~* + 1. Once this is 
done, the rest of the proof for A^-partite case follows as in proposition 4.4.1. 
Thus we have 

Proposition 4.4.1a: If a A^-partite state p acting on 7i = Ti'^^ Ti'^'^ 
• • • H'^^, di < d2 < ■ ■ ■ < dM with Bloch representation Eq.(4.8), where 
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all T'^^\ k > 2 have the completely orthogonal Kruskal decomposition, 
satisfy 



{li,iy,K} 



i,k2,-,kM} ' ' ' 

^^^Si^l|r(^'lkF<l, (4.216) 



then p is separable. ■ 
For a A^-qubit system theorem 4.3.2 and proposition 4.4.1a together imply 
Corollary 4.4.2 : Let a A^-qubit state have a Bloch representation 

and let the tensor in the second term have the completely orthogonal 
Kruskal decomposition. Then p is separable if and only if < l.B 



4.5 Examples 

We now investigate our separability criterion (4.18) for mixed states. We 
consider A^-qubit state 

pily = ^i+ma^i o<p<i (4.25) 
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where \ip) is a A^-qubit W state or GHZ state. We test for = 3,4, 5 
and 6 qubits. We get. 



Table 4.7 The values of p above which the states are entangled. 



\GHZ) 


\W) 


N 


p > 


p > 




0.35355 


0.3068 


3 


0.2 


0.3018 


4 


0.17675 


0.30225 


5 


0.1112 


0.3045 


6 



Entanglement in various partitions of W noisy state Eq.(4.25) is ob- 
tained by using (A^ — n) qubit reduced W noisy state 

(4.26) 

For = 6 and n = 2 we found that the state is entangled for 0.491 < 
p<l. 

For A^ qutrits {d = 3) we test for 

pi''Jsy = ^l + PmW (4.27) 

where \ip) = -^Ylt=i\kkk . . .) is the maximally entangled state for A^ 
qutrits. 

For N = 3 and N = 4 (qutrits) the state pl^-gy in Eq. (4.27) is entangled 

for 



0.2285 <p<l, A^ = 3 

0.2162 <j9<l, A^ = 4 (4.28) 
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The state 

Pno^sy = ^/ + ^ 0^1 (4.29) 

where = |(|112) + |123) + |214) + |234)) in the space is 
found to be entangled for 0.24152 < p < 1. 

All of the above examples involve NPT states. Now we apply our crite- 
rion to PPT entangled states for which PPT criterion is not available. 



We apply our criterion to the three qutrit bound entangled state con- 
sidered by L. Clarisse and P. Wocjan |180j . given by pc where pc 



is the chess-board state given in jl80] and is an uncorrelated ancilla. 
Our criterion detects the entanglement of this state as ||T*^^^)|| = 3.75 > 3. 
Further, the four qutrit state p = {1 — [3)pc® Pc + /5//81 considered by the 
same authors yields entanglement for < < 0.2, after tracing out either 
subsystems 1 and 2 or subsystems 3 and 4. 



Now we consider the important example of the Smolin state |171l 1121 
which is a four qubit bound entangled state given by 

1 4 

PabS'd = -J2\'l'AB){i^AB\ ® \i^hD){i^hD\ (4-30) 



4 



where 1^^^) and IV^^^?) -^^^^ states. Pabcd Bloch repre- 

sentation pfif^l = ^(/®^ + Eti (^?^) so that Corollary 4.4.2 apphes (note 
that the requirement of completely orthogonal Kruskal decomposition is 
trivially satisfied). We find for this state ||T'^'^)||xi? = 3 > 1 confirming its 
entanglement. 

Our last example is the four qubit bound entangled state due to W. Diir 



where is a 4-party (GHZ) state , Pi is the projector onto the state 
which is a product state equal to |1) for party i and |0) for the rest , and 
Pi is obtained from Pi by replacing all zeros by ones and vice versa. We 
get I IT^^^^I = 1.4 > 1 confirming the entanglement of this state. 
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4.6 Summary 

In conclusion we have presented a new criterion for separability of N partite 
quantum states based on the Bloch representation of states. This criterion 
is quite general, as it applies to all A^-partite quantum states living in 
7i = T-f^^ Ti'^'^ • • • 7i^^ , where, in general, the Hilbert space dimen- 
sions of various parts are not equal. Most of the previous such criteria had 
restricted domain of applicability like the states supported on symmetric 
subspace jll6j or, are, in general, restricted to bipartite case. In propo- 



sition 4.4.1, we have given a sufficient condition for the separability of a 
tripartite state under the condition that the tensors occurring in the Bloch 
representation of the state have completely orthogonal Kruskal decompo- 
sition. This result can be generalized to the A^-partite case. Via corollary 
4.4.2, we give a necessary and sufficient condition to test the separability 
of a class of A^-qubit states which includes A^-qubit PPT states. Smolin 
state Eq.(4.30) is an important example in this class. The key idea in our 
work is the matrization of multidimensional tensors, in particular, Kruskal 
decomposation. We have defined a new tensor norm as the maximum of 
the KF norms of all the matrix unfoldings of a tensor, which is easily com- 
puted. We have also shown that this norm can be calculated even more 
efficiently for a A^-qudit state supported in the symmetric subspace. It 
will be interesting to seek a relation of this tensor norm with other en- 
tanglement measures. Again, the entanglement measures like concurrence 
known so far are successfully applied to pure states, bipartite or multipar- 
tite, while our tensor norm can be easily computed for arbitrary A^^-partite 
quantum state. Finally, our result on full separability (proposition 4.3.1) 
of A/^-partite pure states can be easily moulded for the /c-separability of 
an A^-partite pure state. In fact it is straightforward to construct an algo- 
rithm giving the complete factorization of the A^-partite pure state (see the 
paragraph following the proof of proposition 4.3.1). It is also easy to see 
that theorem 4.3.2 can be applied to any partition of a A^-partite system 
via the Bloch representation in terms of the generators of the appropriate 
SU groups. Most important is the observation that all the tensors in the 
Bloch representation can be computed using the measured values of the 
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basis operators {Xak} so that our detectiblity criterion is experimentally 
implement able . 



Chapter 5 

Experimentally accessible geometric 
measure for entanglement in N-qubit 
pure states 



First we guess it. Then we compute the consequences of the guess to see what would be implied if the law 
WE guess is right. Then we compare the result of the computation to nature, with experiment or experience, 

COMPARE IT directly WITH OBSERVATION, TO SEE IF IT WORKS. If IT DISAGREES WITH EXPERIMENT IT IS WRONG. In THAT 
SIMPLE STATEMENT IS THE KEY TO SCIENCE. IT DOES NOT MAKE ANY DIFFERENCE HOW BEAUTIFUL YOUR GUESS IS. It DOES 
NOT MAKE ANY DIFFERENCE HOW SMART YOU ARE, WHO MADE THE GUESS, OR WHAT YOUR NAME IS-IF IT DISAGREES WITH 
EXPERIMENT IT IS WRONG. THAT IS ALL THERE IS TO IT. Richard Feynman 



5.1 Introduction 



In the previous chapter, we dealt with the problem of detecting separable 
(entangled) states based on the Bloch representation of states. In this 
chapter, we deal with the problem of quantification of entanglement of 
multipartite quantum states, that is, to find a measure of entanglement 
(see chapter 1) for these states. 

Quantification of entanglement of multipartite quantum states is fun- 
damental to the whole field of quantum information and in general, to the 
physics of multicomponent quantum systems. The principal achievements 
regarding this problem are in the setting of bipartite systems. Among 
these, one highlights Wootter's formula for the entanglement of formation 
of two qubit mixed states [33] , which still awaits a viable generalization to 
multiqubit case. Others include corresponding results for highly symmet- 
ric states [321 SH [US]- The issue of entanglement in multipartite states is 
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far more complex. Notable achievements in this area include applications 
of the relative entropy [36] , negativity [601 US] Schimidt measure [102] and 
the global entanglement measure proposed by Meyer and Wallach [27] . 

A measure of entanglement is a function on the space of states of a mul- 
tipartite system, which is invariant on individual parts. Thus a complete 
characterization of entanglement is the characterization of all such func- 
tions. Under the most general local operations assisted by classical com- 
munication (LOCC), entanglement is expected to decrease. A measure of 
entanglement that decreases under LOCC is called an entanglement mono- 
tone. On bipartite pure states the sums of the k smallest eigenvalues of the 
reduced density matrix are entanglement monotones. However, the num- 
ber of independent invariants (i.e., the entanglement measures) increase 
exponentially as the number of particles N increases and complete char- 
acterization rapidly becomes impractical. A pragmatic approach would be 
to seek a measure which is defined for any number of particles (scalable), 
which is easily calculated and which provides physically relevant informa- 
tion or equivalently which passes the tests expected of a good entanglement 
measure [1], 1181 . 

In this chapter, we present a global entanglement measure for A^-qubit 
pure states which is scalable, which passes most of the tests expected of a 
good measure and whose value for a given system can be determined ex- 
perimentally, without having a detailed prior knowledge of the state of the 
system. The measure is based on the Bloch representation of multipartite 
quantum states (see chapter 4). 

The chapter is organized as follows. In section 5.2 we give the Bloch 
representation of an A-qubit quantum state and define our measure E-r- In 
section 5.3 we compute Er for different classes of A-qubit states, namely, 
the Greenberger-Horne-Zeilinger (GHZ) and W states and their super- 
positions. In section 5.4 we prove various properties of E-r, including its 
monotonicity, expected of a good entanglement measure. In section 5.5 we 
extend E^ to A-qubit mixed states via convex roof and establish its mono- 
tonicity. In section 5.6 we introduce a related measure which is additive 
and shares all other properties with Eq-- Finally, we conclude in section 
5.7. 
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5.2 Bloch Representation of An A^-qubit State and 
The Definition of The Measure 

Consider the generators {I^dx, Cy, Cz} = {co, c"b c"2, C3} of SU{2) group 
(Pauli matrices). These Hermitian operators form a orthogonal basis (un- 
der the Hilbert-Schmidt scalar product) of the Hilbert space of operators 
acting on a single qubit state space. The times tensor product of this 
basis with itself generates a product basis of the Hilbert space of operators 
acting on the A^-qubit state space. Any A^-qubit density operator p can be 
expanded in this basis. The corresponding expansion is called the Bloch 
representation of p (see chapter 4). 

In order to give the Bloch representation of a density operator acting 
on the Hilbert space • • • of an A^-qubit quantum system, 
we introduce the following notation. We use k, ki {i = 1, 2, • • • ) to denote 
a qubit chosen from A^ qubits, so that A;, /c^ (z = 1, 2, • • • ) take values in 
the set A/" = {1, 2, • • • , A^}. The variables at or for a given k or ki span 
the set of generators of the SU (2) group for the kth or kith qubit, namely, 
the set {/fc., (7i^ , (72^., cr^k } ^i^^ qubit. For two qubits ki and k2 we 

define 

(^^4^ = (/2 ® ^2 • • • (Jak, <^h<^---<^l2) 
= (/2 /2 (S) • • • aa,^ (g) /2 (S) • • • (g) /2) 
^S;^^?^ = (/2 O /2 • • • (Ja,, /2 O • • • ® (Ja,, h) (5.1) 

where Oa^^ and (Jq.^,^ occur at the kiih and A;2th places (corresponding 
to /cith and /c2th qubits, respectively) in the tensor product and are the 
ak^ih and ak^^h generators of SU{2), (a^i = 1,2,3 and au^ = 1,2,3), 
respectively. Then we can write 

{fc}c^ Oik {A;i,fc2} "fciQ!fe2 

V V t (j^k,) {k,) {kM) ^...^ 

^Oik^Oik^-akM^ak^ ^ak2 OikM ^ ^ 
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where s*^^) is a Bloch vector (see below) corresponding to kth subsystem, 
s*-^'' = [Safc]afe=i5 which is a tensor of order 1 defined by 

Sa, = Tripai^^^] = Tripkaa^i (5.3) 

where pk is the reduced density matrix for the kth qubit. Here 

{kuk2r-- ,kM}, 1<M<N, 

is a subset of J\f and can be chosen in (^) ways, contributing (^) terms 
in the sum X^j^^ /^.^ ... ^m} Eq.(5.2), each containing a tensor of order M. 
The total number of terms in the Bloch representation of p is 2^. We 
denote the tensors occurring in the sum X^j^^ ... ^m}' ^ ^ ^ by 
j-ikiM,- M} = [tak^ak^-akM]^ wMch are defined by 



where Pkik2...kM is the reduced density matrix for the subsystem {/ci/c2 • • • ku}- 
We call The tensor in last term in Eq. (5.2) T^^\ 

From Eq.(5.4) we see that all the correlations between M out of N- 
qubits are contained in T^'^^''^^' " '^^^ and all the A^-qubit correlations are 
contained in T^^^ . If p is a A/'-qubit pure state we have 

1 ^ 

k=l {kiM} {kiM,-- ,kM} 

+ ... + lir(^)||2) = i (5.5) 

Any state p = |'?/^)('?/'| existing in a d^-dimensional Hilbert space of operators 
acting on a (i-dimensional Hilbert space of kets can be expanded in the basis 
comprising d? — 1 generators of SU(d) and the identity operator. The set 
of coefficients in this expansion, namely, {rr(pAi)}, z = 1, 2, • • • , — 1, is 
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a vector in and is the Bloch vector of p. The set of Bloch vectors and 

the set of density operators are in one-to-one correspondence with each 
other. The set of Bloch vectors for a given system forms a subspace of 
denoted B{W^ ~^). The specification of this subspace for d > 3 is an 
open problem jl36l 1137] . However, for pure states, following results are 
known |138j : 

M. = ^^ (5.6) 



where D,, and Dn are the balls of radii r = 2{d-\) ^ ~ \l 

respectively, in M'^^"^. 

We propose the following measure for an A/^-qubit pure state entangle- 
ment 

i?r(W) = (||r("'||-l), (5.7) 

where T^^^ is given by Eq.(5.4) for [M = N) in the Bloch representa- 
tion of p = |?/')('?/^|. The norm of the tensor T*-^-* appearing in definition 
(5.7) is the Hilbert-Schmidt (Euclidean) norm ||r(^)|p = (r(^),r(^)) = 
^aia2---aN ^aia2 -aAf' Throughout this chapter, by norm, we mean the Hilbert- 
Schmidt (Euclidean) norm. We comment on the normalization of E-[^{\il))) 
below. 



5.3 GHZ and W States 

Before proving various properties of Et{\iI^))^ we evaluate it for states in 
the A-qubit GHZ or W class. A general A-qubit GHZ state is given by 

1^) = ^|000 • • • 0) + - • • • 1), A>2 (5.8) 

A general element of T^^^ is given by ti^i^...i^ = (V^lcr^^ ■ ■ ■ 
CijvlV^)' ^fc = 1?2,3, /c = 1, 2, • • • , A. The nonzero elements of T^^^ are 
tii...i = 2^yp{l — p), ^33... 3 = p + (—1)^(1 — p). Other nonzero elements of 
r(^) are those with 2ka2S and (A - 2A;)cri's, /c = 0, 1, • • • , [f J , where [x\ 
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is the greatest integer less than or equal to x (e.g., for A'^ = 3, ^122 etc.). 
These are equal to (-1)^ 2y^p{l-p). This gives 



I irw I p = 4p(i - p) + b + (-1)^(1 - P)? + 4p(i - P) E 2 J • ^^-^^ 

k=l ^ ^ 

Thus we get, for £^7- (| ?/;)), 

£;r(|</'» = l|rW||-i 



1^1 

^ 4p(l - p) + b + (-1)^(1 - P)? + 4p(l - p) E j ~ ^ 



10) 



Equation (5.8), with N = 2, represents a general two-qubit entangled 
state in its Schmidt decomposition. 



Thus Eq.(5.10) gives the entanglement in a two-qubit pure state. Using 
Eq.(5.10) it is straightforward to see that Eq-di/j)) for an arbitrary two- 
qubit pure state is related to concurrence by 



where concurrence C for such a state is 2^yp{l — p). 

Figure 5.1 plots Erii^)) in Eq. (5.10) as a function of p for N = 3. For 
the A^-qubit GHZ (maximally entangled) state p = 1/2, so that 



Rn = Et{\GHZ)) = 



LfJ 

Mi + Jii + (-i)T+Egj-i (5.11: 



We see that, as a function of N ^ Er{\GH Z)) increases as a polynomial of 
degree [f J. Figure 5.2 plots Er{\GHZ)) as a function of N. Er{\GHZ)) 
increases sharply with as expected. Note that Er{\GH Z)) > for GHZ 
class of states. Whenever appropriate, we normalize the entanglement of 
an iV-qubit state |V^), ^^rdV^)), by dividing by Rn = Et{\GHZ)). 
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P 



Figure 5.1. Variation of E^rdV')) [Eq. (5.10)] for N = 3, express in units of i?3, with 

parameter p. 

The A^-qubit W state is given by 

3 

where the jth term has a single 1 at the jth bit. The state \W) = 
<^f=i^fW) is given by \W) = Ej |11 • • • 0^-1 • • • 11), N>S, and has a 

single at the jth bit. We note that \W) is locally unitarily connected to 
\W) so that their entanglements must have the same value. The general 
element of T^^) for the state p=\W){W\ is 
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N 

Figure 5.2. Variation of Rn = Er{\GHZ)) [Eq. (5.11)] with number of qubits N . 



tkh-i^ = ^{00 • • • Ij • • • 00\ai^ (8) (g) • • • (g) ai^lOO • • • Ij • • • 00) 

1 ^ 

+^ (00- ••lr--OO|(7i,(g)cri,(g)---0^7i^|OO- -00) 



Only the first term contributes to ^33.. .33 = —1. Other nonzero elements 
have the form t3...3i,3 -3ii3 -3 = jj = h-S2jS -32i3 S- 
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There are (^) elements of each of these two types, so that 




Et{\W)) = Wr^'^^W - 1 = V 1 + 4— - 1 (5.13) 



It is straightforward to check that Et{\W)) = Et{\W)) as expected. 
Note that EtHW)) > 0. _ 

Next we consider a superposition of \W) and \W) states, \ips,(i)) = 
V^|VF) + y/1 — se^^\W). It is clear that the entanglement of \ips,(f>) cannot 
depend on the relative phase (j), as |^s,(/)) is invariant under the local unitary 
transformation {|0), |1)} — ^ {|0),e*^|l)} upto an overall phase factor. As 
we shall prove below, Et is invariant under local unitary transformations. 
Figure 5.3 shows the entanglement of \ips,(j)) as a function of s, calculated 
using our measure. 

An important example of a.W state and its generalizations is the one- 
dimensional spin- 1 Heisenberg antiferromagnet, on a lattice of size N, with 
periodic boundary conditions, given by the Hamiltonian 

N 

H^ = Y^ + y,y,+i + (5.14) 

where the subscripts are mod N and X, Y, Z denote Pauli operators a^, dy, cr^, 
respectively. Hn commutes with Sz = ^Zj, so the eigenstate of is a 
superposition of basis vectors |6i • • • 6n) where s of bi • • • are ones and 
N — s are zeros for some fixed < s < N. When s = 1, the translational 
invariance of implies that the eigenstates are 

l*5^*'> = -^E«*|00---1^0---0) (5.15) 

where the jth summand has a single 1 at jth bit just like the W state and 
the wave number k = for some integer < m < — 1. The state 
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0.2 0.4 0.6 0.8 1 

s 

Figure 5.3. Variation of Er{\ips,4,)), expressed in units of Rn, witli tlie superposition 

parameter s, for N = 3. 



) is locally unitarily transformed to the W state so that it has the 
same value of Et{\W)) or Er{W)). 

For s > 2 the eigenstates of have the form 



\^lJJ,{s)) = -== ^ |00---l,-0---l,^0), (5.16) 

vis) 



where 1 occurs at ji • • • js, {ji • • • js} C J\f = {1, 2, • • • , N}, and can be 
chosen in (^) ways. We see that for \i/jn{s)), ^33.. .3 = (— 1)'^ For even A^, 
^i - i2 --23 --3 with X Vs and y 2's, corresponding to the average of x (7xS, y 
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ay's, and [N 



X 



y\ cr^'s, we get, for even x and even 



t 



1--12--23--3 






Since \iIjn{s)) is a symmetric state, any permutation of its indices does not 
change the value of an element of T*^^^ (see chapter 4), so that 

2 



(N) 



2s 

E 

x,y even 





-i 2 




Figure 5.4 shows the variation of £'r(|^iv(s))) 



(N) 



see that it is maximum at s = y, which is a characteristic of the ground 
state of Hn, as expected. Note that £^T(|V^iv(y))) for the ground state 
(s = y) rises far more rapidly than the entanglement of the A^-qubit GHZ 
state Rn = Er{\GH Z)) [Eq.(5.11)] with the number of spins (qubits) A^. 
This can be understood by noting that |'?/^7v(y)) for s = y can be written 
as a superposition of \ A^-qubit GHZ states. For example, |'?/^4(2)) can 
be written as the superposition of three four-qubit GHZ states, 
1 r 1 1 1 



1 with s. We 



IV^4(2)) 



(|0011)+|1100))+^(|0101)+|1010)) 



V2 



|1001) + |0110)) 



As A^ increases, initially £'r(|V^iv(y))) is comparable to Rn, but after A^ = 

increases very rapidly, reaching 10^ for 100 qubits. 



16 the ratio 



Rn 



Also, as A^ increases, E'r(|V^7v(s))) falls off more rapidly as s deviates from 
y. We are presently trying to understand this behavior. 

Finally, in this section, we consider the superpositions of W and GHZ 
states, 

\i^w+GHzis, 0)) = v^\GHZ) + v^r^ e'^\W), (5.17) 

also considered in [99]. For three qubits, A^ = 3, a direct calculation gives, 
for this state, 

||r(^)||' = 4s2 + 6s(l- 



11/ 



ly (o<s<i). 
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Er{\^w^GHz{s. m = ||^;,(,)) II - 1, (5.18) 

which coincides with the corresponding values oiW {s = 0) and GHZ (s = 
1) states. Note that Et{\'i1^w+ghz{s, (/)))) is independent of the phase 0, in 
contrast to the entanglement measure used in [22]. Figure 5.5 shows the 
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dependence of Erdi/Jw+GHzis, (j)))) on s. 



LU 




Figure 5.5. Variation of Erd'ipw+GHzis, (p))), expressed in units of Rn, with the 

superposition parameter s, for = 3. 



5.4 Properties of Er{\^}) 

To be a valid entanglement measure, E'rdV^)) must have the following 
properties jlS, 122] • 

(a) (i) Positivity : E-rdi^)) > for all A^-qubit pure state I?/'). 

(ii) Discriminance: Er{\'4')) = if and only if is separable 
(product) state. 
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(b) LU invariance : Er{\ip)) must be invariant under local unitary (LU) 
operations. 

(c) Monotonicity : local operators and classical communication do not 
increase the expectation value of Eq-dij)). 

We prove the above properties for E-rilip)). We also prove the following 
additional properties for Erdi/^))'. 

(d) ConUnmty : ||(|^)(^| - |0)(0|)|| ^ \e{\iI;)) - ^ 0. 

(e) Superadditivity : Er(|V^) |0)) > £^t(|V^)) + ^T(|0))• 
We need the following result, which we have proved in chapter 4. 
Proposition 5.4.0 : A pure A^^-partite quantum state is fully separable 
(product state) if and only if 

where s*^^^ is the Bloch vector of kth subsystem reduced density matrix. 
The symbol o stands for the outer product of vectors defined as follows. 

Let u(i), u(2), . . . , u(^) be vectors in W^^^-\ W^^^~\ • • • , M^^i. The outer 
product u*^^) ou*^^) o • • • ou'^^^) is a tensor of order M, (M-way array), defined 
by 

^m2--ni — • • • ^ ^ ^ CLk i, /c — i, ^, , iw. 

Proposition 5.4.1 : Let be an A-qubit pure state. Then, ||tJ^^|| = 1 
if and only if is a separable (product) state. 

Proof : By proposition 5.4.0, is separable (product) if and only if 

r(^)=sWos(2)o...os(^), 
As shown in jl58l 1156] , see also section 1.5, 

{0£i«<",0£i^«) = n£:(.<'',.">). (5.20) 
where (, ) denotes the scaler product. This immediately gives, for qubits, 

lir(^)||2^ (rw^r(^)) = nf^i(s(^),s(^)) = n,||s(^)||2 = L 

Proposition 5.4.1 immediately gives the following proposition. 
Proposition 5.4.2 Let an A-qubit pure state. Then Erdi^)) = if 
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and only if is a product state. 

Proposition 5.4.3 : Let \ip) be an A^-qubit pure state. Then ||T^^)|| > 1. 
It is instructive to show this result by direct computation of T^^^ for the 
cases of two- and three-qubit states. First, consider a general two-qubit 
state 

\^P) = ai|00) + a2|01) + asllO) + a4|ll), l^^l^ = ^' 

k 

by direct computation we get 

||r(2)||2 = 1 + 8(|a2a3| - > 1. (5.21) 

This means, via proposition 5.4.2, that is a product state if |a2a3| = 
|aia4|. Next, consider a three-qubit state in the general Schmidt form [89] 

lijj) = AolOOO) + Aie'^^llOO) + A2IIOI) + A3IIIO) + A4IIII) (5.22) 

where A^ > 0, i = 0,1,2,3,4 and Af = 1. By direct calculation of 
||T(^)|| we get 

I |r 1 1^ > 1 + 12X1X1 + 8X1X1 + 8X1X1 + 8(Ai A4 - A2A3)^ > 1. (5.23) 

Here the conditions for product state become A1A4 = A2A3 and Aq = 0. We 
now prove proposition 5.4.3 for a general A^-qubit state 

If is not a product of single-qubit states (i.e., is not A^- 
separable) then it is (A^ — A;)-separable, k = 2,3,-- - ,A^. Viewing the 
A/^-qubit system as a system comprising N — k qubits, each with Hilbert 
space of of dimension 2, and k entangled qubits with the Hilbert space of 
dimension 2^, we can apply proposition 5.4.0 to this separable system of 
A^ - A; + 1 parts in the state IV'). We get = (s^^)) o (s^^)) o • • • (s^^-^)) o 

This implies, as in proposition 5.4.1, via Eq. (5.20) and Eq. (5.6) that 
11^11^ = n»l'-+'||{.«)||^ = ^^(^^>l (d = 2'). (5.24) 

If k = N we attach an ancilla qubit in an arbitrary state |0) and apply 
proposition 5.4.0 to (A^+l)-qubit system in the state |0) where I?/') is 
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the A^-qubit entangled state. This result, combined with proposition 5.4.1, 
completes the proof. Proposition 5.4.3 immediately gives the following 
proposition. 

Proposition 5.4.4 : Er{\i^)) > 0. 

We now prove that E'r{\^JJ)) is nonincreasing under local operations 
and classical communication. Any such local action can be decomposed 
into four basic kinds of operations [22] (i) appending an ancillary system 
not entangled with the state of original system, (ii) performing a unitary 
transformation, (iii) performing measurements, and (iv) throwing away, 
i.e. tracing out, part of the system. It is clear that appending ancilla 
cannot change HT*^^)!!. We prove that E-rdi^)) does not increase under the 
remaining three local operations. 

Proposition 5.4.5 : Let Ui {i = 1,2, N) be local unitary operator 
acting on the Hilbert space of subsystems Ti''^^ (i = 1,2, ■ • • , N), respec- 
tively. Let 

P = i®l,U,)p'{®l,Ul) (5.25) 

for density operators p and p' acting on 7-^ = ®f^{h&''' and let T^^^ and 
T'^^^ denote the A^-partite correlation tensors for p and p' , respectively. 
Then 

^ ||r(^)||, so that Er{p) = Er{p'). 
Proof : Let U denote a one-qubit unitary operator; then it is straight- 
forward to show that UaJJ'^ = ^p^apcrp, where [Oap] is a real matrix 
satisfying 00^ = I = 0^0. It is an element of the rotation group 0(3). 
Now consider 

= Tr[p{(^l,U,)a,, a,, ® • • • <^^A(^^lU!)] 
= Tr{pUi(T,^Ul (g) U2(7,^Ul (g)---(g) UN(J^^Ul) 
= J2 TT{pa^,®cj^,®---®<Ja.)0'^aOfX^---Of^l 

ai---ai\f 



= (r(^)xiO(i)x2 0(2)...x^o(^) 
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where is the /c-mode product of a tensor T^^^ G m3x3x- -3 ^^iq orthog- 
onal matrix O^''^ G R^""^ [1581 [IMl [157| , see section 1.5. Therefore, 



r'(^) = r(^) X 1 0^X2 0(2)... x^.o(^) 

By proposition 1.5.3 in chapter 1, we get 

lir'(^)|l = iirw xi 0(1) X2 0(2) . . . o(^)|| = ||r(^)||. 



Proposition 5.4.6 : If a multipartite pure state is subjected to a local 
measurement on the kth qubit giving outcomes ik with probabilities and 
leaving residual A/^-qubit pure state \4>i,.), then the expected entanglement 
^ij^PikET{\4'ik)) of residual state is not greater then E-rdij)), 

J^P^^MIM) < Erm- (5-26) 

ik 

Proof : Local measurements can be expressed as the tensor product matrix 
D = L>(i)(g)L>(2)(g). . .(g)L>(^) on the expanded coherence vector T |IH2j. The 
expanded coherence vector T is the extended correlation tensor T (defined 
below) viewed as a vector in the real space of appropriate dimension. The 
extended correlation tensor T is defined by the equation 

1 ^ 

^=2^ '^^l^2■■■^N^^l®^^2®■■■®^^N^ (5.27) 

where oi^, G {/, cr^;, cr^, cr^} is the ^th local Pauli operator on the kih 
qubit (do = /) and the real coefficients Ti^i^...ij^ are the components of 
the extended correlation tensor T. Equation (5.2) and Eq.(5.27) are equiv- 
alent with 7^oo--o = 1 , ^iOO--o = sf^, • • • , '^ii2--«MOO---o = '^^^'^f^^ ' ' ' 
and = T^fJ,,^^, n,^2, ••• , ^tv 7^ 0. L'(^); /c = 1, 2, • • • n are 4 x 4 

matrices. Without losing generality, we can assume the local measure- 
ments to be positive operator valued measures (POVMs), in which case 
£)(^) = diag{l, D(^)) and the 3x3 matrix D^") is contractive, D^^^^ D^^^ < I 
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|TS2j- The local POVMs acting on an A^-qubit state p correspond to the 
map p I — > A4 {p) given by 

iii2---iN 

where L^^^ are the linear, positive, trace-preserving operators satisfying 
L-^''^Lp'' = / and \L^'^\ L^'^^^] = 0. The resulting correlation tensor of 

/-^Ik ik Ik ^ Ik ^ Ik ^ 

Ai{p) can be written as 

r'(^) = r(^) xi X2 X • • • x^. D(^) 

where D'^''^ is 3 x 3 matrix and D^^^^ D^^^ < I. 

The action of POVM on /cth qubit corresponds to the map A1fc(|^)('0|) = 
T.^,M,,pMl where M,, = / • • • • • • /, T.^,LfHf^ = I and 

[Lf^^hf^"^] = 0, with the resulting mixed state YliikPik\'i^ik){4'ik\'' where 
l^i^) is the A/^-qubit pure state which results after the the outcome ik with 



Y.Pi.Er(\4,i,)(4,i,\)^Y.P^,\\r;;P' 



probability The average entanglement of this state is 

r(^)||_l 

= 5:pJ|^(f)x,D(^)||-l 

ik 

= ^Kj|D»)T,„(|¥i))||-l 



^k Ik 



where, by proposition 1.5.2 in chapter 1, D^^^Tf^j^-jd^jj)) is the kth matrix 
unfolding (see chapters 1 and 4) of Xfc D^'^\ Therefore, from the def- 
inition of the Euclidean norm of a matrix, ||74|| = ^yTr{AA^), jl66j we 
get 

ik ik 

= ^p,,{Tr[DW^D(')r,,,(|*))Tt,,(|V'»]}^ - 1 

ik 
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ik 

= ll'^^'l|-l=£^r{W), 

because < /, and ^i^Pi^ = 1- We have also used the fact that 

Euchdean norm of a tensor equals that of any of its matrix unfoldings. ■ 
As an example, we consider the four-qubit state fl83j 

^j^CD = ^(lOOOO) + 10011) + 10101) + 10110) + 11010) + 11111)). (5.28) 



A POVM {741,742} is performed on the subsystem A, which has the form 
Ai = Uidiag{a, /3}V and A2 = U2diag{Vl - o?, ^/l^^}V. Due to LU 
invariance of HT*^^-*!! we need only consider the diagonal matrices in which 
the parameters are chosen to be a = 0.9 and (3 = 0.2. After the POVM, 
two outcomes |0i) = Ai\iIj) / ^/pi and |02) = ^2|V^)/V^ obtained, with 
the probabilities as pi = 0.5533 and p2 = 0.4467. We find 

Erim = 0.7802, ^r(|0i)) = 0.0725^, Erilh)) = 0.0436/^2- 
This gives 

ETim-[piEri\M+P2Eri\h))] =0.6641 > 0. 



This is to be contrasted with the similar calculation in jl83j , with the same 



state in Eq.(5.28) and the same POVM given above. 

Proposition 5.4.7 : Let be an A^^-qubit pure state. Let p denote 

the reduced density matrix after tracing out one qubit from the state 

Then 

1 1^(^-1) 1 1 <- 11^(^)11 
iKp II ^ ll^i^) II 

with equality only when = |0) Ix)? where \x) is the state of the qubit 
which is traced out. 

Proof : we prove this for a special case whose generalization is straight- 
forward. Let \^JJ) = a\bi ■ ■ ■ b]\i) + b\b'i ■ ■ ■ 6'^), |ap + |6p = 1. Here \bi) and 
1 6^) are the eigenstates of a^z ^ operating on the ith qubit. Now consider 
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sets of A^-fold tensor products of qubit operators {cTq,}, q; = 1, 2, 3, namely, 
S = {(Ja^ CTaa • • • CTc^}, tti • • • ttiv = 1, 2, 3. Choosing , c^at = 3 

we get (73 0-3 • • • cr^lbi • • • b^) = ±|6i • • • 67V). We can choose an opera- 
tor from S, denoted B, such that 5|6i • • • 6iv) = • • • b'j^)- If B contains 
q < N Gx operators we can replace k < q oi them by Gy operators. We 
denote the resulting tensor product operator by Bj. {Bq = B). We have, 
Bk\bi---bN) = ±{if\b'^---b'j,). Then, 

(61 • • • hnWi • • • (g) o-3|6i • • • 67v) = ±1 = (6i • • • fe'^vks • • • G2\h\ • • • 6^) 
(61 • • • hN\B\h\ • • • 6^) = ±1 = (6; • • • h'^\B\bi ■■■hN) 
{b[---h'^\Bk\hi---hN) = ±{if 
{bi---bN\Bk\b\---b'^) = ±{-if 

Now, 

tav-aN = (V'kai • • • (JaM 

= \af{bi ■ ■ ■ bNlcTar®- ■ • • • bN)+\bf{b'i ■ ■ ■ 6jvkai<^- • -(^^ajvl^i • • • ^n) 

+a*b{bi ■ ■ ■ bN\(Jai®- ■ ■'S)Gaj^\b[ ■ ■ ■ b'^)+ab*{b[ ■ ■ ■ b'^lGa^^- ■ ■'S)GaJbi • • • ^tv) 

The nonzero elements of tai-ajv ^33-3 = i I^P? = ±a*6 = 
±2\a\\b\cos{(l)a - (l)h), 



tB, = ±(itab' ± (-ita'b = { ^l^rt^ - *| I f 

^ ^ ^ ' \ ±2\a\ \b\sin{(t)a - (pb) if k is odd. 

We get Ylk=o Qk) elements with cos{(l)a — 4>b) and (2fc+i) elements 

with sin{(f)a — (t)b)- If q is odd (for the given state]'?/')) the number of cosines 
and the number of sines are equal. When q is even the number of cosines 
exceeds by 1. Finally we get 



k=0 



4\af\bfsin\^a-^b)Yl 



k=0 



2k + I, 
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Note that, using |ap + |6p = 1, it is easy to see that HTj^^ ^|| > 1, showing 
that Eq- > 0. Next we consider 

= \af\bi . . . . . . bN\ + • • • b'j,){b[ ■ ■ ■ b'j,\ + 

ab*\bi • • • bN){b[ • • • 6^1 + a*b\b'i • • • ' ' ' &a^I 

and trace out the A^th qubit to get the {N — l)-qubit reduced density 
matrix 

p = la^foi • • • bN-i){bi • • • bN-i\ + \b\^\b[ • • • b'j,_,){b[ • • • b'j,_,\ 

+ab*\bi ■ ■ ■ bN-i){b[ ■ ■ ■ b'j,_,\{bNK) + a*b\b[ ■ ■ ■ b'^_,){bi ■ ■ ■ bN-iUb'^lbN) 
Now 

tav-aN-i = Tr{p(Ja, CT^.^ (g) • • • Cr«^_J 

= \af{bi ■ ■ ■ 6iv-i|o-ai • • • cTajv-il^i • • • ^iV-l) 

+ \b?{b[ ■ ■ ■ b'j^_-^\(Ta, (8) • • • (g) (Tar^JK ■ ■ ■ 
+a*6(6i • • • bN-lWai • • • craN-i\b'i ■ ■ ■ &iv-i)(^iv|^']v) 
+a6*(6i • • • 6'7v-ikai • • • cra;v-il^l " " " &iV-l)(&ivlM- 

We have for AT — 1 tensor product operators (73 (g) (73 (8) • • • (8) (73I61 • • • 67V-1) = 
±|6i • • • bN-i). We construct the operators D and corresponding to B 
and acting on N — 1 qubits. We then get D\bi - • • 67V-1) = ±|6'j • • • fe'^v-i) 
and Dk\bi---bN-i) = ±(i)^|6'i • • • 6^_i). Now, the nonzero elements of 
r;^-^)are hs-s = ^ \b\^ 

to = ±ab*{bM\b'N) ± a*b{b'^\bN) = 2\a\\b\\{b'j,\bN)\cos{^a - (t>b - a), 
to, = ±{i)''ab''{bN\b'j^) ± (-z)V6(6'^|6^> 

_ J ±2|a||6| \{b'js^\bN)\cos{(j)a — 4>b — ot) ii k is even, 
\ zb2|a||6| |(6^|67v)|sm(0a — cj)}, — a) \i k is odd. 

Finally we get 

I l^^-'^l I' = ± + 4|ap|6|2|(6'^|6^)| W(0, - 06 - 1^) E 

fc=0 
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+A\a\^\hWh'^\hN)\hin\(l)a -(j)b-a)J2 



k=0 



q' 

2A; + 1 



where q' < q is the number of ai operators in D. Since |(&7vl^-/v)P < 1 we 
see that 

equality occurring when l^^v) = W^), in which case = |0) (g) l^^v)- It is 
straightforward, but tedious to elevate this proof for the general case 



a, 



baj^), ai = 0, 1. 



Basically we have to keep track of (2) B type of operators, where r is the 
number of terms in the expansion of \tJj), in order to obtain all nonzero 
elements of When Nth particle is traced out, the corresponding 

elements of Tp^~^^ get multiplied by the overlap amplitudes, which leads 
to the required result. ■ 
Continuity of Et'- We show that for A/'-qubit pure states — |<;z!))(0|)|| - 







Er{m-Erm 



0. 



Proof: ||(|V^)(V^|-|0)(0|)||^O 



{N) 



Butl|r(^^-Tf)| 



> 



{N)\ 



- \\r 



1^) 
(^)l 



Therefore I IT^^^ - rl^^^l 



Er{m-Er{m 







\% 



^\<i>) I 



_ 11^(^)1 



0. 







0. 



5.4.1 Entanglement of Multiple Copies of a Given State 

LU invariance. We show that Eq- for multiple copies of A^-qubit pure 
state is LU invariant. Consider a system of x A; qubits in the state 
Ix) = IV^) ® IV^) (8) • • • IV^) {k copies). It is straightforward to check that 
(see chapter 4) 



(5.29) 
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This implies, in a straightforward way, that 

ll^lx) II ~ II^W II • 

Since by proposition 5.4.5 ||'7j'^^'*|| is LU invariant, so is ||'7j'^^^||. 

Let be a A^-qubit pure state and |x) = \ip) <S) Then Eq-{\x)) is 
expected to satisfy 

EtUx)) > Mm- 



We again use the fact that 
which gives 



11^(^)11 _ 11^(^)112 

ll^lx) II ~ II^IV-) II • 
Smce||^^)||>lweget||^^)||>||^^)||or, 

et{\x)) > Mm- 

Superadditivity : We have to show, for A^qubit states {i/j) and \(/)) that 

Erm m > Erim + Erim- (5.30) 
We already know that for |x) = IV^) ^ 1^) 

||^(iV)|| _ ||^(iV)|| ||o-(iV)|| 

ll^lx) II ~ II^IV') II ll^l<^) II 
Thus Eq. (5.30) gets transformed to 

||^(iV)|| ii^Wii _ -I ^ 11^(^)11 , ll-r^^)!! _ 9 

ll^l^) II ll^l^) II 1 ^ ll^l^) II -f- ll^l^) II ^ 

which is true for ||'7j'^'^^|| > 1 and ||'7j'^^^|| > 1. ■ 
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5.4.2 Computational Considerations 

Computation or experimental determination of Et involves 3^ elements of 
T^^^ so that it increases exponentially with the number of qubits A^. How- 
ever, for many important classes of states, Er can be easily computed and 
increases only polynomially with A^. We have already computed Eq- for the 
class of A^-qubit W states, GHZ states, and their superpositions. We have 
also computed Eq- for an important physical system like one-dimensional 
Heisenberg antiferromagnet. For symmetric or antisymmetric states, T^^^ 
is supersymmetric, that is, the value of its elements are invariant under 
any permutation of its indices (see chapter 4). This reduces the problem 
to the computation of ^{N + 1)(A^ + 2) distinct elements of T''^\ which is 
quadratic in A^. 



5.4.3 Entanglement Dynamics : Grover Algorithm 

We show that Ej- can quantify the evolution of entanglement. We consider 
Grover's algorithm. The goal of Grover's algorithm is to convert the initial 
state of A^ qubits, say |0- • -0), to a state that has probability bounded 
above | of being in the state \ai - ■ ■ a^), using 

Ua\bl---bN) = {-lf'''^'^\bl---bN) 

the fewest times possible. Grover showed that this can be done with 
0{V2^) uses of Ua by starting with the state 

2^-1 



2^ n 

X={j 



\x) = //«^'|0---0), 



where 

H 



1 1 



and then iterating the transformation H®^UaH®^Ua on this state [27] 



The initial state is a product state as is the target state, but intermediate 
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states ip{k) are entangled for A; > iterations. Figure 5.6 shows the de- 
velopment of E'r{\'ip{k))) with number of iterations /c, for six qubits. The 
values of k for which Ej- vanishes are the iterations at which the probabil- 
ity of measuring |ai • • • a^v) is close to 1. Thus Eq- can be used to quantify 
the evolution of a A^-qubit entangled state. 



0.4 




20 40 60 80 100 



k 

Figure 5.6. Entanglement in Grover's algorithm for six qubits as a function of number of 

iterations. 



5.5 Extension to Mixed States 

The extension of Ej- to mixed states p can be made via the use of the 
convex roof or (hull) construction as was done for the entanglement of 
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formation [33]. We define E^ip) as a minimum over all decompositions 



p = '^iPilipiji'ipil into pure states i.e., 

Er{p)= mmT^P^Erm)- (5-31) 

The existence and uniqueness of the convex roof for Ej- is guaranteed 
because it is a continuous function on the set of pure states [42]. This 
entanglement measure is expected to satisfy conditions (a), (b) and (c) 
given in section 5.4 and is expected to be (d) convex under discarding of 
information, i.e., 

J2p^Mp^) > Er{J2p^P^)- (5-32) 

i i 

The criteria (a)-(d) above are considered to be the minimal set of re- 
quirements for any entanglement measure so that it is an entanglement 
monotone [15] . 



Evidently, criteria (a) and (b) are satisfied by E-rip) defined via the 
convex roof as they are satisfied by Ej- for pure states. Condition (d) fol- 
lows from the fact that every convex hull (roof) is a convex function filj . 
We need to prove (c), which is summarized in proposition 5.5.1. 
The proof follows from the monotonicity of E-rHip)) for pure states, that 
is, propositions 5.4.5, 5.4.6 and 5.4.7. Bennett et al. prove a version of 
proposition 5.5.1 in [32], which applies to any measure satisfying propo- 
sitions 5.4.5, 5.4.6 and 5.4.7. Thus the same proof applies to proposition 
5.5.1. However, we give it here for the sake of completeness. 
Proposition 5.5.1: If an A^-qubit mixed state p is subjected to a lo- 
cal operation on ith qubit giving outcomes ki with probabilities pk^ and 
leaving residual A^-qubit mixed state pk^, then the expected entanglement 
^^^k^Ph^TiPh) of the residual state is not greater than the entanglement 
Er{p) of the original state, 

5]pfc.^T(pfcJ < Et{p). 



(If the operation is simply throwing away part of the system, then there 
will be only one value of /c^, with unit probability.) 
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Proof : Given mixed state p there will exist some minimal-entanglement 
ensemble 

e = te, IV^)} (5.33) 

of pure states realizing p. 

For any ensemble S' realizing p, 

E{p) < Et{8'). (5.34) 

Applying the propositions 5.4.5, 5.4.6 and 5.4.7 to each pure state in the 
minimal-entanglement ensemble we get, for each j, 

Y.Pk.\jEr{Pjk.)<ETm). (5.35) 

hi 

where pj^. is the residual state if pure state is subjected to iih partite's 
operation and yields result /c^, and pj.^\j is the conditional probability of 
obtaining this outcome when the initial state is 

Note that when the the outcome ki has occurred the residual mixed 
state is described by the density matrix 

Ph = ^Pj\hPjh- (5.36) 

3 

Multiplying Eq.(5.35) by pj and summing over j gives 
By Bayes theorem 

Pj,ki = PjPki\j = PhPj\kv (5-38) 

Equation (5.37) becomes 

Y^PkmhEriPj^) < Er{p). (5.39) 

Using the bound Eq.(5.34), we get 

Y^PhE^Ph) < J^PkiJ^Pjlki^-^M < Erip). (5.40) 
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Note that any sequence of local operations comprises local operations 
drawn from the set of basic local operations (i)-(iv) above, so that proposi- 
tion 5.5.1 applies to any such sequence. Thus we can say that the expected 
entanglement of a A^-qubit system, measured by Eq-ip)-, does not increase 
under local operations. 



5.6 A Related Entanglement Measure 

We consider the following entanglement measure. Consider 
Srim = ^o^2||r(^^|| = log2[ETm + 1], 

where T^^^ is the N-w&y correlation tensor occuring in the Bloch rep- 
resentation of p = \ip){i/j\. 

Proofs of propositions 5.4.1-5.4.7 and 5.5.1 easily go through for Sril'ip))- 
We prove continuity as follows. 

Continuity of StU'^P))- We have to show, for two A/'-qubit states and 

^0. 



10), that - i0)(0i)ii ^ erim - STim 



We have ||(|^)(^| - ^ ^ - r(f)|| ^ 0. 



Bntic-^rii^iii^riiHi^riii- 

Further, whenever 1 1 |7j'^'^'*| | > 1 and |^^)^'*|| > 1 
we have |||7fg"|| - ||7f^'||| > - ^adl-Jj^' 

Thus II ^0 ^ |||^<^>||-|Clll - 0^ I'^S^dC 

iofl2(||I[^'||)H ^ |fr(W) - M|.#>»| - 0. 

However, ETili^)) has the added advantage that it is additive [ while 
Eq-{\ijj)) is superadditive] . Indeed, from section 5.4.1 we see that for k 
copies 

Similarly (g) |0)) = Srim + Srim. 
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The extension of £'r(| '?/')) to mixed states via convex roof construction is 
similar to that of E-rdij)). Thus Srili^)) has all the properties of Erdi')), 
with an additional property that Sr{\ip)) is additive, while ErUip)) is su- 
peradditive. 



5.7 Conclusion 

In conclusion, we have developed an experimentally viable entanglement 
measure for A^-qubit pure states, which passes almost all the tests for being 
a good entanglement measure. This is a global entanglement measure in 
the sense that it does not involve partitions or cuts of the system in its 
definition or calculation. This measure has quadratic computational com- 
plexity for symmetric or antisymmetric states. Computational tractability 
is not a serious problem if N is not too large, and the measure can be easily 
computed for systems comprising small number of qubits, which can have 
many important applications such as teleportation of multiqubit states, 
quantum cryptography, dense coding, distributed evaluation of functions 
1184] etc. However, finding other classes of states for which Er can be com- 



puted polynomially will be useful. It will be very interesting to seek appli- 
cations of this measure to situations like quantum phase transitions [ I20| . 
transfer of entanglement along spin chains [ jl85] . NOON states in quantum 
lithography |186j etc. Finally, we have extended our measure to the mixed 



states and established its various properties, in particular, its monotonic- 
ity. We may also note that neither its definition nor its properties depends 
in an essential way on the fact that we are dealing with qubits, so that 
this measure can be defined and applied to a general A^^-partite quantum 
system. 
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Chapter 6 



Summary and Future Directions 

This is not the end. Nor is this a beginning of the end. 
This may at most be the end of a beginning. 
-Sir Winston Churchill^ 

In this thesis, I have tried to enhance understanding of the foUowing 
two questions: 

A- Given a multipartite quantum state (possibly mixed), how to find out 
whether it is entangled or separable? (Detection of entanglement.) 
B- Given an entangled state, how to decide how much entangled it is? 
(Measure of entanglement.) 

Answers to both these questions are known for bipartite pure states. 
For multipartite states, general answers to both these questions are not 
known. Many separability criteria are proposed. Example: Generaliza- 
tions of Peres-Horodecki criterion. The genuine entanglement of pure mul- 
tipartite quantum state is established by checking whether it is entangled 
in all bipartite cuts, which can be tested using Peres-Horodecki criterion. 

For mixed states this strategy does not work because there are mixed 
states which are separable in all bipartite cuts but are genuinely entangled. 
A direct and independent detection of genuine multipartite entanglement 
is lacking. We have explored two approaches. 

I- In the first approach, we assign a weighted graph with multipartite quan- 
tum state and address the question of separability in terms of these graphs 
and various operations involving them (Chapters 2 and 3). 

II- In the second approach, we use the so called Bloch representation of 

^From his last address to the British parliament as the prime minister of U.K. ( 3rd June 1946) 
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multipartite quantum states to establish new criteria for detection of mul- 
tipartite entangled states (Chapter 4). We further give a new measure for 
entanglement in A^-qubit entangled pure state and formally extend it to 
cover A^-qubit mixed states (Chapter 5). 

In the following, I give some of the key results obtained in the this the- 
sis. 

Chapter 2. 

1- We have given rules to associate a graph with a quantum state and a 
quantum state to a graph, with a positive semidefinite generalized Lapla- 
cian, for states in real as well as complex Hilbert space. 

2- We have shown that projectors involving states in the standard basis 
are associated with the edges of the graph. 

3- We have given graphical criteria for a state being pure. In particular, 
we have shown that a pure state must have a graph which is a clique plus 
isolated vertices. 

4- We have given an algorithm to construct graph corresponding to a convex 
combination of density matrices, in terms of the graphs of these matrices. 

5- We have defined a modified tensor product of two graphs in terms of 
the graph operators £, r/. A/", Q and obtained the properties of these opera- 
tors. We have shown that this product is associative and distributive with 
respect to the disjoint edge union of graphs. 

6- We have proved that the density matrix of the modified tensor product 
of two graphs is the tensor product of the density matrices of the factors. 
For density matrices, we show that a convex combination of the products 
of density matrices has a graph which is the edge union of the modified ten- 
sor products of the graphs for these matrices. Thus we can code werner's 
definition of separability in terms of graphs. 

7- We have generalized the separability criterion given by S. L. Braunstein 
et al. to the real density matrices having graphs without loops. 

8- We have found the quantum superoperators corresponding to the basic 
operations on graphs, namely addition and deletion of edges and vertices, it 
is straightforward to see that all quantum operations on states result in the 
addition / deletion of edges and/ or vertices, or redistribution of weights. 
However, addition / deletion of edges / vertices correspond to quantum 
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operations which are irreversible, in general. Hence it seems to be difficult 
to encode a unitary operator, which has to be reversible, in terms of the 
operations on graphs. Further, graphs do not offer much advantage for 
quantum operations which only redistribute the weights, without changing 
the topology of the graph, as in this case the graph is nothing more than 
a clumsy way of writing the density matrix. 

10- Finally, we have given several graphical criteria for the positive semidef- 
initeness of the generalized Laplacian associated with a graph. This char- 
acterizes a large class of graphs coding quantum states. 
Chapter 3. 

1- We settle the so-called degree conjecture for the separability of multipar- 
tite quantum states, which are normalized graph Laplacians, first given by 
Braunstein et al. [76]. The conjecture states that a multipartite quantum 
state is separable if and only if the degree matrix of the graph associated 
with the state is equal to the degree matrix of the partial transpose of this 
graph. We call this statement to be the strong form of the conjecture. In 
its weak version, the conjecture requires only the necessity, that is, if the 
state is separable, the corresponding degree matrices match. We prove the 
strong form of the conjecture for pure multipartite quantum states using 
the modified tensor product of graphs defined in chapter 2, as both neces- 
sary and sufficient condition for separability. 

2- Based on this proof, we give a polynomial-time algorithm for completely 
factorizing any pure multipartite quantum state. By polynomial-time al- 
gorithm, we mean that the execution time of this algorithm increases as a 
polynomial in m, where m is the number of parts of the quantum system. 

3- We give a counterexample to show that the conjecture fails, in general, 
even in its weak form, for multipartite mixed states. 

4- Finally, we prove this conjecture, in its weak form, for a class of multi- 
partite mixed states, giving only a necessary condition for separability. 
Chapter 4. 

1- We give a new separability criterion, a necessary condition for separa- 
bility of A^-partite quantum states. The criterion is based on the Bloch 
representation of a A^-partite quantum state and makes use of multilinear 
algebra, in particular, the matrization of tensors. Our criterion applies to 
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arbitrary A^-partite quantum states in H = H^^ ® H'^^ • • • ® H^^^ . The 
criterion can test whether a A/'-partite state is entangled and can be apphed 
to different partitions of the A^-partite system. 

2- We provide examples that show the ability of this criterion to detect 
entanglement. We show that this criterion can detect bound entangled 
states. 

3- We prove a sufficiency condition for separability of a three-partite state, 
straightforwardly generalizable to the case > 3, under certain condition. 

4- We also give a necessary and sufficient condition for separability of a 
class of A^-qubit states which includes A^-qubit PPT states. 

Chapter 5. 

1- We present a multipartite entanglement measure for A^-qubit pure states, 
using the norm of the correlation tensor which occurs in the Bloch repre- 
sentation of the state. 

2- We compute this measure for several important classes of A^-qubit pure 
states such as GHZ states, W states and their superpositions. We compute 
this measure for interesting applications like one dimensional Heisenberg 
antiferromagnet. We use this measure to follow the entanglement dynam- 
ics of Grover's algorithm. 

3- We prove that this measure possesses almost all the properties expected 
of a good entanglement measure, including monotonicity. 

4- Finally, we extend this measure to A^-qubit mixed states via convex roof 
construction and establish its various properties, including its monotonic- 
ity. 

5- We also introduce a related measure which has all properties of the 
above measure and is also additive. 

Here are some of the interesting research problems emerging from our 
work. 

(i) The principal achievement of first two chapters, apart from giving 
a new formulation is the proof of the degree criterion for separability of 
N-partite pure states and their factorization into entangled parts. One of 
the open problems of this new formulation is to find graphical criteria for 
the non-negativity of the generalized Laplacian associated with a graph. 
As we have seen, degree criterion fails, in general, for the mixed states. 
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We could prove this criterion only for states with density matrices with 
real weighted graphs without loops. Such matrices have all elements real 
so there are no coherences [ jl87 ] and elements in every row and every col- 
umn sum up to zero |188j . Can we then use graph topology to classify 
quantum states based on separability and seek (possibly different) criteria 
for separability of different classes of states? Is it possible to code LOCC 
operation in terms of operations on graphs? It may be a good idea to 
use Jamiolkowaski isomorphism jl89j between states and quantum opera- 
tions. If we combine these two questions, we can seek the classification of 
N-partite entanglement in terms of classes of states not inter-convertible 
via SLOCC [1]. These are some of the interesting questions on the basis of 
chapters 2, 3, but we feel that its a long way to get there, if at all we can. 

(ii) There is a variety of questions emerging from chapters 4 and 5. 
It is interesting to look for lower bound on say, concurrence jll3j of three 
partite state to the violation of separability condition based on the criterion 
stated in chapter 5. Further, it will be interesting to seek a new PPT 
entangled state which is detected by the criterion in chapter 4 but not 
by any other criterion. It will be interesting to generalize the measure in 
chapter 5 to d-level systems instead of qubits. Using our measure, can 
we get tight upper or lower bounds on the entanglement of superposition 
of multipartite states jl90j ? A very interesting question is whether we 
can obtain the entanglement dynamics of a multipartite system in terms 
of our entanglement measure? For this, we will have to get the effect on 
the correlation tensor T^^^ of the action of a SLOCC operator or a local 
SL{N, C) group |191j on the state |177j . If this programme is successful, we 
hope to classify the multipartite entanglement into classes that are SLOCC 
inequivalent. All this may have applications in thermal entanglement of 
many particles systems [ jl92] . Finally, it will be very interesting to seek 
applications of our entanglement measure to situations like quantum phase 
transitions fl93j . transfer of entanglement across spin chains fl94j etc. 
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